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Abstract. When the soil water balance method is applied 
at a field scale, estimation of the spatial variability and 
confidence interval of actual evapotranspiration is rare, 
although this method is sensitive to the spatial variability 
of the soil, and thus to the sampling strategy. This work 
evaluated the effect of soil sampling strategies for soil 
water content and water flux at the bottom of the soil 
profile on the estimation of the daily and cumulative 
evapotranspirations. To do that, according to the statisti- 
cal properties of daily measurements in a field experiment 
with a soybean crop, the water content and flux through 
the base to the soil profile in space (field scale) and time 
(daily scale) were simulated. Four different sampling 
strategies were then compared, and their effects on daily 
and seasonal cumulative evapotranspirations quantified. 
Strategy 1 used ten theoretical sites randomly located in 
the field. The daily water content estimates were assumed 
to be available each day from these same ten locations, 
which were located from 0.15 m to 1.55 m in depth, with 
space steps of 0.10 m. Strategy 2 assumed that daily wa- 
ter content estimates combined two sources: in the 0.00— 
0.20 m soil layer, ten theoretical sites were selected but 
changed every day, with thin soil layers for soil moisture 
sampling, from 1 to 5 cm in thickness. In the 0.20—1.60 m 
soil layer, the daily water content estimates were assumed 
to come from the same ten locations (the first soil mois- 
ture estimate was located at 0.25 m, and the others were 
located every 0.10 m until 1.55 m). Strategy 3 used ten 
theoretical sites located in the field, as in strategy 1, how- 
ever the water content estimates in the 0.00—0.20-m soil 
layer were assumed to come from accurate water content 
measurements (soil layers from 1 to 5 cm in thickness), 
while for the 0.20—1.60 m soil layer, the strategy was 
similar to strategies 1 and 2. Strategy 4 used 10 new theo- 
retical locations of measurement every day. Precise water 
content estimates for thin layers were assumed to be 
available in the 0.00—0.20 m soil layer as in strategy 2. 
The layers for water content estimates in the 0.20—1.60 m 
were similar to those of strategies 1, 2, and 3. Results 
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showed that the spatial variability of the daily actual 
evapotranspiration may not be negligible, and differences 
from approximately +1.0 mmd~! to +3.0 mmd~! were 
calculated between the four sampling strategies. Strate- 
gy 1 gave the worst results, because variations in the wa- 
ter content of the top soil layers were neglected, and thus 
the daily evapotranspiration was underestimated. Strate- 
gy 2 led to a considerable variability for estimating daily 
evapotranspiration which was explained by the effect of 
the spatial variability due to the daily site sampling for 
the top soil layers (0 to 0.2 m). Strategy 3 appeared to be 
the best practical compromise between practical field 
considerations and the necessity to obtain accurate eva- 
potranspiration measurements. The accuracy of daily 
evapotranspiration could reach +0.5 mmd~', and could 
be further improved by increasing the number of mea- 
surement sites. The best results were obtained with strat- 
egy 4, although such a destructive and time-consuming 
strategy is not likely to be practical. 





The estimation of evapotranspiration under natural con- 
ditions at different space scales (local, field, regional), 
and time steps (instantaneous, daily, seasonal) is crucial 
for water conservation and management, in particular 
planning the efficiency of irrigation and plant water use 
(Stone et al. 1973; Molz 1981), or field validation of water 
transport models (Camillo etal. 1983; Witono and 
Bruckler 1989). Available methods for measuring evapo- 
transpiration are based on soil water balance estimation 
(Daudet and Vachaud 1977; Vachaud et al. 1978), direct 
measurements using lysimeters, micrometeorological 
methods, including portable chambers, eddy correlation 
or Bowen ratio (Perrier et al. 1975; Reicosky and Peters 
1977; Itier 1981; Pal Arya 1988), or plant physiological 
methods (Sakuratani 1981; Valancogne and Nasr 1989). 
But none of these methods of measuring evapotranspira- 
tion is applicable under all conditions (Sharma 1985). 
The plant physiological methods based upon sap flow 
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measurements relate to a single plant and only measure 
transpiration, the soil water balance measurements 
provide local estimates of evapotranspiration, whereas 
the micrometeorological methods refer to field or region- 
al scales for continuous and homogeneous crops. 

When applied at a field scale, the major drawback of 
the soil water balance method is that it lacks sensitivity to 
measure daily evapotranspiration values, and that it is 
sensitive to the spatial variability of the soil, and thus to 
the sampling strategy (Peyremorte et al. 1972; Dunin and 
Aston 1981; Rambal et al. 1984; Vauclin et al. 1984). Nu- 
merous experimental data on the spatial variability of the 
soil water content have been published (Nielsen et al. 
1973; Gajem et al. 1981). However, estimation of the spa- 
tial variability of actual evapotranspiration are rare, al- 
though such data are crucial for the development of 
stochastic water transport models. This work evaluated 
the effect of soil sampling strategies and measurement 
methods (neutron probe measurements and/or destruc- 
tive mass based soil samples) on the estimation of the 
daily and cumulative evapotranspiration, and to quantify 
the spatial variability of evapotranspiration at field scale 
when soil water balance methods are used. To do this, 
water content and deep flux variations in space (field 
scale) and time (daily scale) were simulated, according to 
observed statistical properties of daily water content and 
the water flux at the bottom of the soil profile measure- 
ments of a field experiment. Different sampling strategies 
were then simulated and their effects on daily and season- 
al cumulative evapotranspiration were quantified. 


Simulation procedure for the soil water balance variability 
Objectives and procedure 


The purpose of the proposed procedure was to simulate 
a realistic spatial variability of daily water contents and 
deep fluxes (at 1.55 m), for estimating the evapotranspi- 
ration at field scale and its accuracy according to differ- 
ent sampling strategies. Every day, the distributions of 
the daily water content profiles and the flux at 1.55 m 
(qZmax) Were generated at 1000 theoretical sites according 
to a statistical procedure. This procedure assumed that: 


i) Experimental mean and variance of water contents 
and water fluxes were unbiased estimates of the true 


mean and variance at the field scale. 


ii) The variability of data was isotropic and spatially in- 
dependent, and could be described by a normal distribu- 
tion. The variability of soil data may be such that the 
isotropic and random assumption are not valid, and the 
magnitude of the error with the assumption of random- 
ness as opposed to some trends in the variation that is 
related to soil properties should be separated. Neverthe- 
less, in our case, data collected on the same site (Bertuzzi 
et al. 1987) showed that at all depths in the 0.00—0.20 m 
soil layer, the histogram of dry bulk density was normally 
distributed and the experimental semivariogram did not 
show any trend (pure nugget effect). Moreover, results 
with gravimetric water contents measured in the 0.00- 


0.20 m soil layer were similar, and purely random distri- 
butions of water contents were observed. 


For the water content at each depth, and the deep flux 
(qZmax)» the simulation procedure started with the sam- 
pling of 1000 probabilities [p()] distributed on the unit 
normal distribution. For each day, the mean (1) and the 
standard deviation (s) were estimated from field data, 
and it was then possible to simulate the distribution of 
1000 spatially distributed values [x(i)] using the relation- 
ship: 
x(i)=p(i)-s+p. (1) 


Moreover, for the water content profiles, we observed 
from experimental data that the drier series stayed drier 
and the wetter sites stayed wetter, i.e. the order of the sites 
(more or less dry or wet) was approximately maintained 
versus time (see below in the next section). This was taken 
into account in the statistical procedure, in that the same 
probability p(i) generated by the statistical computation 
was used every day at a given site to compute the distribu- 
tion of water content whatever the depth. We applied the 
same procedure for the deep flux (qz,,,,)- This maintained 
the physical coherence observed on daily field data, i.e. 
the drier sites at the start of the procedure stayed drier 
and the wetter sites stayed wetter. 

The soil water balance could then theoretically be cal- 
culated at each theoretical site. Thus, for the 3000 m? 
field, it was assumed that each was representative of the 
soil water balance of a 3 m? soil surface. Nevertheless, 
such theoretical work could be applied whatever the 
space scale and time step. The mean water balance com- 
ponents calculated on 1000 points were considered as a 
‘reference’, i.e., the ‘true’ value of the water balance com- 
ponents. Data sets coming from different sampling 
strategies in the simulation were then sampled and com- 
pared to this theoretical ‘reference’, to finally compare 
their relative performances. 


Comparing four sampling strategies 


The characteristics of the four sampling strategies are 
indicated below 


i) Strategy 1 used ten theoretical sites randomly located 
in the field. The daily water content estimates, and deep 
water fluxes, were assumed to come every day from these 
same ten location (neutron probe measurements for ex- 
ample). The first estimate began at 0.15 m, and the others 
were located every 0.10 m until 1.55 m (0.25; 0.35; 0.45; 
0.55; 0.65; 0.75; 0.85; 0.95; 1.05; 1.15; 1.25; 1.35; 1.45; 
1.55 m). 


ii) Strategy 2 assumed that daily water content estimates 
combined data of two origins: in the 0.00—0.20 m soil 
layer, ten theoretical sites were selected for soil moisture 
estimates, but these theoretical sites changed every day 
(destructive soil samples for example). In the 0.20-1.60 m 
soil layer, 10 theoretical locations were also selected, but 
the daily water content estimates and deep water fluxes, 
were assumed to come every day from these same ten 
locations (neutron probe measurements for example). 





The water content estimates in the 0.00—0.20 m layer 
used soil moisture of nine soil layers: 0.00—0.01, 
0.01—0.02,, 0.02—0.03, 0.03—0.04, 0.04—0.05, 0.05—0.07, 
0.07—0.10, 0.10—0.15, 0.15—0.20 m. For the 0.20—1.60 m 
soil layer the first soil moisture estimated was located at 
0.25 m, and the other were located every 0.10 m until 
1.55 m. 


iil) Strategy 3 used ten theoretical sites located in the 
field, as strategy 1. Water content estimates in the 0.00-— 
0.20 m soil layer were assumed to come from accurate 
and non-destructive water content measurements (using 
capacitive probes or Time Domain Reflectrometry mea- 
surements for example). Thus, the estimation of water 
storage in the 0.00—0.20 m soil layer was from soil mois- 
ture estimates of nine soil layers: 0.00—0.01, 0.01—0.02, 
0.02—0.03, 0.03—0.04, 0.04—0.05, 0.05—0.07, 0.07—0.10, 
0.10—0.15, 0.15—0.20 m. For the 0.20—1.60 m soil layer, 
the first soil moisture estimate was located at 0.25 m, and 
the other were located every 0.10 m until 1.55 m (similar 
to those of strategies 1 and 2). 


iv) Strategy 4 used 10 new theoretical locations of mea- 
surement evey day (destructive soil samples from the soil 
surface to 1.55 m for example). Precise water content es- 
timates for thin layers were assumed to be available in the 
0.00—0.20 m soil layer as in strategy 2. The layers for 
water content estimates in the 0.20—1.60 m were similar 
to those of strategies 1, 2, 3, i.e. every 0.1 m, starting at 
0.25 m. 


For these four strategies, the sites were sequentially 
selected randomly and the procedure prevented the same 
site being selected more than once. 


Materials and methods 
Test site 


The field experiment was conducted at the Soil Science Laboratory 
of the Institut National de la Recherche Agronomique (INRA) at 
Montfavet (France) from July to the beginning of October 1989. 
Fifty-two soil samples (0O—0.2 m depth) were taken on a hexagonal 
grid (0.3 ha) to evaluate the spatial variability of the particle size 
distribution. The soil is a clay loam with about 27% clay and 62% 
silt and the standard deviation for particle size classes is between 0.9 
and 1.2% (Table 1), and thus the top soil layers of the field were 
regarded as ‘homogeneous’. Clay contents were approximately con- 
stant until 0.6 m, and then increased (Table 2). No quantitative data 
about spatial variability of particle size distribution below 0.2 km 
were available, but field observations showed that spatial variability 
was low until 0.8 m in depth and increased after 0.8 m. The exper- 
imental field was subdivided into three subplots (Fig. 1). Soybean 
(Glycine max, (L.), cv “Weber”’) was sown on 3 July in rows 0.30 m 
apart on each subplot, and harvested on 4 October. Data were 
collected from Day 16 until Day 93 after sowing. The three subplots 
were irrigated with a sprinkle irrigation equipment (approximately 
7mm h~') on 23 (77.5mm) and 30 (57.5 mm) July, and on 12 
(45.5mm) and 31 (17mm) August. It rained on 14 August 
(2.6 mm), and on 8 (1 mm), 10 (35 mm), 13 (1.6 mm), and 18 (1 mm) 
September. The total water supply (rainfall and irrigation) was 
238.7 mm, of which 17.3% was accounted for by rainfall. 


Soil and crop measurements 


Measurements combining the use of neutron probe and tensiome- 
ters were made at four sites on subplot B (Fig. 1). Each site was 


107 


Table 1. Basic statistics of particle size distribution (clay, fine silt, 
coarse silt, sand, and organic matter) for the 0.00—0.20 m soil layer, 


over a 0.3 ha field (52 locations) 





Clay 


(%) 


Fine 
silt 
(%) 


Coarse 
silt 
(%) 


Sand 


(%) 


OM 





Minimum 
Maximum 
Mean 

SD 

CN 


25.1 
29.6 
2h 
41 
4.0 


43.8 
48.3 
45.9 
$1 
2.6 


13.8 
17.8 
15.8 
0.9 
5.6 


9.0 
13.1 
11.0 

0.9 

8.3 





OM: organic matter 


SD: 


standard deviation 


CV: coefficient of variation (%) 


Table 2. Mean particle size distribution (clay, fine sild, coarse silt, 
sand) for the 0.00—1.6 m soil layer, over the 0.3 ha field 





Depth (m) 


Clay 


(%) 


Fine 
silt 
(%) 


Coarse 
silt 
(%) 


Sand 


(%) 





Pj i. 
24.3 
25:5 
28.9 
36.8 
30.5 
33.8 
Bo 


45.9 
46.1 
43.4 
34.5 
30.7 
24.2 
23.0 
30.0 


15.8 
17.4 
17.7 
IPA 
9.7 
8.3 
72 
10.7 


11.0 
12.2 
13.4 
21.3 
22.8 
40.0 
36.0 
23.8 





SUBPLOT A 


SUBPLOT B 


SUBPLOT C 





=35m 


length 





FETE HEHEHE HEH HHH 
FEEFE HEHE PEPE HH 
FEEEE HEHE HEHE ET 
FEE HEHE HEE E THT 
FEET HEFE TEETH ET 
FETE HEHEHE HE+ EHH 
FEREE HEHE EPE TT 
FETE EET HEHEHE +H FHF 
FETE HE tHE HHH HHH 
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length = 60 m 


© Neutron measurement 


+ Gravimetric measurement 


N Neutron acces tube 
T Tensiometers 


Fig. 1. The experimental setup 
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equipped with an aluminium access tube (length 2.0 m) surrounded 
with 15 tensiometers (mercury manometers) which were arranged in 
a circle of 1 m diameter around the access tube. The neutron probe 
(model SOLO 40) an the tensiometers (model DTM 5000, 15 mm in 


external diameter, 900 mb in maximum bubbling pressure) were 
manufactured by NARDEUX (Saint Avertin, France). Tensiometer 
measurements were made at 0.5 m, 0.10 m and every 0.10 m from 
0.15 to 1.55 m below the soil surface. Neutron probe measurements 
were made at the same depths except the first two depths (0.05 m, 
0.10 m). Measurements were made between 10.00h and 12.00h 
GMT every day until Day 75 after sowing, then at increasing time 








interals from two to five days until harvest. Near the neutron sites, 
a gamma-ray transmission probe was used for bulk density mea- 
surements of the 0.00—0.20 m soil layer (Bertuzzi et al. 1987), and 
a retrodiffusion gamma ray probe manufactured by NARDEUX 
(Saint Avertin, France) for the 0.20—1.60 m soil layer, using the 
same access tubes as for the neutron probe measurements. Dry bulk 
densities varied approximately from 1200 kg m~? to 1700 kg m=? 
(Fig. 2) between the soil surface and the subsoil. The spatial vari- 
ability between the four locations was rather high, especially in the 
top [0—0.25 m] soil layers and this may be due to previous wheel 
track compaction, as well as the tillage effects. Field observations 
showed that the depth of tillage could vary from 0.12 m to 0.22 m, 
and therefore the variability in the top soil layers was partly due to 
variations in the thickness of the tillage layer. The variability was 
also high at the bottom of the soil profile (1.30 to 1.60 m), and soil 
sampling with auger showed that it was probably due to variations 
in soil texture. 

The neutron probe was calibrated in the field at the beginning 
of the experiment and during the observation period (347 data pairs 
of count rates and volumetric water contents were collected). Cali- 
bration was made using the corresponding dry bulk densitiy profile 
measured at each site to compute the volumetric water content. 
“ : Different calibration relationships were grouped into account for 
a ’ ’ 7 r * ’ three soil layers, 0.15 m, 0.25—0.85 m, and 0.95—1.55 m (Table 3). 

4100 1200 1300 1400 1500 1600 1700 1800 Thus, the residual standard deviation varied from 0.016 m? m3 to 
Dry bulk density (Kg/m3) 0.018 m* m~?. Concurrently with neutron probe and tensiometer 

measurements in subplot B, gravimetric cores using an auger was 

Fig. 2. Experimental dry bulk density profiles (sites 1 to 4) taken on subplots A and C to characterize every day more precisely 
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Fig. 3. Mean values and standard deviations of the volumetric water content at three selected depths (0.005 m, 0.250 m, 1.050 m) 
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Fig. 4. Measured water storage at four sites versus time (Period 2: 


from day 28 to day 40 after sowing, Period 3: from day 41 to day 
59 after sowing) 


Table 3. Fitted parameters of the calibration of the neutron probe 
for 3 soil layers (0.00—0.15, 0.25—0.85, 0.95—1.55 m) and for the 
whole soil profile (0.00—1.55 m) 





Depth n a 
(m) 


0.00—0.15 32 
0.25—0.85 194 
0.95-1.55 


0.00—1.55 347 





0.00075 
0.00085 
0.00087 


0.00076 





number of data 
slope of the regression line 
intercept of the regression line 
coefficient of correlation 

: standard deviation (m? m7?) 


the contribution of near-surface soil layers to water storage varia- 
tions. At eight randomly selected locations, nine core samples were 
taken at predefined depths (0.00—0.01, 0.01 —0.02 , 0.02—0.03, 0.03- 
0.04, 0.04—0.05, 0.05—0.07, 0.07—0.10, 0.10-0.15, and 0.15- 
0.20 m) to determine the gravimetric water content by oven-drying 
for 24 hat 105°C. The corresponding volumetric water content was 
calculated by using the mean dry bulk density profile (from 0 to 
0.20 m). Root distribution and depth were evaluated using the pro- 
file wall method (Tardieu, 1988). One to four vertical soil profiles 
(0.30 m in width) were exposed mechanically and soil carefully re- 
moved every week or more to expose roots visible at the wall of the 
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profile. Roots present were recorded with a 0.02 m square grid. 
Observations were always made at a depth between the soil surface 
and the maximum visible root depth. 

Figure 3 shows the variation of the volumetric water content 
versus time for three selected depths (0.005 m, 0.250 m, 1.050 m). 
The first depth (0.005 m) corresponds to destructive soil samples (8 
samples per day), whereas the last two depths (0.250 m and 1.050 m) 
correspond to water contents derived from the neutron probe mea- 
surements. The order of magnitude of the standard deviation was 
approximately between 0.01 and 0.06 m? m7 °. Standard deviations 
were greater for the top soil layer (0.005 m), where water supply 
and/or evaporation processes generally induced a high spatial vari- 
ability. Moreover, the variability of water content with time could 
not be considered constant, but as dependent on water supply 
events (rain or irrigation). For deeper layers, 1.050 m for example, 
water content variability was low until approximately 55 days after 
sowing, and increased after 55-72 days after sowing as root growth, 
which was irregular in the field, induced a more or less rapid water 
depletion in the deeper soil layers. 

Figure 4 showed the measured water storages at the four sites 
versus time for period 2: from day 28 to day 40 after sowing: and 
period 3: from day 41 to day 59 after sowing. Two periods were 
selected in Fig. 4, but similar results were obtained for all the exper- 
imental periods. Figure 4 indicated that the drier sites stayed drier 
and the wetter sites stayed wetter, i.e. the order of the sites (more or 
less dry or wet) was approximately maintained with time. As indi- 
cated, this experimental fact was taken into account in the statistical 
procedure, in that the same probability p(i) generated by the statis- 
tical computation was used every day at a given site to compute the 
distribution of water content whatever the depth, and this main- 
tained the physical coherence because the drier sites at the start of 
the procedure stayed drier and the wetter sites stayed wetter. Similar 
results obtained under trickle irrigation were presented by Rolston 
et al. (1991). 


Water balance estimates 


When the soil tensiometric profiles showed a zero flux plane z, 
below the maximum root depth, evapotranspiration was calculated 
from the water storage variations between the soil surface and the 
zero flux plane z, (Vachaud et al. 1978). When zero flux plane was 
no longer visible on the soil tensiometric profiles, evapotranspira- 
tion was calculated from the water storage variations between the 
soil surface and a soil depth (z,,,,=1.55 m) below the maximum 
root depth, and from the estimated water flux at this given depth 
(Zmax)> depending on water potential gradients and unsaturated hy- 
draulic conductivity. Due to the development of the root system 
(Fig. 5), the zero flux plane was always located under root depth, 
from the beginning of the experiments to days 71, 60, 59, and 71 
after sowing for sites 1 to 4, respectively. After these dates, the 
position of the zero flux plane for each site was unknown, because 
it was deepeer than the deepest tensiometer, maximum depth 
1.55 m. 

The unsaturated hydraulic conductivity — volumetric water con- 
tent relationship in the bottom of the soil profile (1.55 m) was given 
by Van Genuchten (1980) according to Mualem’s model (Mualem 
1976): 


K (0)=K,,, : [9]"? - [1—[1-[0]""I"? (2) 
with 

0 =[0—6,]/[0, —9,] (dimensionless) 

K,,,: Saturated hydraulic conductivity (m s~*) 

0,: Saturated volumetric water content (m* m7 *) 


6.: Residual volumetric water contents (m* m~ *) 
m: Fitting parameter, 


where the terms 0,, 6, and m come from the water potential — 
volumetric water content relationship (Van Genuchten, 1980) given 
by: 


w(0)=1/a - [[@)"/"—1)""" (3) 
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with 
=1-1/n 
n, a: Fitting parameters. 


The parameters 6,, 0,, n, a were obtained by fitting Eq. (3) using 
water potential-volumetric water content data at depth z,,,, 
(Table 4). The procedure for estimating the unsaturated hydraulic 
conductivity-volumetric water content relationship then used ten- 
siometric and neutron probe water content data. For the days when 
a zero flux plane Z, existed, the unsaturated hydraulic conductivity- 
water content relationship [K (@)] was estimated using a non-linear 
procedure (Bard 1974) by optimizing the term K,,, of the K(@) curve 
(Eq. (2)) to verify the equation [4], given as follow: 


[S(Z.z —S (Zo, Zmax)jl/At;, j+1 
= —K(6) Liane 5 [(dy//dz) Zmax 1] (4) 


av 


with 


S(Zo,Zmax)j+1: Water storage from Z, to Z,, 

S(Zo, Zmax);: Water storage from Zp to Z,,,, at time j(m) 

At, 541. Time step from time j to time j+1(s) 

K(@): ° unsaturated hydraulic conductivity (m s~') 
: water potential (m) 

: depth (m). 


ax at time j+1(m) 


The terms S(Z 9, Zmax)j+1 ANd S (Zp, Zmax); Were derived from neu- 
tron probe measurements, and the water potential gradients (dy/ 
6Z)Zmax {rom tensiometric data. Thus, only the term K,,, was to be 


estimated when the water balance between z,. and Zn... had been 
established (Table 4), and once the K (@) relationship was identified, 
the daily evapotranspiration could be estimated when the zero flux 
plane was no longer available. The functional relationships that 
express conductivity as a function of water content and the water 
potential as a function of water content do not account for hys- 
teresis. Moreover, based on the variation of the bulk density as 
measured with the gamma probe, there should be significant differ- 
ences in hydrodynamic properties between the four experimental 
locations. Thus, Table 4 indicates that different K,,, values were 
estimated for each location. The spatial variability for the parame- 
ter K,,, between the four locations was high (mean of 2.31 1077 m 
s~', standard deviation of 2.45 10-7 ms~'). The calculated water 
fluxes at the bottom of the root zone, and daily evapotranspi- 


sites) 


Table 4. Fitted parameters for the water potential-water content 
and for the hydraulic conductivity-water content relationships at 
four sites 





Site 1 2 





Water content — water potential relationship 


6, (m3 m7 3) 0.392 0.418 0.341 
6, (m3 m~ 4) 0.231 0.063 0.142 
a (m~') 0.727 1.013 0.825 
n i es 1.192 1.143 


Water content — hydraulic conductivity relationship 


Ka(90-mrs: *) 0.39 5.39 3.19 
Mean (1077 ms~') 2.31 
Standard deviation 2.45 
(1077 ms~*) 





ration for the four sites showed great variations between days 
(Fig. 6). Positive fluxes at 1.55 m corresponded to downward fluxes 
(drainage), whereas negative fluxes corresponded to upward flux 
(capillary raise). Spatial variability on daily evapotranspiration was 
essentially due to the spatial variability of the soil water content in 
the soil layers, and to measurement errors coming from neutron 
probe measurements. Residual standard deviation of the calibration 
relationships could have a not negligible effect on evapotranspira- 
tion estimates, when standard deviation is expressed as a fraction of 
daily evapotranspiration. For example, an error of 0.010 m* m=? 
for a soil layer of 0.1 m depth corresponds to an error of 1 mm in 
evapotranspiration. Thus, errors from soil moisture measurement 
may result in the estimated actual evapotranspiration being greater 
than the potential evaporation. On the other hand, calculations 
showed that accuracy of weighting in the gravimetric method in- 
duced relative errors in gravimetric water content of less than 0.2%. 
Thus, errors in gravimetric measurements were considered as negli- 
gible. As a matter of fact, water storage calculations for the 0.00- 
0.20 m soil layer between two consecutive days varied greatly from 
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one day to another, depending on soil moisture sampling. Water 
fluxes at 1.55 m depth showed great variability between sites for 
both downward and upward fluxes (Fig. 6), and the mean flux 
varied approximately from +1 to —1 mm d~! during the experi- 


ment. 
Results and discussion 
Comparing the four strategies 


An example of the estimated variability of simulated wa- 
ter contents, water fluxes at 1.55 m and daily evapotran- 


Water content (m3.m-3) - Layer 0.00 - 0.005 m 
Mean=0.32 Std=0.01 Min=0.28 Max=0.36 


Deep water flux (mm.d-1) 
Mean=0.32 Std=0.17 Min=-0.08 Max=1.06 


Evapotranspiration (mm.d-1) 
Mean=8.74 Std=0.40 Min=7.63 Max=9.84 


transpiration versus time 


spiration at the field scale is shown in Fig. 7. Here spatial 
variability in water contents (0.28 m? m~? to 0.36 m? 
m~ 3), and water fluxes at the bottom of the soil profile 
(—0.08 to 1.06 mm d~! at 1.55 m) resulted in a variabil- 
ity in the evapotranspiration (from 7.63 to 9.84 mm d™') 
which was not negligible. In Fig. 8 and 9, daily and cumu- 
lative evapotranspiration values according to strategies 1 
to 4 compared to the reference daily and cumulative 
evapotranspiration. Daily evapotranspirations according 
to the four strategies are statistically compared to refer- 
ence daily evapotranspirations in Fig. 10, where the solid 
line represents a slope do unity, and the dotted line the 
regression line fitted to the data: 
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Fig. 8. Comparison between the mean simulated daily evapotranspiration (reference from 1000 sites) and the estimated daily evapotranspi- 


ration using four different sampling strategies 


— Strategy 1 gave the worst results, because variations in 
the water content of the top soil layers were neglected, 
and thus the daily evapotranspiration was underestimat- 
ed (Fig. 8). For example, evaporation after moderate 
rewetting of the top soil layers cannot be taken into ac- 
count using this strategy. 


— Strategy 2 gave a quite high variability for estimating 
daily evapotranspiration and gave less satisfactory re- 
sults than strategy 3. This was explained by the effect of 
the spatial variability due to the daily site sampling for 
the top soil layers (0 to 0.2 m) resulting from the gravi- 
metric method. 


— Strategy 3 appeared to be the best practical compro- 
mise between restrictions in the field and the necessity to 
obtain accurate evapotranspiration measurements. Nev- 
ertheless, it would require the use of automatic sensors 
for soil moisture measurements in the top soil layer which 
have a fine vertical resolution such as capacitive probes. 
The accuracy of daily evapotranspiration could reach 
+1.0mm d™', and this could be increased further by 
increasing in the number of measurement sites. 


— The best results were obtained with strategy 4. In terms 
of the theory, this strategy simply means increasing the 
sample population, which in turn should increase the 
accuracy of estimation. Nevertheless, such a strategy is 
not likely in practical cases: Given the methods of field 


water content measurement now available this strategy 
will be very destructive and time-consuming. 


Influence of the number of sites 


In the case of strategy 3, the effect of the number of sites 
on the mean daily or cumulative evapotranspiration is 
shown on Fig. 11. The daily evapotranspiration was de- 
termined at field scale with an increasing number of sites: 
3, 5, 20, 50, 100, 250, 500, 1000 where the soil water 
balance was calculated every day. Sites were selected ran- 
domly and added to the sites selected at the previous step. 
Ratios between simulated mean evapotranspirations and 
the corresponding reference values were computed every 
day for daily and cumulative evapotranspiration. The 
challenge here is determining the minimum number of 
samples that will give a reasonable answer from the prac- 
tical standpoint, as increasing the number of values in the 
average will always improve the average. Under our field 
conditions and assuming that the range 0.90—1.10 for the 
ratio between simulated mean evapotranspirations and 
the corresponding reference values was satisfactory, ac- 
curate estimates of daily evapotranspiration at field scale 
required a minimum number of about 30 sites. This num- 
ber was smaller, about 10 sites, for cumulative evapotran- 
spiration, which is less sensitive to spatial variability. 
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Fig. 9. Comparison between the mean simulated cumulative evapotranspiration (reference from 1000 sites) and the estimated cumulative 
evapotranspiration using four different sampling strategies 


STRATEGY 1 STRATEGY 2 








10 12 14 
ET reference (mm.d-1) 
10 12 14 


ET reference (mm.d-1) 
Oo 2 4° 6 & 


024 6 8 




















0 2 4 6 8 10 12 14 0 2 4 6 8 10 12 14 
ET (mm.d-1) ET (mm.d-1) 


STRATEGY 3 STRATEGY 4 








10 12 14 
10 12 14 


r =0.988 


ET reference (mm.d-1) 


ET reference (mm.d-1) 
02 4 6 8 


0 24 6 8 


Fig. 10. Comparison between four sampling 
strategies to estimate the simulated daily evapo- 
transpiration 
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Fig. 11. Influence of the number of sites on the estimation of the simulated daily and cumulative evapotranspirations using strategy 3 


These results confirm the fact that the determination at 
field scale of the daily actual evapotranspiration can be 
very inaccurate in the case where a limited number of 
measurement sites is used. Cumulative evapotranspira- 
tion estimates were less sensitive to the spatial variability 
of soil properties than the daily estimates. 


Conclusions 


A procedure for data simulation was established and per- 
formed to study, at a field scale, the effects of different 
sampling strategies on the actual daily and cumulative 
evapotranspiration. On the basis of realistic simulated 
data, this study shows that the spatial variability of the 
daily actual evapotranspiration is not negligible, even in 
small fields, 0.3 ha. Comparison of four soil moisture 
sampling strategies shows that a daily random sampling 
of sites provides the most accurate and unbiased daily 
evapotranspirations. Strategy that uses automatic sen- 
sors having a fine vertical resolution for soil moisture 
measurements in the top soil layers, such as capacitive 
probes, appeared to be the best practical compromise 
between field restrictions and the necessity to obtain ac- 
curate evapotranspiration measurements. In such a case, 
the accuracy of daily evapotranspiration could reach 
+1.0 mm d~', and could be further improved by increas- 
ing the number of measurement sites. These results were 


obtained assuming random spatial variability of the soil 
characteristics and it is necessary to extend them in the 
case of non-random spatial variability. The simulation 
procedure is able to provide a random field distribution 
of actual evapotranspirations. Simulation techniques, 
such as the turning band method, can be used to generate 
spatially correlated fields of evapotranspiration and then 
to compare the influences of various sampling strategies 
on the estimates of evapotranspiration at a field scale. 
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Abstract. Growth and yield responses of developing 
almond trees (Prunus amygdalus, Ruby cultivar) to a 
range of trickle irrigation amounts were determined in 
1985 through 1987 (the fifth through seventh year after 
planting) at the University of California’s West Side Field 
Station in the semi-arid San Joaquin Valley. The treat- 
ments consisted of six levels of irrigation, ranging from 50 
through 175% of the estimated crop evapotranspiration 
(ET.), applied to a clean-cultivated orchard using a line 
source trickle irrigation system with 6 emitters per tree. 
ET, was estimated as grass reference evapotranspiration 
(ET,) times a crop coefficient with adjustments based 
upon shaded area of trees and period during the growing 
season. Differential irrigation experiments prior to 1984 
on the trees used in this study significantly influenced the 
initial trunk cross-section area and canopy size in the 
50% ET, treatment and 125% ET, treatment. In these 
cases, treatment effects must be identified as relative ef- 
fects rather than absolute. The soil of the experimental 
field was a Panoche clay loam (nonacid, thermic, Typic 
Torriorthents). The mean increase in trunk cross-section- 
al area for the 3-year period was a positive linear function 
(r? = 0.98) of total amounts of applied water. With in- 
creases in water application above the 50% ET, treat- 
ment, nut retention with respect to flower and fertile nut 
counts after flowering, was increased approximately 
10%. In 1985 and 1987, the nut meat yields and mean 
kernel weights increased significantly with increasing wa- 
ter application from 50% to 150% ET.. Particularly in 
the higher water application treatments, crop consump- 
tive use was difficult to quantify due to uncertainty in 
estimates of deep percolation and soil water uptake. 
Maintenance of leaf water potentials higher than 
—2.3MPa during early nut development (March 
through May) and greater than —2.5 MPa the remainder 
of the irrigation season (through August) were positively 
correlated with sustained higher vegetative growth rates 
and higher nut yields. 


Correspondence to: R. B. Hutmacher 


Trickle irrigation of almond (Prunus amygdalus) orchards 
in California’s San Joaquin Valley has increased signifi- 
cantly since the mid-1970’s partly because of expected 
economic benefits related to improved water and nutrient 
use efficiency (Micke and Kester 1978; Fereres et al. 1981; 
Krauter 1985). However, the potential water conserva- 
tion benefits have not generally been realized because 
research and field experiences have not developed the 
appropriate relationships between the trickle irrigation 
system, its operation, the soil hydraulic properties, and 
the irrigation requirement of almond. To further compli- 
cate concepts of appropriate trickle irrigation manage- 
ment for almonds, some of this expansion in planted 
areas has occurred in regions with heavier clay and clay 
loam soils in the western San Joaquin Valley rather than 
in the coarser-textured soils typical of much of the al- 
mond production area in the eastern San Joaquin Valley. 

Spatial variation in soil water distribution due to dif- 
ferences in flow geometry and surface soil infiltration 
characteristics, combined with uncertainty in defining ac- 
tive root zones to monitor for water uptake, makes soil 
moisture-based irrigation scheduling difficult for many 
tree crops (Jury and Earl 1977; Fereres et al. 1981). In 
previous reports on the project described in this paper, 
Rolston et al. (1991) and Nightingale et al. (1991) dis- 
cussed spatial variability of soil water distribution with 
respect to distance from the tree and trickle lateral. They 
described the very complicated wetting pattern and con- 
siderable potential for deep percolation losses and leach- 
ing with the line-source trickle system used in this exper- 
iment. 

Increasing availability of weather stations and evapo- 
ration pans for meteorological-based evapotranspira- 
tion (ET) predictions may provide the basis for greater 
use of reference ET-based irrigation scheduling, provided 
accurate information is available to relate water use by 
developing crops to standard measures of reference ET 
(Aljibury et al. 1974; Doorenbos and Pruitt 1977; Fereres 
et al. 1981). Trickle irrigation systems are particularly well- 
suited to automated operation to apply some percentage 
of reference ET at any desired irrigation frequency. 
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Under advective conditions of an immature orchard, 
Fereres et al. (1981) demonstrated a nonlinear relation- 
ship between percent groundcover and percent of ET of 
a mature almond orchard (assumed to have 70% ground- 
cover), while Aljibury et al. (1974) assumed a linear rela- 
tionship. Fereres et al. (1981) concluded that the nonlin- 
earity of this relationship was further accentuated by re- 
duced soil evaporation and sensible heat transfer associ- 
ated with bare, dry soil between tree rows under trickle 
irrigation, and developed a crop coefficient based upon 
percent ground cover for use with developing almond 
orchards. 

As with many other species, almond tree responses to 
water deficits depend on the rate of water stress develop- 
ment, method of irrigation, and available soil water as 
influenced by root distribution, wetted soil volume, and 
soil water holding capacity (Castel and Fereres 1982; 
Ruiz-Sanchez et al. 1988). Almond studies under trickle 
irrigation have concentrated on tree responses to long 
periods (weeks, months) of complete cessation of irriga- 
tion (El-Sharkawi and El-Monayeri 1976; Castel and 
Fereres 1982; Krauter 1985). Tree responses to continu- 
ing, high frequency irrigation in amounts ranging from 
deficit to full ET are the focus of this study, and little 
information is available on almonds under these condi- 
tions. 

This investigation reports the effects of the quantity of 
applied water, by a single-line source trickle irrigation 
system, on plant water status and selected growth and 
nut production parameters for almond grown in a clay 
loam soil. Tree responses are described in relation to ap- 
plied water plus rainfall rather than a complete water 
balance. This study has focused on the early nut-bearing 
years (4 to 7 years after planting) for the Ruby cultivar. 


Experimental procedures 


The experimental field was at the University of California’s West 
Side Field Station 9.6 km south of Five Points, California, on the 
west side of the San Joaquin Valley. The soil is a Panoche clay loam 
(nonacid thermic family of Typic Torriorthents) derived from re- 
cently deposited alluvium from calcareous sandstones and shales. 
Water retention properties within the tree root zone (upper 2 m of 
soil profile) are relatively uniform with depth (Nielsen et al. 1964), 
and the soil has a high plant-available water-holding capacity, aver- 
aging over 425 mm in a 2.5 m deep soil profile (Nielsen et al. 1973; 
Grimes et al. 1975). 

The experiment utilized almond (Prunus amygdalus) trees plant- 
ed in February 1980. The 3.4ha orchard had 63 rows of trees, 
14 trees per row. Tree spacing was a 6.1 by 6.1 m square pattern. 
Each treatment plot had 3 east-west-oriented rows of 7 trees each, 
2 border rows (Butte cultivar) and a central experimental row of the 
Ruby cultivar. One tree at each end of a row within each plot and 
the two rows of Butte cultivar served as the border trees to the five 
experimental trees. Both cultivars were grafted to Nemaguard root- 
stock which has shown moderate to good late-season drought toler- 
ance (Micke and Kester 1978; Micke and Kester 1988). The Ruby 
cultivar is a medium-size tree with moderate vigor, an upright 
growth habit, and a large number of spurs, and is classified as late 
blooming and late nut maturity relative to other California-adapted 
cultivars (Kester et al. 1985). 

No weeds were allowed to develop within the orchard, resulting 
in a clean orchard floor for the duration of the experiment. The area 
immediately surrounding the orchard on three sides was in row- 


Table 1. Monthly total reference evapotranspiration and accumu- 
lated heat units (base temperature 12°C) in March through August 
of consecutive years at West Side Field Station near Five Points, CA 





Month Accumulated heat 


units (base 12°C) 


Grass-reference evapo- 
transpiration (mm) 





Year Year 





1985 1986 1987 1985 1986 1987 





March 115 94 100 27" © TiS 70 
April 187 151 174 198 146 §8=214 
May 210 216 192 226 269 281 
June 219 234 = 214 379 361 346 
July 236 = 218 241 436 401 428 
August 210 207 228 370 386 = 437 





* Heat units calculated based on daily average temperature (deter- 
mined from daily maximum and minimum temperatures) minus a 
base temperature of 12°C 


crop production (cotton, processing tomatoes, cantaloupe) and in 
the downwind direction was bordered by table grapes. During the 
high ET periods of each year (June through September), surround- 
ing crops were not water-stressed and had leaf area indices ranging 
from a low of 2 in early June through 4 or more in July and August. 
Monthly grass reference ET and heat unit accumulation each year 
during the period from opening of new leaves in March through 
irrigation cutoff prior to nut harvest are shown in Table 1. The 
annual rainfall at the orchard for consecutive years was 136, 183, 
and 216 mm in 1985, 1986, and 1987. The long-term (15-year) aver- 
age annual rainfall averages 140 mm. 

A trickle irrigation experiment was initiated following the 1980 
planting and until the end of 1983 had 6 treatment combinations 
(T1 through T6) coded W,N,, W,N,, W3N,, W,N,, W,N,, and 
W,N,, respectively. The treatment designations W,, W,, and W, 
corresponded to trickle irrigation levels of 67%, 100%, and 133% 
of a calculated crop evapotranspiration (ET,), and N, and N, 
corresponded to either a single level of nitrogen (N) fertilizer ap- 
plied in the fall after nut harvest (subscript 1) or the same total 
amount of N applied continuously in small increments during the 
growing season (subscript 2). In the spring of 1984, the original 
experiment was discontinued and the number of water application 
treatments was increased to the 6 levels of trickle irrigation consid- 
ered in this paper (50, 75, 100, 125, 150, and 175% of the computed 
ET,, designated with the new treatment codes T1, T2, T3, T4, T5, 
and T6, respectively). The original 6 treatments became, respective- 
ly, the new T2, T3, T6, T1, T4, and TS treatments. The 1984 grow- 
ing season was considered transitional and only limited observa- 
tions from that year are reported. Tree trunk diameters were mea- 
sured in December of 1983, just prior to the transition to the new 
irrigation treatments, in order to assess carry-over effects of old 
treatments on initial tree size. The experimental design was a ran- 
domized complete block, with six irrigation treatments and seven 
blocks. 

The irrigation water was of relatively constant inorganic chem- 
ical composition, with a mean electrical conductivity (EC) of 
1.54 dS m~', and 3-year mean B, Na*, SO2~, and Cl~ concentra- 
tions of 0.08, 8.3, 5.3, and 3.3 mol 1~!. One drip irrigation line was 
used for each row of 7 trees, with the line placed parallel to the tree 
row direction within 10 cm of the tree trunks. There were 3 emitters 
on both sides of each tree at distances of 0.5, 1.25, and 2.0 m from 
the tree, with emitters delivering 4, 8, and 81h~', respectively. The 
full daily irrigation amounts were applied during either one or two 
continuous application periods. One m? of water per tree was equiv- 
alent to a depth of 26.9 mm of water if applied uniformly over the 
37.2 m? of land area allotted to each tree. 

The 6 levels of irrigation used in 1984 through 1987 were based 
on percentages of the calculated crop evapotranspiration (ET,) 





requirement derived from a nearby (100 m distance) weather station 
(data processed to determine grass reference evapotranspiration, 
ET,), and a crop factor based on the time of year and the percent 
of ground area shaded by the tree canopy (Micke and Kester 1978; 
Fereres et al. 1981). 

The shaded area (measured between 11:30 and 12:45 h, Pacific 
Standard Time) was expressed as a percent of the plot ground area 
allocated to each tree (37.2 m”) based on planting density. Shaded 
ground area was estimated according to the methods of Fereres 
et al. (1981) on three trees within each of four field replicates of a 
treatment by marking off the total 37.2 m? ground area per tree in 
sub-plots and classifying them as either shaded by tree leaves and 
branches or not shaded. Since treatment T3 was designated to re- 
ceive irrigation to replace 100% of estimated ET,, shaded ground 
area measurements of that treatment as a function of time were used 
in the calculation of ET,. Treatment-specific shaded ground areas 
were not used in the ET, calculations for other treatments. Irriga- 
tion volumes of other treatments (T1, T2, T4, T5, T6) were calculat- 
ed to apply 50, 75, 125, 150, and 175% of ET, calculated for treat- 
ment T3. 

The temporal persistence of soil-water contents and associated 
soil water potential as a function of location within the root zone 
was investigated as part of this experiment by Rolston et al. (1991). 
Changes in soil water content were determined by neutron attenua- 
tion from 12 to 25 neutron access tubes per tree distributed across 
a wide range of positions on one side of trees in all treatments other 
than TS, and soil water potentials at specific locations were moni- 
tored using tensiometers, as discussed by Rolston et al. (1991). Their 
data from this experiment was used to calculate a soil water balance 
for treatment T3. This data was used to modify crop coefficients in 
1986 and 1987. 

The irrigation season was divided into two-week periods. For 
each two-week period, the mean daily required quantity of water to 
meet 100% of estimated ET, was calculated. Post-harvest irriga- 
tions were applied in uniform amounts to all treatments during 
October and November of each year. No water was applied from 
December until the beginning of irrigation each spring. 

Trees received liquid nitrogen-fertilization (urea, 28% N by 
weight) through the trickle line at a rate of 0.34kg N tree~' just 
after harvest in 1985; 0.45 kg in 1986; 0.23 kg in May and 0.45 kg 
just after harvest in 1987. Early-summer (May) and mid-summer 
(July) leaf blade samples (50 leaves per plot) were analyzed for leaf 
blade N, P, K, Ca, Mg, Na, B, Cu, Mn, and Zn. The concentrations 
of these elements in all leaf samples were sufficient and not growth- 
or yield-limiting according to University of California guidelines 
(Micke and Kester 1978). 

The trees were pruned annually to encourage new growth for 
sustained production of fruiting spurs, following University of Cali- 
fornia pruning guidelines (Micke and Kester 1978). The effects of 
irrigation practices on tree growth and canopy size were evaluated 
by measuring: trunk circumference at the same point each year after 
leaf abscission; scaffold branch diameter at a distance ranging from 
20 to 30cm of attachment to the primary trunk; tree shade area 
from March through mid-July. Branch diameter measurements 
were made on two primary scaffold branches per tree, three trees per 
treatment, in each of three blocks during 1986 and 1987. 

Leaf water potentials (LWP) were determined with a pressure 
chamber (Soil Moisture Model 3000, Santa Barbara, CA) on 
fully-illuminated (photosynthetic photon flux density greater than 
1400 pmoles m~? s~') spur branch leaves collected from sites with- 
in 30 cm of branch tips. In preparation for measuring LWP, leaves 
were placed in plastic bags before excision. Water potentials were 
determined on four separate trees in three separate field blocks of 
each treatment. Leaf conductance to water vapor was determined 
on two fully-illuminated spur branch leaves on each of three trees 
in three field blocks of each treatment using a steady-state porom- 
eter (Li-Cor Model 1600, Lincoln, NE). 

Nut retention, as a percent of the number of flowers counted 
during peak bloom, was determined 5 to 7 weeks after the rate of 
embryo abortion had significantly declined and initial fruit drop 
had subsided, and again just before harvest in all three years. 
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These counts were made on marked 1 m long sections of branches, 
starting at the branch tip. Branches 1.5 to 2 m above the soil surface 
were measured in north, south, east and west quadrants of each of 
12 trees per treatment. 

Nut yields and nut characteristics for each irrigation treatment 
were based on 5 trees harvested in each of the 7 blocks per treat- 
ment. A modified commercial tree shaker, windrow and sweep were 
used for the harvest. A 5 to 7 kg random sample was collected from 
each plot to determine wet and air dry (45 C) weights, trash percent- 
age, and shelling percentage. Average kernel weights were deter- 
mined from three separate 100-kernel samples from each treatment 
block. 


Single-degree-of-freedom orthogonal contrasts were used to de- 


termine irrigation treatment effects on growth and yield responses 
(SAS Institute 1985). 


Results and discussion 
Irrigation quantity 


Applied water plus rainfall for the six irrigation treat- 
ments for 1985 through 1987 are shown in Table 2. Spring 
and summer trickle irrigation for the treatment designed 
to receive 100% of estimated ET, (T3) increased from a 
total of 13.1 m? tree~! in 1985 to 19.2 m? tree~! in 1986 
to 26.7 m? tree~ ! in 1987 due to modifications in the crop 
coefficient to account for both increases in tree canopy 
size and changes in the calculated soil water use. 

Data in Fig. 1 illustrate the range of volumetric water 
content in the soil profile across several treatments at 
neutron access tube positions 1.5 m east of the instru- 
mented trees in the trickle lateral direction and 1 m south 
of the lateral. Higher water applications associated with 
increases in tree shaded area in 1986 and 1987 and more 
winter rainfall than 1985 resulted in higher early- to mid- 
summer soil water contents with depth in all treatments 
than in 1985. In specific areas of the root zone, similar 
treatment differences in patterns of soil water content 
persisted during much of each year due to low annual 
rainfall and consistent treatment differences in applied 
water. Soils at the experiment site have been previously 
identified as having average volumetric water contents of 
about 0.35 and 0.13 at —0.03 MPa and —1.5 MPa soil 
matric potential, respectively (Nielsen et al. 1973; Sharples 
1986). 

Soil water depletion and deep percolation in this ex- 
periment were discussed by Nightingale et al. (1991) and 
Rolston et al. (1991). Rolston et al. (1991) described the 
1986 soil water depletion between days 100 and 270 (just 
after harvest) as less than 2000 | tree~' and 3000 | tree~! 
in the 50% ET, and 175% ET, treatments, respectively. 
These amounts represent about 20% and less than 10% 
of the total applied water for the year in the 50% and 
175% ET, treatments, respectively. Season-total soil wa- 
ter use estimates were of similar magnitude in treatments 
measured in 1985 and 1987. 

Fereres et al. (1981) procedure for calculating imma- 
ture tree water use, which was used in this study, is based 
on ET, rates determined under flood irrigation for 
mature trees with a correction based on the percent of 
shaded ground area. In flood irrigation the water is ap- 
plied over a large ground area per tree, while the single 
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Fig. 1a,b. Mean soil volumetric water content as a function of depth in the soil profile and irrigation treatment for trickle irrigated 
almonds in (a) July, 1985; and (b) July, 1986 


Table 2. Rainfall and trickle-applied water as a function of irrigation treatment in 1985, 1986, 1987 for trickle irrigated almonds at the West 
Side Field Station near Five Points, CA 





Treat- Irrigation Year 
ment rate 
code (as % 1985 1986 1987 

of ET,) 














Postharvest Main season Postharvest Main season Postharvest Main season 
irrigation irrigation” irrigation irrigation irrigation irrigation 
plus rain* (mm) plus rain (mm) plus rain (mm) 

(mm) (mm) 





50 Loy 171 267 440 368 
75 136 261 395 431 546 
100 141 347 511 431 716 
125 140 441 647 414 897 
150 135 525 761 415 1085 
175 t35 597 885 403 1258 





* Irrigation water applied after the previous year’s harvest (usually in October-November of previous year) plus rainfall through September 
harvest of the year shown; post-harvest irrigation amount averaged 70 mm, 85 mm, 159 mm in 1984-1985, 1985-1986, and 1986-1987, 
respectively 


> Water applied between February and September of each year 


trickle line per tree row in this study only resulted in wet —_in considerable deep percolation losses, as suggested by 
surface areas extending 0.5 to 1.25 m beyond the trickle Nightingale et al. (1991). The relative size of the wetted 
line location. The mismatch between the small infiltra- zone and infiltration area could be changed with modifi- 
tion area of the single trickle line per tree row and the "cations of emitter application rates and irrigation fre- 
relatively high water application rates could have resulted § quency (Jury and Earl 1977). Such modifications may be 





Table 3. Almond tree trunk cross-sectional area determined from 
measured diameters as a function of irrigation treatment for the 
Ruby variety at the West Side Field Station in 1983 through 1987. 


Orthogonal comparisons are shown for major treatment compari- 
sons 





Irrigation 
treatment 
code Year 


Mean trunk cross-sectional area* (cm?) 








1983 





89.9 
100.3 
103.9 
113.1 
102.1 
103.9 


Comparisons” 

Ti VS. T2 * ** OK ** and 
Ft VS. 53 ** ** ** ** ** 
T2 vs. T3 NS * bales * * 
F2 Vs. T4 * ** ** ** ** 
T3 vs. T4 NS +¢ ++ + +% 
T3 Vsi ES NS iy ed =% + 
T3 vs. T6 NS *¢ +* +* ” 
T4 vs. TS NS NS NS NS NS 
T4 vs. T6 NS NS NS NS NS 
T5 vs. T6 NS NS NS * NS 





* Mean of 35 trees (7 blocks, 5 per block) after fall leaf drop 
(between day 307 and day 337, depending on the year) 

> Orthogonal comparisons within a year designated as NS, *, or ** 
were not significant, significant at the P<0.10 level, or significant 
at the P<0.05 level, respectively 


worth considering in attempts to limit deep percolation 
losses and better match short-term water storage to crop 
water use. 


Tree growth 


Just prior to initiation of the 1984 treatments, trunk 
cross-sectional areas in December, 1983 were significant- 
ly smaller in new treatment T1 trees than in all other 
treatments (Table 3). During 1981 through 1983, treat- 
ment T1 was a low water application treatment (W,N,), 
and the smaller trunk cross-section area in December, 
1983 indicated a carry-over effect of the prior treatment 
on the 1984 through 1987 treatment. Initial trunk cross- 
section area of treatment T4 was also significantly greater 
than T2, but not significantly different from other treat- 
ments. December 1983 trunk cross-section areas of the 
other treatments were not significantly different. 

The long-term relative effect of specific treatments can 
be assessed by looking at the consistency of treatment 
effects relative to the 100% ET, treatment (T3). The ratio 
of trunk cross-section areas of low water application 
treatments (T1, T2) to the 100% ET, treatment (T3) gen- 
erally declined from 1983 through 1987 (calculated from 
Table 3). The ratio for high water application treatments 
(T4 through T6) started near 1.0 to 1.1 in 1983 and in- 
creased to 1.2 to 1.3 in subsequent years. In treatments 
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Fig. 2. Increase in trickle-irrigated almond tree trunk cross-sec- 
tional area of all irrigation treatments during the December, 1983 
through December, 1987 period as a function of cumulative post- 
harvest plus irrigation season trickle-applied water. Linear regres- 
sion equation: y=0.035 x—8.39; r?=0.98, n=30 


T2 through T6, increases in trunk cross-section areas in 
1984 through 1987 were roughly in proportion to applied 
water, and the ratio of treatment cross-section areas to 
treatment T3 areas were also positively correlated with 
applied water. The potential carry-over effects of the 
pre-1984 treatment cannot be separated from post-1984 
responses in treatment T1. 

Other than treatment T1, tree trunk cross-section areas 
in 1983 were not significantly affected by the pre-1984 
nitrogen subtreatments within each irrigation level (data 
not shown). 1984 is a difficult year to assess in terms of 
treatment effects. Even though 1984 was described as 
transitional between the two experiments, treatments in 
1984 were irrigated according to the same 50 to 175% 
ET, rates scheduled for the 1985 through 1987 seasons. 
Trunk growth rates in 1984 were roughly proportional to 
applied water (Table 3). 

As shown in Fig. 2, a linear relationship (r? =0.99) 
was observed between the increase in trunk cross-section- 
al area from 1985 through 1987 and the total amount of 
irrigation water. Although this linear relationship implies 
a direct response to increases in applied water, it is prob- 
ably related to a combination of factors including the 
increased wetted-soil volume created by the higher ap- 
plied water volumes, increased root activity, and perhaps, 
access to more soil nutrients within a larger rooting vol- 
ume. 

Scaffold branch growth expressed as a mean increase 
in branch diameter for five two-month periods in 1987 
is presented in Table 4. The highest rate of growth 
(um day~') occurred from April 16 to June 3 for all irri- 
gation treatments. As water stress intensified in mid-sum- 
mer, branch growth rates declined more in treatments T1 
and T2 than in higher water treatments. 

The increases in shaded ground area during the three 
growing seasons are illustrated in Fig. 3 for the T3 (100% 
ET.) treatment. The within-season changes in scaffold 
branch diameter were not linear, following a different 
seasonal pattern than changes in shaded area (Table 4, 
Fig. 3). Since shaded ground area as measured cannot 
distinguish between intensities of shading nor detect mul- 
tiple layers of leaves, it is not necessarily a good estimate 
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Fig. 3. Mean shaded ground area (expressed as percentage of total 
ground area of 37.16 m? per tree) of treatment T3 (100% ET,) 
almond trees as a function of day of year in 1985, 1986, and 1987. 
Linear regression equations were as follows: (1985) y=0.024 x+ 
17.4, r?=0.91, n=18; (1986) y=0.045 x+17.4, r?=0.89, n=18; 
(1987) y=0.042 x +24.3, r?=0.88, n=18 


Table 4. Average daily increase in the diameter of the primary scaf- 
fold branches of almond trees as a function of trickle irrigation 
treatments for selected time periods in 1987 





Date Increase in mean diameter (um day~ ') 








From Trickle irrigation treatments 





T1 12 “3s “i4 "TS RS 





Feb. 5 Apr. 16 SOP EO ASD DS ES 
™ ®@ ®@ @® OD @) 
Jun. 3 44 52 46 47 48 59 
6 @” @& @ ©) (Qi) 
Aug. 12 13 30 42 43 47 53 
6 @3) @d2) © 4 @©® 
Sep. 28 2 9 17 15 17 19 
4 @ @® ® ® & 
Dec. 18 ao. Seo, SS: eS Vee Se 
og ©) @®@ @ @ @® 


Apr. 16 
Jun. 3 
Aug. 12 


Sep. 28 





* Mean of 2 primary scaffold branches per tree, 3 trees per treat- 
ment plot in each of Blocks 3, 4, and 5; n=18 per treatment for each 
measurement date; standard error of the mean shown in parenthe- 
ses under the mean for each treatment and date 


of within-season changes in total tree leaf area. Rather, it 
is more closely related to changes in canopy diameter. 

Between April of 1984 and July of 1987, ground areas 
shaded by trees (as a percent of the 37.2 m? of land allo- 
cated each tree) changed from 16% to 23% in treatment 
T1, from 17% to 27% in T2, from 17% to 35% in T3, 
from 19% to 36% in T4, from 20 to 40% in T5, and from 
18% to 41% in T6. The only significant difference in 
initial 1984 shaded ground areas was between T1 versus 
T4 and TS. As with trunk cross-section areas, the lower 
initial shaded ground area of treatment T1 indicated a 
significant carry-over effect of the 1980 through 1983 
treatment. 

Under the conditions of this experiment, significant 
leaf drop (resulting in more than a 5% reduction in tree 
shaded area) was not observed prior to the cessation of 
irrigation in any year even in the T1 or T2 treatments. 


Table 5. Mean leaf water potential for each of the months of May, 
June, July and August of fully-illuminated spur branch almond 
leaves on the southwestern quadrant of trees within 30 cm of branch 
ends as a function of irrigation treatment in 1985, 1986, and 1987 





Month _Irri- 

of gation 

year treat- Year 
ment 


Leaf water potential (MPa) 








1985 1986 1987 





Mean* S.E.*. Mean S.E. Mean _ S.E. 





—2.58 0.05 0.12 
—2.40 0.04 ; 0.08 
—2.17 0.08 A 0.06 
—2.08 0.06 : 0.09 
—2.02 0.02 . 0.04 
—1.94 0.08 ; 0.06 


—3.07 0.04 4 0.08 
—2.96 0.09 : 0.04 
—2.74 0.13 f 0.05 
—2.51 0.11 : 0.04 
—2.42 0.10 , 0.04 
—2.43 0.07 : 0.04 


—3.12 0.06 : 0.11 
—3.07 0.14 : 0.03 
—2.68 0.09 : 0.04 
—2.33 0.06 i 0.09 
—2.22 0.10 —2.31 0.07 
—2.15 0.09 —2.27 0.13 


—3.20 0.15 —2.81 0.08 
—3.02 0.06 —2.63 0.04 
—2.65 0.06 —2.42 0.09 
—2.35 0.06 —2.35 0.07 
—2.26 0.12 —2.29 0.11 
—2.21 0.08 —2.27 0.06 


—2.61 0.10 
—2.57 0.06 
—2.39 0.03 
—2.21 0.04 
—2.09 0.07 
—2.16 0.08 


—2.85 0.11 
—2.69 0.05 
—2.52 0.06 
—2.34 0.05 
—2.22 0.07 
—2.17 0.06 


—2.92 0.12 
—2.85 0.08 
—2.59 0.05 
—2.31 0.07 
—2.19 0.04 
—2.22 0.06 


—3.08 0.15 
—2.89 0.12 
—2.64 0.07 
—2.49 0.02 
—2.29 0.05 
—2.33 0.05 





* Mean and standard error of the mean (S.E.) for leaves sampled on 
from two to three separate days per month. On each treatment 
within each sample date, one spur branch leaf per tree on each of 
four separate trees was sampled in field blocks 3, 4, and 5; n=12 per 
treatment 


However, at harvest (12 to 19 days after irrigation cut- 
off), leaf drop was significantly higher (25 to 40% of total 
leaves) in treatments T1 and T2 than in other treatments 
(<10% of total leaves). 


Plant water status 


Average leaf water potentials (LWP) of fully-illuminated 
leaves across irrigation treatments exhibited similar 
trends during each year of the study (Table 5). At water 
applications below the 150% ET, level, LWP declined 
progressively (Table 5). LWP of treatments T5 and T6 
were not significantly different on any measurement date 
during the entire study. March and April LWP values 
were higher than —2.5 MPa in all treatments (data not 
shown). Late-season (July, August) leaf water potentials 
were higher in treatments T1, T2 and T3 in 1986 than in 
the other two years. Average soil water contents of these 
treatments were 5 to 8% higher in late July and mid- 
August in 1986 than during similar periods in the other 
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Fig. 4. Leaf conductance to water vapor (g,) of fully-illuminated 
almond leaves as a function of day of year and irrigation treatment 
in 1987. Standard error of the means are shown as vertical lines only 
for treatments T1 and T6 


Table 6. Mean almond nut retention for three levels of irrigation in 
1985, 1986, 1987 as a function of flower counts and early nut set at 
the West Side Field Station near Five Points, CA 





Year Day Nut retention as 
of percent of initial 
year flower count* (%) 


Nut retention as 
percent of initial 
nut count” (%) 





Treatment Treatment 





T1 4 TS EE ys TS 





24 32 35 56 62 64 
(2) =) (3) 3) @& @ 
20 28 30 46 55 55 
(2) (3) (3) 3) @ © 
30 ND‘! 26 75 75 73 
(4) “ (3) 6) @) @ 
232 28 23 26 71 72 74 
(3) (2) (3) 4) 38) @ 
160 20 23 27 67 73 75 
(4) (2) (3) (7) (4) 6) 
19 22 26 63 70 73 
(3) (2) (1) 56) 8) @® 
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* Mean of flower counts made on calendar days 58 to 68 in 1985, 
55 to 70 in 1986, and 61 to 73 in 1987. Mean of counts on two 
marked branches on each of 3 trees per treatment in blocks 3, 4, 
and 5; n=18 branches 

> Mean of nut counts on two marked branches on each of 3 trees 
per treatment plot in blocks 3, 4, and 5 (n=18 branches) measured 
initially on day of year 112, 91, and 105 in 1985, 1986, and 1987, 
respectively 

© Values shown in parentheses below mean values are standard 
error of the mean 

4 ND=no data 


two years, perhaps reflecting higher winter rainfall in 
early 1986. Lower prevailing afternoon vapor pressure 
deficits during most July and August measurement dates 
in 1986 may have also contributed to lower leaf water 
potentials than in similar periods in 1985 and 1987. Mid- 
season LWP vaues of high water application treatments 
in this study were 0.1 to 0.3 MPa lower than reported for 
well-watered almonds in several trickle irrigation studies 
under conditions of coarser soils and lower ET, (Ruiz- 
Sanchez et al. 1988; Castel and Fereres 1982). 
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The slowly-developing water deficits produced with 
continuous deficit irrigation in the low water treatments 
of this study are in contrast to the long-term (one to two 
month duration) periods without irrigation investigated 
in studies by Uriu et al. (1970), El-Sharkawi and El-Mon- 
ayeri (1976), Castel and Fereres (1982), and Krauter 
(1985). Plant responses to even the most severe continu- 
ous deficit irrigation treatments in this study were princi- 
pally reduced overall plant size and growth, with little 
leaf shedding or changes in leaf orientation during the 
irrigation season, and maintenance of at least partially 
open stomata in all treatments (Fig. 4). Leaf conductance 
to water vapor (g,) of treatments T1 and T2 were signif- 
icantly lower than other treatments on most measure- 
ment dates. Leaf conductance declined linearly with de- 
creasing LWP, similar to Torrecillas et al. (1988), al- 
though no consistent threshold LWP associated with re- 
duced g, could be identified. During the month of July in 
all three years, the minimum g, measured in illuminated 
spur branch leaves in treatment T1 and T6 averaged 0.22 
and 0.62 mol m~? s~', respectively. 

Tree growth and nut yields are a composite response 
to plant water status over a long period, and while several 
water-deficit sensitive growth stages exist for almond 
(Micke and Kester 1978), a critical LWP or g, value below 
which yield or long-term growth will be affected has not 
been identified (Castel and Fereres 1982; Ruiz-Sanchez 
et al. 1988). In this study, however, maintenance of LWP 
higher than —2.3 MPa during early nut development 
(March through May) and greater than —2.5 MPa the 
remainder of the irrigation season was correlated with 
sustained higher vegetative growth rates and nut yields. 
Correlations between LWP and yield or trunk cross-sec- 
tion areas within any year were positive and significant, 
with all r? values in excess of 0.9. 


Almond nut retention 


In 1985 and 1987, nut retention increased with increases 
in applied water, although no significant differences oc- 
curred across treatments in 1986. Additional nut loss 
during the period from June through a time just before 
harvest in any year ranged from near zero to 18% 
(Table 6). 

Nut retention was also evaluated with respect to the 
initial nut counts taken in April of each year to determine 
nut loss during the summer. In 1985 and 1987, the 50% 
ET, irrigation treatment (Table 6) and 75% ET, treat- 
ment (data not shown) had significantly higher nut losses 
than the other treatments, while no significant differences 
occurred in 1986. Trickle irrigation at levels above 100% 
ET, reduced nut loss by only 2% to 4% over the three- 
year study. 


Nut meat yield and kernel weight 
In 1985 and 1987, nut meat yields increased significantly 


with increases in applied water from 50% ET, through 
125% ET, (Table 7). Results were mixed in 1986, with 
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Fig. 5. Almond nut meat yields as a function of total applied water 
plus rain in 1985, 1986, 1987. Second-degree polynomial equations 
fitted to data for individual years were as follows: (1985) y=51.22+ 
2.45 x—0.0013 x”, r?=0.993; (1986) y=793.4+ 1.04 x—0.0003 x?, 
r?=0.958; (1987) y=277.4+1.99 x—0.0007 x”, r?=0.993. Stan- 
dard error of the means are shown as vertical lines for each data 
point 


Table 7. Mean almond (Ruby variety) nut meat yield and kernel 
weight as a function of trickle irrigation rate. Orthogonal compari- 
sons within each year are shown for major treatment comparisons 





Irrigation 
treatment 
code 


Mean nut meat yield 
(kg ha~*)* 


Mean kernel weight 
(mg kernel ~')” 





Year Year 





1983 1985 1986 1987 





T1 340 594 862 787 
T2 374 677 901 881 
T3 392 742 901 918 
T4 388 781 892 953 
12 472 828 940 946 
T6 405 864 946 974 


Comparisons * 


Ti Vs. T2 ** 2 * OR ** 
Ti VS. 3. ** ** aK OK ** 
T2 vs. T3 nS * - NS * 

T2 vs. T4 Ne ee we» ONS. ee 
T3 vs. T4 nS, * NS NS NS 
T3 vs. TS ~ ony ** NS NS 
T3 vs. T6 i SiN SNS SNS 
T4 vs. TS NS NS NS NS NS 
T4 vs. T6 - . “4 NS NS 
TS vs. T6 NS NS NS NS NS 





* Based on the mean nut yield of 5 trees harvested in each of 7 blocks 
per treatment. Planting density equivalent to 268 trees per ha 

> Based on the average kernel weight of 100 kernels from a 4 to 7 kg 
subsample of harvested nuts from 5 trees per treatment in each of 
7 blocks per treatment 

© Orthogonal comparisons within a year designated as NS, *, or ** 
were not significant, significant at the P<0.10 level, or significant 
at the P<0.05 level, respectively 


higher yields in 150% and 175% ET, treatments (T5, T6) 
than at lower ET, levels. The second degree polynomial 
prediction equations for nut yield as a function of total 
applied water plus rain during each year are presented in 
Fig. 5. Although the polynomial curve fits (Fig. 5) seem 
to indicate that nut yields continue to increase with water 


applications greater than 150% ET,, there were no con- 
sistent, significant increases in indices of tree growth or 
nut yields at water applications in excess of 150% ET, 
(Tables 3, 7). 

In terms of water use, mature almond trees were de- 
fined by Micke and Kester (1978) as those shading more 
than 70% of the ground area, and they reported an aver- 
age annual ET, in the San Joaquin Valley as 40.6 m? 
tree” '. Water applications to treatment T3 trees with a 
shaded ground area of 30 to 32% in this study were 65 to 
70% of estimated ET, for a mature orchard, as described 
using the crop coefficient procedures of Fereres et al. (1981). 
Even though they did not approach the 70% ground area 
shading of a mature orchard, trees irrigated at the 150% 
and 175% ET, levels received applied water equal to or in 
excess of mature almond ET, (Micke and Kester 1978). 

Tree growth, nut yields and plant water status mea- 
surements all indicate that the 100% ET, treatment was 
deficit-irrigated, with reductions in yield and growth rel- 
ative to the 150% and 175% ET, treatments. These find- 
ings appear contradictory to previously reported work of 
Nightingale et al. (1991) from this same study showing 
that the leaching fraction in the upper 1 m of the soil 
beneath the trickle lateral out to about 1m from the 
lateral was 10 to 22% of applied water. The existence of 
significant leaching in the 100% ET, treatment, where 
plant water deficits resulted in reduced tree growth and 
yields, is somewhat difficult to reconcile. It is possible 
that irrigations were not correctly scheduled to meet 
plant water needs over time, resulting in over-irrigation 
and leaching during some periods and under-irrigation 
during critical periods where water deficits influenced 
growth and yield. LWP data does not support this con- 
clusion, however, since treatments T1 through T3 had 
consistently lower LWP than higher water application 
treatments during all measured periods from May through 
August (Table 5). Another possibility is that significant 
leaching losses in some areas of the wetted soil volume 
resulted from the mismatch of the small infiltration area 
of the single trickle line and the relatively high trickle 
application rates in this study, as suggested by Nightin- 
gale et al. (1991). This second possibility is more compat- 
ible with the coexistence of plant water deficits in treat- 
ments T1 through T4 and significant leaching losses iden- 
tified by Nightingale et al. (1991). In future studies it may 
be worthwhile to determine if changes in irrigation fre- 
quency, application rates, or addition of more laterals to 
wet a larger soil volume can more evenly distribute water 
applications over time and space and avoid some of the 
leaching losses suggested in this study. 

Nightingale et al. (1991) previously described the in- 
fluence of long-term use of this quality of water on soil 
accumulations of boron and other irrigation water con- 
stituents. Problems of accumulation of salt and boron at 
the perimeter of the wetted soil volume were identified, 
particularly in the low water application treatments. The 
irrigation water in this study was not found to cause a 
serious salinity or boron hazard in the short-term (several 
years) (Nightingale et al. 1991). However, if minimal 
leaching fractions (less than 0.05 to 0.10) were achieved, 
salt and boron could accumulate within the root zone in 





Table 8. Water use efficiency (determined as kg of kernels per m? of 
applied water) in trickle irrigated almonds at the West Side Field 
Station in 1985 through 1987. Orthogonal comparisons are shown 
for major treatment comparisons 





Irrigation 
treatment 
code Year 


Water use efficiency (kg m~*)* 








1985 





0.257 
0.249 
0.225 


Comparisons” 


T1 vs. ss » 

T1 vs. 7 + 
T2 vs. NS * 

T2 vs. ss “7 
T3 vs. = NS 
T3 vs. _ * 

T3 vs. = +* 
T4 vs. NS NS 
T4 vs. bes = 
T5 vs. NS NS 





* Mean of 35 trees 

> Orthogonal comparisons within a year designated as NS, *, or ** 
were not significant, significant at the P<0.10 level, or significant 
at the P<0.05 level, respectively 


detrimental amounts over a period of several years. The 
threshold salinity for growth reductions in almond is 
1.5dS m~' (Maas 1986). Accumulations of boron and 
salinity within the root zone in the low applied water 
treatments (T1, T2) could eventually interact with the 
effects of deficit irrigation in reducing nut yields and size. 
Hoffman etal. (1989) observed both short-term and 
long-term growth and yield reductions in plums, another 
Prunus spp. (Prunus salicina), at average root zone salin- 
ities lower than those observed in this study. 

The mean almond kernel weights as a function of year 
and irrigation treatment are shown in Table 7. Kester 
et al. (1985) stated that the Ruby cultivar produces medi- 
um-sized kernels, 945 to 1135 mg kernel~ '. During 1985, 
each incremental (25% of ET.) increase in applied water 
resulted in a significant increase in kernel weight. Since all 
1985 kernel weights were lower than the characteristic 
values stated by Kester et al. (1985), the combination of 
a mid-summer 20-day period with peak daily tempera- 
tures 3 to 5°C above normal and more severe water 
deficits in early summer (Table 5) may have reduced ker- 
nel weights in all treatments. In 1986 and 1987, kernel 
weights were significantly smaller only in the low ET, 
treatments (T1 in 1986, Ti and T2 in 1987). While a 
variety of factors could be involved in the response of 
kernel weights to applied water, even moderate water 
stress during the embryo growth phase has been associat- 
ed with reductions in the expansive growth of almond 
kernels (Micke and Kester 1978). 
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Water use efficiencies (WUE), kg nut meat per m? of 
trickle applied water, were calculated for each of the 
3 years (Table 8). For the 3-year means, WUE ranged 
from a high of 0.296 for the lowest water application 
treatment (T1) to a low of 0.163 for the highest irrigation 
level (T6), with significantly higher WUE in the lower 
three water application treatments (T1, T2, T3) than in 
higher applied water treatments in at least two out of 
three years. 


Conclusions 


During the 1985 through 1987 period of this study, leaf 
conductance, leaf water potential, and soil water content 
data indicated relatively consistent differences in soil and 
plant water status across irrigation treatments. Deficit 
irrigation resulted in significantly smaller trees and re- 
duced production of leaves and fruiting sites. Tree growth, 
nut yield, and plant water status indicators improved 
significantly with increases in applied water up to the T4 
(125% ET,) level, and did not generally increase signifi- 
cantly at the 150% and 175% ET, levels. 

Differential treatments prior to 1984 resulted in signif- 
icantly smaller trees (trunk diameter, canopy size) in the 
50% ET, treatment at the initiation of the new experi- 
ment when compared to other treatments. These carry- 
over effects of the prior experiment prevented identifica- 
tion of the absolute responses to the 1984 through 1987 
50% ET, treatment. The initial trunk cross-section area 
of T4 was also significantly greater than treatment T2. 
No other treatments were significantly different in mea- 
sured growth parameters at the time of initiation of the 
new experiment. 

Annual tree growth and nut meat yield were positively 
correlated (r7>0.9) with averages of afternoon LWP 
measured during May through August. Sustained LWP 
of higher than —2.3 MPa during early nut development 
(March through May) and higher than —2.5 MPa the 
rest of the irrigation season were correlated with the 
highest growth rates and nut meat yields. 

Water applied in the 100% ET, treatment, determined 
as described by Fereres (1981) with additional adjust- 
ments according to soil water depletion data, was not 
adequate to maximize yields and maintain non-water- 
stressed trees under the conditions of this trickle irriga- 
tion experiment. The water application amounts in the 
125% or 150% ET, treatments may more closely repre- 
sent irrigation needed to maximize growth and yield un- 
der conditions (climate, soil, irrigation system manage- 
ment) at the experiment site. Rolston et al. (1991) and 
Nightingale et al. (1991) previously defined the uncer- 
tainty involved in quantifying actual ET, and deep perco- 
lation for this experiment. If the irrigation system and 
management could be improved to better store applied 
water within the active root zone, the high leaching losses 
described by Nightingale et al. (1991) would be avoided. 
With such improvements, yield and growth might be 
maximized at lower irrigation amounts (closer to 125% 
of ET, calculated by methods of Fereres (1981)). 
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Abstract. Total salinity and K/Na ratios in the irrigation 
water were varied for two different types of plants; sor- 
ghum, a moderately salt tolerant crop plant whose 
growth is inhibited by high NaCl levels, and Atriplex, an 
extreme halophyte whose growth is stimulated by high 
NaCl levels. Sorghum growth was significantly reduced 
at all salinity levels from 50 to 150 mM, but Atriplex 
growth was not reduced until salinity exceeded 100 mM. 
In both species, growth reductions were increased by in- 
creasing the ratio of K to Na from 1/100 to 1/1 in the 
irrigation water. The amount of K and Na accumulated 
in the leaves of Atriplex reflected the relative amounts in 
the irrigation water, but in sorghum most of the increase 
in the concentration of K + Na was due to the increase 
in K regardless of the ratio of K/Na in the irrigation 
water. Nitrate levels in leaves of sorghum were little af- 
fected by salinity but were decreased in Atriplex by the 
combination of high salinity and K/Na ratio. The lower 
amount of salt accumulation in sorghum was compensat- 
ed for by greater accumulation of soluble organic com- 
pounds, such as proline and soluble carbohydrates, that 
presumably were osmotically active and could contribute 
to osmotic adjustment. All of the results supported the 
general conclusion that, despite differences between gly- 
cophytes, such as sorghum, and halophytes, such as 
Atriplex, in sensitivity to the total soluble salt concentra- 
tion in irrigation water, in both species growth was re- 
duced more by K/Na ratios of 1/1 and 1/10 than by a 
ratio of 1/100 at those salinity levels that significantly 
reduced growth in each species. 


* Present address: Soil Salinity Laboratory, Alexandria, Bacos, 
Egypt 

** Present address: Bioresources Research Facility, University of 
Arizona, 250 E. Valencia Road, Tucson, AZ 85706, USA 
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The solutes responsible for increasing salinity in a plant’s 
environment usually are Na salts. Plants growing in such 
environments are basically of two types, glycophytes and 
halophytes. Glycophytes are relative excluders of the Na, 
and growth is reduced by relatively low levels of Na in the 
leaves. In contrast, halophytes are relative accumulators 
of Na, and many of them show enhanced growth at rela- 
tively high levels of Na in the leaves. In halophytes, Na is 
the major solute contributing to decrease in leaf osmotic 
potential in response to increased salinity in the root en- 
vironment. Since glycophytes exclude most of the Na, K 
is the primary inorganic solute that is accumulated in 
leaves in response to increased external salinity. Howev- 
er, in spite of the similarity between Na and K, they do 
not elicit the same responses in plants (Hylton et al. 1967; 
Aslam 1975; Helal et al. 1975; Eshel 1985; Matoh et al. 
1986). Potassium is the only monovalent cation that is 
essential for growth of all higher plants. Sodium has only 
been shown to be essential for some species, but appears 
to stimulate the growth of others (Brownell 1967; Brownell 
and Crossland 1972; Rush and Epstein 1981; Eshel 1985). 
In some glycophytes, such as sorghum (Weimberg et al. 
1984) and mung bean (Salim and Pitman 1983), high 
levels of K were found to be more harmful than high 
levels of Na, but in other cases addition of K salts to the 
nutrient medium reduced the harmful effects of high Na 
levels (Hylton et al. 1967; Helal et al. 1975). In halo- 
phytes, high K levels have been speculated to cause simi- 
lar responses (Flowers et al. 1977; Flowers and Lauchli 
1983). However, since the response of both a glycophyte 
and a halophyte to high K levels at different salinities has 
not been determined under similar conditions, it is not 
valid to conclude that high external levels of K have the 
same effect on both types of plants. 

This study was designed to compare the effects of 
different salinity levels and K/Na ratios on growth and 
solute content of two different types of plants, a moder- 
ately salt tolerant crop plant (Sorghum), whose growth is 
inhibited by high NaCl levels, and an extreme halophyte 
(Atriplex), whose growth is stimulated by high NaCl 
levels. The hypothesis tested was that high levels of K in 
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the root environment will not reduce growth in the halo- 
phyte as much as in the glycophyte, as long as Na is 
present in equal or greater amounts. 


Materials and methods 
Experimental procedure 


Seeds of Sorghum bicolor L. ““BTX 62” were sown in 10 liter plastic 
pots containing 12 kg air-dried and sieved sand and irrigated with 
a modified Hoagland solution (see below). Rooted stem cuttings of 
Atriplex amnicola were transplanted directly into such pots, one per 
container, and irrigated with the same solution. After 15 days, the 
sorghum seedlings were thinned to 5 plants per pot, and the treat- 
ments began. 

Plants were grown in an evaporately-cooled greenhouse in 
Tucson, Arizona under high natural radiation levels (ca. 800- 
1050 Wm~?), high temperatures (ca. 20—24°C nights/26-30°C 
days) and moderate humidity levels (70-80% RH). See Glenn 
(1984) for details of seasonal environmental conditions in a similar 
greenhouse at the same location. The modified K-free Hoaglands 
solution was prepared by replacing KH,PO, with NaH,PO, and 
KNO, with additional Ca(NO,),. The Ca and NO, concentrations 
were 5 mM and the PO, concentration was 2 mM. Known amounts 
of KCl and NaCl were added to give final concentrations of 
(K + Na) of 4 (control), 100, 200, and 400 mmol for Atriplex and 4 
(control), 50, 100, and 150 mmol for sorghum. At each level of 
salinity there were three different K/Na ratios of 1/1, 1/10, and 
1/100, so the total number of treatments was 12 for each species. 
The treatments were laid out in a randomized block design with 4 
replicates. 

Pots were weighed every other day to determine the amount of 
evapotranspiration. The pots were watered with sufficient amounts 
of the treatment solutions to restore the soil to field capacity plus 
an additional amount equivalent to twice the amount of evapotran- 
spiration to insure adequate leaching. Drainwater from each pot 
was quantitatively collected and the electrical conductivity mea- 
sured to insure that the desired salinity levels in each pot were 
maintained. 


Plant harvest 


Plant shoots were harvested 50 days after the initiation of the treat- 
ments. Fresh and dry weights were determined, and subsamples 
(1 g) of fresh leaves were frozen and later used for organic analyses. 
After harvesting, the experiment was repeated for Atriplex with the 
same treatments and experimental measurements. 


Inorganic analysis 


The oven-dried material was dry ashed at 550°C for 5 h, the ash was 
dissolved in dilute HCl (Chapman and Pratt 1961), and Na, K, Ca, 
Mg were assayed by atomic absorption spectroscopy. Oven-dried 
material was ground and extracted for 0.5 h with hot distilled water, 
and nitrate was determined colorimetrically with a Technicon auto- 
analyzer. 


Organic analysis 


Thawed fresh material was extracted with methanol:chloroform: 
water, 12:5:1 (v:v:v) (Singh et al. 1973). Proline was determined 
with the acid ninhydrin method of Troll and Lindsley (1955), except 
that toluene was used instead of benzene. Free amino acids were 
determined with ninhydrin reagent (Rosen 1957). Glycine was used 
as the standard, so the results are presented as glycine equivalents. 


Total soluble carbohydrate was determined with a mixture of phe- 
nol and H,SO, (Dubois et al. 1956). Since glucose was used as the 
standard, the results are presented as glucose equivalents. 


Water potential 


Water potentials of leaf discs were determined psychrometrically 
with a Wescor Model HR33T microvoltmeter. 


Results 
Fresh and dry weights 


The effects of salinity and K/Na ratios on shoot fresh and 
dry weights of Atriplex were significant at the 5% level, 
but the interactions were not (Table 1). Regardless of the 
K/Na ratio, the 100 mmol treatment did not significantly 
affect fresh weight in experiment I, but it did have a sig- 
nificant effect in experiment II. However, further increas- 
es in salinity reduced fresh weight significantly in both 
experiments. At all salinity levels, decreasing the K/Na 
ratio increased the fresh weight in both experiments. 

Similar to fresh weight, the dry weights were not sig- 
nificantly reduced by the 100 mmol treatment (regardless 
of K/Na ratio), but further increases in salinity signifi- 
cantly reduced dry weight. Decreasing the K/Na ratio 
from 1/1 to 1/100 within each salinity level significantly 
increased dry weight in both experiments. 

Fresh and dry weights of sorghum decreased signifi- 
cantly with all increases in salinity level, regardless of 
K/Na ratio (Table 2). There was a significant interaction 
between salinity level and K/Na ratio in sorghum. Conse- 
quently, decreasing the K/Na ratio did not have as clear 


Table 1. Effect of salinity and K/Na ratio on shoot fresh and dry 
weights of Atriplex amnicola 





Salinity treatment K/Na 


Fresh weight 


Dry weight 


i* I 





(K+ Na) mmol/l g/plant 


4 (control) 240 
270 
278 


242 
263 
266 


161 
176 
197 


75 
93 
116 





L.S.D. (0.05) 
For salinity (S) 27 
For K/Na ratios (R) 24 





* Experiment I 
** Experiment II 





Table 2. Effect of salinity and K/Na ratio on fresh and dry weights 
of Sorghum bicolor 
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Table 3. Effect of salinity and K/Na ratio on leaf water potential of 
Atriplex and Sorghum 





Salinity treatment K/Na Fresh weight Dry weight 





(K + Na) mmol/l g/plant 


4 (control) 1/1 96.0 18.4 
1/10 93.0 18.1 
1/100 104.1 19.7 


1/1 74.9 14.9 
1/10 66.5 13.0 
1/100 64.8 14.6 


1/1 42.8 8.6 
1/10 48.7 9.8 
1/100 329 10.3 


1/1 21.9 4.4 
1/10 27.6 5.5 
1/100 6.6 





L.S.D. (0.05) 

For salinity (S) 

For K/Na ratios (R) 
ForSxR 





an effect with sorghum as with Atriplex, except at the two 
highest salinity levels where decreases in the ratio signifi- 
cantly increased both fresh and dry weights. 


Leaf water potential 


In both species, the leaf water potential decreased with 
increased salinity (Table 3). The decline was greater in 
Atriplex than in sorghum, both in absolute value and 
percent decrease. For example, at the 100 mM level, the 
absolute value was over 50% lower than at the control 
level in Atriplex while in sorghum the difference was less 
than 40%. Water potential was not significantly affected 
by the K/Na ratio, but rather seemed to respond to total 
salinity. 


Ash content 


At all salinity levels, including the relatively non-saline 
control level, the ash content of Atriplex was consider- 
ably higher than that of sorghum (Tables 4 and 5). In 
both species, the ash content increased significantly with 
increases in salinity level of the irrigation water. The 
K/Na ratio of the water did not significantly affect the 
ash content of Atriplex, at any salinity level, but in sor- 
ghum, at the higher salinity levels (100 and 150 mM), ash 
content was significantly decreased as the ratio of K to 
Na in the irrigation water declined. 


Cation and nitrate concentration 
Cation and nitrate concentrations in the leaf tissue of Atri- 


plex are from Experiment I only (Table 4). The data from 
Experiment II were very similar. As expected from the 


Treatments K/Na Water potential (MPa) 





Atriplex Sorghum 





Salinity Experiment Experiment 
{mM (K+ Na)] I II 





4 (control) —1.4 —1.4 
—1.3 —1.2 
—1.3 —1.5 





L.S.D. (0.05) 

For salinity (S) 

For K/Na ratio (R) 
For SxR 





ash contents, the concentrations of Na and K in the leaf 
tissues increased with salinity in both Atriplex (Table 4) 
and sorghum (Table 5). In Atriplex, Na and K increased 
in relation to the increases in salinity when they were 
present in equal concentrations in the irrigation water, 
but when the ratio of K to Na was less than 1, most of the 
increase in salt content was accounted for by Na. This is 
indicated by the K/Na ratios in the tissue (Table 4). The 
K concentration changed very little with increases in 
salinity beyond the 100 mM level at K/Na ratios of 1/10 
and 1/100. 

In Sorghum, in contrast, most of the salinity increase 
in the plant tissue was accounted for by K. Sodium did 
not increase much until the external salinity exceeded 
100 mM (Table 5). In fact, both Na and K increased dra- 
matically at the 150 mM level. Even when the K/Na ratio 
in the irrigation water heavily favored Na, the plants still 
accumulated large quantities of K. As a result, the K/Na 
ratio in the tissue remained high, even at external ratios 
of 1/100. If the dramatic increase in tissue salt concentra- 
tion had been the result of a complete loss of the ability 
of the root membranes to control salt uptake, the ratio of 
K/Na in the tissue might be expected to more closely 
reflect the ratio of those ions in the external solution. 
Since that was not the case, the membranes apparently 
still retained some differential permeability, even at the 
external salt concentration of 150 mM. 





130 


Table 4. Effect of salinity and K/Na ratio on ash and mineral content of Atriplex amnicola 





Salinity treatment K/Na Ash Na K Ca 
%DW 





meq/100 g DW 





(K+ Na) mmol/l 


4 (control) 





L.S.D. (0.05) 

For salinity (S) 

For K/Na ratios (R) 
For SxR 





Table 5. Effect of salinity and K/Na ratio on ash and mineral content of Sorghum bicolor 





Salinity treatment K/Na Ash Na K Ca 
%DV 





meq/100 g DW 





(K+ Na) mmol/l 


4 (control) ; ; 50.8 


S12 
43.3 


Tok, 
63.7 
61.0 


128.7 
70.7 
63.8 





L.S.D. (0.05) 

For salinity (S) 

For K/Na ratios (R) 
ForSxR 





The divalent cation content of the tissue did not salinity levels the nitrate concentration was significantly 
change much with increases in external salinity. In both less at the higher K/Na ratios. In sorghum, however, the 
species, the changes in Ca concentration were not signif- nitrate concentration of the tissue did not change signifi- 
icant (Tables 4 and 5). The changes in Mg concentration cantly with either salinity or K/Na ratio (Table 5). 
were not significant in sorghum, but in Atriplex the Mg 
concentration did decline with salinity. At the higher 
salinities there also was a significant effect of K/Na ratio Soluble carbohydrate 
in the irrigation water. When K represented half of the 
salinity, the Mg concentration was significantly lower Soluble carbohydrate content of plant tissue increased 
than when K was only 1/10 or 1/100 of the total salinity. much more in sorghum (Table 7) than in Atriplex (Table 6) 

The nitrate concentration of the Atriplex tissue de- in response to salinity. The K/Na ratio of the irrigation 
clined with increasing salinity (Table 4), and within all — water did not have much influence in Atriplex, the differ- 





Table 6. Effect of salinity and K/Na ratio on soluble carbohydrates, 
free amino acids and proline contents in leaves of Atriplex amnicola 





Salinity 
treatment 


K/Na Soluble Free 


carbohydrate amino acids 


Proline 


Glucose 
equiv. 


Glycine 
equiv. 





re oi I II 





(K + Na) 
mmol/l 


4(control) 1/1 Za ¥: 6.5 
1/10 22. & 6.5 
1/100 26 & 7.4 


1/1 2¢ 3 6.5 
1/10 23 2 73 
1/100 2.6 3. 6.7 


1/1 27 3 7A 
1/10 2.7 3. 7.2 
1/100 28 3. 73 


1/1 29 3. 6.9 
1/10 2.6 3. 9.5 
1/100 2.6 3. 8.7 


umol/g fresh weight 





L.S.D. (0.05) 

For salinity (S) 

For K/Na ratios (R) 
ForSxR 





* Experiment I 
** Experiment II 


Table 7. Effect of salinity and K/Na ratio on soluble carbohydrates, 
free amino acids and proline contents in leaves of Sorghum bicolor 





Soluble Free 
carbo- amino 
hydrate acids 


Salinity 
treatment 


K/Na Proline 


Glucose 
equiv. 


Glycine 
equiv. 





(K+ Na) mmol/l umol/g fresh weight 


4 (control) 1/1 & 7) Ta 
1/10 3.8 6.4 
1/100 52 Te 


1/1 3.7 15.6 
1/10 5.5 12.4 
1/100 5.7 14.9 


1/1 48 18.2 
1/10 6.1 18.3 
1/100 7.8 16.8 


1/1 4.9 22.2 
1/10 6.5 20.4 
1/100 8.9 





L.S.D. (0.05) 

For salinity (S) 0.3 
For K/Na ratios (R) 0.3 
For SxR 0.6 





* Experiment I 
** Experiment II 
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ences being non-significant in the second experiment and 
small, albeit significant, in the first experiment. On the 
other hand, soluble carbohydrate content of the tissue in 
sorghum increased considerably as the K/Na ratio of the 
irrigation water decreased at all salinity levels beyond the 
control level. 


Amino acids 


In general, the amino acid content of Atriplex tissue was 
not affected much by the external salinity level nor by the 
K/Na ratio of the irrigation water (Table 6). The proline 
level was increased slightly by the intermediate salinity 
levels, and it was considerably higher at the highest salin- 
ity level. At that high salinity level, the proline content 
was highest at the highest K/Na ratio, but at the other 
salinity levels the proline content was not affected. 

In sorghum, the proline content steadily increased 
with salinity. At the two higher salinity levels, the proline 
content was significantly higher when the K/Na ratio was 
higher, similar to the situation in Atriplex. 


Discussion 


As is characteristic for the two types of plants represent- 
ed, fresh and dry weights of the glycophyte, sorghum, 
decreased with each increment of added salinity, but were 
not affected in the extreme halophyte, Atriplex, until 
salinity exceeded 100 mM. For some reason or other, 
however, the Atriplex did not have significantly higher 
growth at 100 mM salinity than at the relatively non- 
saline level of 4 mM, a response that is not in agreement 
with the results of Aslam et al. (1986) for this same spe- 
cies. In both species, the growth reduction due to salinity 
was increased by increasing the K/Na ratio of the irriga- 
tion water, consistent with results obtained by others for 
these and similar species (Flowers et al. 1977; Weimberg 
et al. 1984). This response probably was due not only to 
the growth inhibiting effects of high levels of K (Munns 
et al. 1983; Salim and Pitman 1983; Weimberg et al. 1984) 
but also to the growth stimulating effects of Na, in halo- 
phytes (Brownell and Crossland 1972; Flowers et al. 
1977; Eshel 1985) as well as glycophytes (Williams 1960; 
Brownell 1967; Devitt et al. 1984; Rathert 1983). 

There were no apparent signs that either of the species 
suffered from water stress under any of the treatments. 
Leaf water potentials declined about as expected with 
increases in salinity in both species. The only salinity level 
that was common to both species was 100 mM, and that 
caused a drop in water potential in sorghum of 0.2 to 
0.3 MPa and in Atriplex of 0.5 to 0.7 MPa. The 100 mM 
increase in K + Na of the irrigation water is equivalent to 
a drop in osmotic potential of about 0.45 MPa, so in one 
species the leaf water potential decrease was slightly less 
than the drop in external water potential and in the other 
slightly greater. Atriplex had a leaf water potential that 
was 0.2 to 0.4 MPa lower than that of sorghum even in 
the non-saline treatment. At 100 mM, the difference be- 
tween the two species was 0.4 to 0.8 MPa. That may 





132 


be a reflection that sorghum is a relative excluder, and 
Atriplex a relative accumulator, of Na and Cl when pres- 
ent in high concentration in the root environment. This 
response apparently was a response to total salinity, since 
there were no significant differences in leaf water poten- 
tial at different K/Na ratios at any salinity. 

The relative difference in exclusion vs accumulation, 
alluded to above, is readily apparent when the tissue ion 
contents are considered. In Atriplex, the ash content of 
the leaves was high even under non-saline conditions 
(4 mM), and increased substantially with increasing salini- 
ty (Table 4). In contrast, the ash content of the leaves of 
sorghum was relatively low at the lowest salinity (4 mM), 
and it increased very little as salinity increased to 50 and 
100 mM. At 150 mM external salinity, the tissue ash 
content increased dramatically, especially when K made 
up half of the salinity (Table 5). That clearly reflects the 
marked uptake of K in those plants. In Atriplex, both Na 
and K contributed equally to the increase in Na+ K con- 
tent of leaves when they were present in equal amounts in 
the irrigation water, and as the proportion of K in the 
irrigation water decreased, the relative contribution of K 
to leaf Na+ K content decreased. In contrast, in sorghum 
most of the Na+K content was accounted for by in- 
crease in K regardless of the K/Na ratio in the irrigation 
water. This has been reported to be the case in numerous 
Graminaceous halophytes (Albert and Popp 1977; Gor- 
ham et al. 1980), and these results for sorghum indicate 
that it also could be true for Graminaceous glycophytes. 

The total amount of K+ Na accumulated in leaves 
was affected differently in the two species by the ratio of 
those two ions in the irrigation water. In Atriplex, the 
amount of K + Na in the leaves increased with decreasing 
proportion of K in the irrigation water, but in sorghum, 
the amount of K + Na tended to decrease with decreasing 
proportion of K in the irrigation water (Tables 4 and 5). 
Lauter et al. (1988) found increasing levels of K + Na in 
leaves of cotton as the proportion of K decreased in solu- 
tions with varying proportions of K and Na, which is 
similar to the results with Atriplex in this study. This may 
be related to cotton’s higher tolerance of salinity than 
most crops. 

Even though nitrate has been reported to play a role in 
osmoregulation in other species (Smirnoff and Stewart 
1985; Veen and Kleinendorst 1985; Gorham et al. 1986), 
it did not seem to be involved in the species studied here. 
It either did not change significantly with salinity, as in 
sorghum (Table 5), or it decreased with salinity, as in 
Atriplex (Table 4). Such a decrease has been observed in 
other species as well (Abdel-Kadir and Paulsen 1982; 
Aslam et al. 1984). It is interesting to note that the salin- 
ity-induced decrease in nitrate content was enhanced by 
increasing the proportion of K in the irrigation water, the 
significance of which is not immediately obvious. 

The lower amount of salt accumulation in sorghum in 
response to salinity was compensated for by greater ac- 
cumulation of soluble organic compounds, which pre- 
sumably were osmotically active and could contribute 
to osmotic adjustment. The increase in amino acids, and 
specifically proline, was much greater than the increase in 
soluble carbohydrates (Table 7). In Atriplex, the situation 


was similar, except that the increase in amino acids, in- 
cluding proline, was not nearly as great as in sorghum 
(Table 6). 

In summary, the hypothesis posed at the outset of this 
study was not supported by the data. In spite of the 
ability of the halophyte, Atriplex, to accumulate and tol- 
erate high levels of Na, K uptake still was high enough to 
cause significant growth reduction even when provided 
with equal concentrations of Na and K in the environ- 
ment. For example, in sorghum at a K/Na ratio of 1/100, 
the DW at the highest salinity (150 mM) was 33% of that 
at the non-saline control, but when the K/Na ratio was 
1/1, the DW at the highest salinity was only 24% of that 
at the non-saline control. In Atriplex, when the K/Na 
ratio was 1/100, the DW at the highest salinity (400 mM) 
was 52% in Experiment I and 60% in Experiment II of 
that at the non-saline control. However, when the K/Na 
ratio was 1/1, the DW at the highest salinity was only 
45% and 46% of that at the non-saline control. Thus, 
even though the absolute amount of growth reduction 
due to the total salinity and to the relative amounts of Na 
and K present in the external solution was not the same 
in the two types of plants, it is apparent that they respond 
similarly to those environmental stresses. 
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Abstract. Direct measurements of the xylem sap flow by 
the stem heat balance technique can be a valuable aid for 
determining the irrigation demand of field crops. In the 
present study, soybean (Glycine max (L.) Merr.) sap flow 
was evaluated under well-watered and water-stressed 
conditions using Dynamax SGA10 sap flow gauges. So- 
lar radiation was measured continuously throughout the 
growing season. Soil water content was measured before 
and after each irrigation. There was a close relationship 
between solar radiation and xylem sap flow. The water 
flux in the soybean stems responded realistically to 
changes in the soil water content. However, the absolute 
values of sap flow were highly questionable. Calculating 
crop transpiration from sap flow measurements, the re- 
sults were up to 4 times as high as calculated transpira- 
tion from soil moisture data and simulated transpiration 
using the locally calibrated soybean crop growth model 
SOYGRO. A sensitivity analysis of the stem heat balance 
technique gave no indications of technique or input er- 
rors. The gauge design was possibly not appropriate for 
the outdoor installation on soybeans. 





The water flux in the stem of intact plants can be mea- 
sured directly using the stem heat balance (SHB) tech- 
nique according to Sakuratani (1981). Since the xylem 
sap flow is closely related to transpirational water losses 
(about 99% of plant water uptake is lost by transpira- 
tion), sap flow is assumed to equal transpiration. The 
relative accuracy of the SHB method compared with di- 
rect weighing of water losses is reported to be about 
+10% (Sakuratani 1981, 1984, 1987; Baker and Van 
Bavel 1987; Steinberg et al. 1990). However, most of the 
studies using the SHB technique were conducted in con- 
trolled environments. The objective of this study was to 
determine soybean crop transpiration by sap flow mea- 
surements under well-watered and water-stressed condi- 
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tions in the field and to compare the results with soil 
moisture measurements and simulated transpiration by 
the locally calibrated crop growth simulation model 
SOYGRO V5.42 (Jones et al. 1989). 


Materials and methods 


The research was conducted in collaboration with the Technical 
University of Lisbon (Instituto Superior de Agronomia, Departa- 
mento de Engenharia Rural) at the irrigation station ‘Antonio Teix- 
eira’ in Monte da Barca near Coruche (38° 57’ N, 8° 32’ W). 
One hectare of soybeans Glycine max (L.) Merr. cv. Kingsoy, 
were sown on May 17, day of year (DOY) 137, in 1991 on a sandy 
soil with a sand portion of 88.3% to a soil depth of 60 cm. Row 
spacing was 50 cm, plant spacing 6 cm, and plant density 32 m~?. 
Water was supplied by sprinkler irrigation. There were two treat- 
ments, (1) well irrigated (Control) and (2) water-stressed after 
bloom (Stress). The Stress treatment covered 0.25 ha and was not 
irrigated for 13 days (DOY 208-220). 

Water flux in the stem of intact soybean plants was measured 
directly using the stem heat balance (SHB) technique after Saku- 
ratani (1981). With the SHB method, the thermal energy balance of 
a stem section which is continuously heated over a short vertical 
distance is determined. The distribution of radial and vertical heat 
is measured and the convective heat transport by the sap stream can 
be calculated. The energy balance of a heated stem segment is 
expressed as: 


Q,=Q,+Q,+Q +8. (1) 


Q,, — heat flux applied to the stem (W), Q, — radial heat flux away 
from the stem (W), Q, — vertical heat flux through the stem (W), Q, 
— heat flux carried by the xylem sap (W), S — rate of heat storage 
(W). 

Heat storage (S) in stems of herbaceous plants is negligible 
because of the small size of the heated stem segment and near steady 
state conditions in the system (Ham and Heilman 1990). At low 
flow rates, however, the heat storage term should be included in the 
calculations (Groot and King 1992). 

The energy input, Q,, follows Ohm’s law 


Q,=U?/R. 
U — input voltage (V), R — heater resistance (2). 


The vertical heat conduction, Q,, can be up and down the 
stem: 


Q,=Q,4+Q,- (3) 
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Q,,, Qy — axial heat conduction up and down the stem (W). 
Fourier’s law for a one-dimensional heat flow describes the axial 
heat conduction: 


Q,=K,-A- dT,/dx (4) 
Q,=K,: A> dT,/dx (5) 


K, — thermal conductivity of the stem tissue (Wm~'K~'), A - 
cross sectional area of the stem (m2), dT,,, dT, — temperature above 
and below the heater (K), dx — spacing between thermocouple junc- 
tions (m). 

The radial heat loss, Q., is computed by 


Q,=K,:E. (6) 


K, — thermal conductance constant for a gauge installation 
(W V~*), E — radial temperature gradient (V). 

Solving equation (1) for Q,, the sap flow rate can be calculated 
as described by Sakuratani (1981) and Baker and Van Bavel (1987): 


F=(Q,-Q,-Q)/(C, : dT). (7) 


C, — specific heat of water (Ws g-' K~'), dT — temperature 
increase of the sap (K). 

The sensor (sap flow gauge) consisted of a thin resistance heater 
completely surrounding the stem and two pairs of thermojunctions 
which measured the temperature of the stem surface above and 
below the heater (dT,, dT, in Eqs. (4) and (5)). An eight junction 
thermopile, wired across the cork backing of the heater, measured 
radial heat loss (E in Eq. (6)). 

In this study, eight of the improved four-channel-gauges (Stein- 
berg et al. 1990) with weather shield kits were used (model SGA10, 
Dynamax, Houston, USA). To prevent wetting and improve ther- 
mal contact with the stem, an electrical insulating compound (mod- 
el G4, Dow corning, Midland, USA) — was applied to the gauges 
(Van Bavel and Van Bavel 1990). Four sap flow gauges were at- 
tached to soybean stems (9 mm diameter) in the Control plot and 
another four in the Stress treatment. The input voltage (4—5 V) to 
the resistance heater of the gauges was supplied by an adjustable 
voltage divider and a datalogger-multiplexer unit (model 21X/ 
AM32, Campbell Scientific, Leicestershire, UK) continuously 
recorded the electrical signals. Solar radiation was also recorded by 
the datalogger using a pyranometer (model SKS 1110, Skye Instru- 
ments, Llandrindod Wells, UK) at a height of 2 m above the soil 
surface. Energy was supplied by a 12 V battery which was constant- 
ly recharged by a 35 W solar panel (model M65, Arco Solar, 
Chatsworth, USA). An execution interval of 30 s was chosen for all 
sensors and every 30 min the averaged values were stored. 

Since maximum rooting depth of soybeans was found to be 
60 cm, soil water content was measured at six depths (0—60 cm) in 
10cm increments using a neutron probe (Didcot Instrument, 
Abingdon, UK). In every treatment three access tubes were used 
and at each depth three replications were measured. The moisture 
content of the surface layer (0Q—10 cm) was determined gravimetri- 
cally and multiplied with the dry bulk density to get the volumetric 
water content. Measurements were made one day before and after 
each irrigation. During stress periods, soil water was determined 
every second day. 

For comparing the sap flow results with simulated crop transpi- 
ration, the soybean model SOYGRO V5.42 (Jones et al. 1989) was 
used. The model was calibrated at the experimental site in 1990 and 
1991 for the cultivar Kingsoy using local weather, soil profile and 
management data (Gerdes 1993). For that purpose, leaf area index 
(LAI), dry matter accumulation and root length density were mea- 
sured during the growing seasons. In SOYGRO, soil water balance, 
root growth and root water uptake are calculated using the model 
developed by Ritchie (1985). Saturated upper limit, drained upper 
limit and lower limit of plant available water for each soil layer were 
derived from laboratory measured pF-curves. Actual plant transpi- 
ration in the SOYGRO model depends on the potential evapotran- 
spiration (E,), calculated after Priestley and Taylor (1972), and the 
LAI. For LAI <3.0, transpiration is E, multiplied by a LAI depen- 
dent factor and when LAI >3.0, plant transpiration is assumed to 


equal E,. The rate of root water absorption depends on plant 
available soil water, root length density and soil water conductivity 
coefficients as defined by the soil profile input data. When the total 
root water uptake is less than actual plant transpiration in the 
model, a water stress factor reduces plant growth (Ritchie 1985). 


Results and discussion 


Diurnal course of sap flow and solar radiation in the well 
irrigated treatment is shown in Fig.1. The sap flow 
curves represent two soybean plants measured with dif- 
ferent gauges. There was a close relationship between sap 
flow and solar radiation as would be expected since solar 
radiation influences stomatal aperture, air temperature 
and vapour pressure deficit. These weather and plant 
factors affect directly and indirectly the rate of transpira- 
tion. Measurable sap flow began at 6:30h and ceased 
between 20:00 h and 21:00 h. In all replications, the max- 
imum sap flow was recorded between 13:00h and 
15:00 h and ranged from 55 to 95 gh~'. The total sap 
flow for DOY 206 (Fig. 1a) was 583.6 g (ml) and for 
DOY 210 (Fig. 1b) 512.3 g. 

Daily accumulated sap flow per soybean plant in the 
well-watered and water-stressed treatments together with 
the corresponding solar radiation data are shown in 
Fig. 2. Though standard deviations of sap flow were con- 
siderable for some days, there were pronounced differ- 
ences between the treatments. In the water-stressed plot, 
the sap flow rate declined after DOY 213, seven days 
after the last irrigation. The gradual reduction of sap flow 
in both treatments after DOY 226 was due to the ad- 
vanced growth stage of soybeans (decreasing LAI). 
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Fig. 1a, b. Diurnal course of soybean sap flow and solar radiation 
on (a) DOY 206 and (b) DOY 210 in the well watered treatments 
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Fig. 2. Daily accumulated sap flow (F,,) per soybean plant in the 
Control and Stress treatment (irrigated DOY 207 and 221, 1991) 


with corresponding solar radiation. The bars indicate standard de- 
viations 


There was a reasonable relationship between daily solar 
radiation and sap flow. 

Soybean crop transpiration was calculated in this 
study by multiplying the daily accumulated sap flow per 
plant by a population of 32 plants m~?. There were no 
major differences in leaf area and stem diameter among 
the soybean plants in the measurement area. The crop 
growth model SOYGRO was used to compute soybean 
transpiration based on weather conditions and simulated 
LAI. It is assumed that the Priestley-Taylor equation is 
appropriate to calculate E, for Central Portugal. The 
calibrated model satisfactorily predicted LAI, soil water 
content and water stress (Gerdes 1993). Therefore, the 
simulated SOYGRO results for crop transpiration are 
considered as reasonable reference values. However, a 
long-term validation of the calibrated model has not yet 
been carried out. In Table 1, the soybean crop transpira- 
tion derived from sap flow measurements is compared 
with the simulation results for crop transpiration ob- 
tained by the locally calibrated SOYGRO model. In the 
Control treatment, transpiration calculated using the 
SHB method on an average was 3.3 times as high as the 
simulated transpiration. In the Stress treatment, mea- 
sured and simulated transpiration decreased uniformly. 
However, crop transpiration based on sap flow measure- 
ments was still about 4 times as high as the SOYGRO 
simulation results. 

Water extraction according to soil moisture measure- 
ments (neutron probe) is shown in Table 2. The mean ET 
(evapotranspiration) values were calculated by determin- 
ing the decrease of water content in the effective soil 
profile (0—60 cm depth) in the period between irrigations. 
Leaching of soil water was negligible and after DOY 200 
(nine weeks after sowing), the soybean canopy totally 
covered the soil so that ET approximated actual plant 
transpiration. The mean ET values for both treatments 
were in good agreement with the simulated transpiration 
in Table 1. A transpiration rate of 4 to 8 mm d™~! at that 
growth stage (pod elongation) under non-water stress 
conditions appears reasonable (Doss and Thurlow 1974). 
The standard deviations for ET in Table 2 indicate that 
the transpiration rate measured by the Dynamax SGA10 
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Table 1. Soybean crop transpiration (T) measured using the stem 
heat balance (SHB) method and simulated using the crop growth 
model SOYGRO for the well irrigated (Control) and water-stressed 
(Stress) treatments in 1991 





Day of 


year 


Measured (SHB) Simulated (SOYGRO) 





Stress* 
T(mmd~') 


Control 
T(mmd~') 


Stress* 
T(mmd~') 


Control 
T(mmd~!) 





203 dee, 5.4 5.4 
204 16.4 16.3 5.3 5.3 
205 19.0 19.1 5.7 5.7 
206 18.7 18.6 39 $9 
207 14.4 15.3 5.9 59 
208 21.3 21.0 5.9 5.9 
209 47.1 the x 5.3 
210 16.4 fe 5 4.5 4.5 
211 13.0 13.6 4.2 4.2 
212 18.1 17.8 5.4 5.4 
213 22k 17.5 6.0 4.2 
214 21.3 11.4 6.3 1.9 
i 23.3 9.0 6.2 4.3 
216 20.1 5.1 5.9 1.0 
217 17.6 3.2 5.0 0.8 
218 15.6 2.4 5.7 0.6 
219 17.6 PIR $2 0.5 
220 20.0 1.9 5.8 0.5 





* Water deficit period from day of year 208-220 


Table 2. Consumptive use (ET) of the soybean crop calculated from 
changes in the soil water content (0—60 cm depth) between irriga- 
tions for the well-watered (Control) and the water-stressed (Stress) 
treatments in 1991 





Treat- 
ment 


Period Num-_ Extracted Mean Standard 
(day of _berof soil water ET deviation 
year) days (cm*cm~*) (mmd~') (mmd“‘) 





171-177 
180-183 
186-189 
213 

216-217 
211-212 


213-216 
217-220 


Control 0.035 3.0 
0.028 4.2 
0.029 4.4 
0.013 7.8 


0.021 6.3 


0.024 Lp: 
0.019 29 
0.006 0.9 


HRN NRF HHA 





* Water deficit period from day of year 208-220 


sap flow gauges (Table 1) was unrealistic. Considering 
the standard deviations for the sap flow measurements 
(Fig. 2) for DOY 203 to 220, the transpiration in both 
treatments was still overrated. 

The calculation of the consumptive use of the soybean 
crop and the simulated transpiration results indicated 
that the flow gauges used in this study overestimated 
plant water loss under the prevailing conditions. To ex- 
amine possible sources of error which can arise from the 
stem heat balance method, a sensitivity analysis was 
made (Table 3). The SHB method was most sensitive to 
changes in the heater resistance (R in Eq (2)) of the gauge. 
An increase in R by 10% resulted in a 17.4% decrease in 
sap flow. However, the sensitivity analysis demonstrated 
that even large errors in a particular input parameter of 
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Table 3. Sensitivity analysis for the sap flow measurements with 
respect to modified input parameters. The analysis is based on 
results from a Dynamax SGA10 gauge installation on day of year 
210 in 1991 





Modified P-change 
parameter (P) (%) 


Sap flow (F) F-change 
(gd~') (%) 





K, (WmV~!)* 
0.745 0.0 
0.820 +10.0 
0.671 —10.0 


Stem area (cm?) 
0.636 0.0 
0.785 +23.4 543.1 —3.9 
0.950 +49.4 526.5 —6.8 
1.131 +77.8 510.0 —9.7 
Heater resistance (Q) 
112.0 0.0 
123.2 +10.0 
100.8 —10.0 
K, (Wm! K~!)° 
0.544 0.0 
0.59 +10.0 
0.49 —10.0 
0.42° —22.2 


565.0 0.0 
518.7 —8.2 
777.0 +37.5° 


565.0 0.0 


565.0 0.0 
466.5 —17.4 
894.2 +58.3° 


565.0 0.0 
554.7 —1.8 
580.7 +2.8 
616.7 +9.2 





*K, — thermal conductance constant for a particular gauge installa- 
tion, ° unrealistic because of calculated sap flow also at night, 
© K, —- thermal conductivity of the stem, * herbaceous plants (Saku- 
ratani 1984), *° woody plants (Sakuratani 1984) 


+10% or more, which are unlikely, did not modify the 
calculated sap flow to more acceptable values. 

Shackel et al. (1992) used the SHB method to measure 
transpirational sap flow on three-year-old peach trees 
under field conditions. They reported that in some cases, 
errors in sap flow approached two orders of magnitude as 
a result of very small and sometimes negative tempera- 
ture differentials (dT) between the stem surface tempera- 
ture above and below the gauge heater. These authors 
suggested that the effect of ambient conditions on gauge 
signals should be critically evaluated before considering 
heat balance estimates of sap flow as reliable under any 
given conditions. Baker and Nieber (1989) also recog- 
nized that dT is critical in the calculation of sap flow. 

A very low dT may cause overestimations of sap flow 
as dT is part of the denominator in Eq. (7). However, 
Fig. 3 demonstrates that dT during the measurement pe- 
riod of this study was generally in the recommended 
range of 0.2°C to 8°C. In the morning, dT should not 
exceed 8°C and during the peak flow period dT should be 
above 0.2°C (Van Bavel and Van Bavel 1990). 

The course of the stem surface temperature dT versus 
sap flow rate is shown in Fig. 4 for two gauges on differ- 
ent soybean plants. As would be expected, dT declined 
with increasing heat transfer by convection. To avoid 
ambient effects on the temperature regime of the gauges, 
a weather shield kit was used. Moreover, gauges were 
well-shaded by the soybean canopy, so that an error in dT 
seems unlikely for the sap flow determinations reported in 
this research. Figure 5 shows the heat balance of a SGA10 
gauge installation on a soybean stem. The distribution of 


@-© July 25, 1991 (gouge 1) 
4 4-4 July 29, 1991 (gauge 2) 
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Fig. 3. Diurnal variations of the stem surface temperature dT 
above and below the heater. Data from two gauge installations 
(model SGA10) on different soybean plants 
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Fig. 4. Temperature difference dT at the stem surface versus sap 


flow rate. Data from two gauge installations (model SGA10) on 
different soybean plants 
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Fig. 5. Diurnal course of heat distribution on a soybean stem. Data 
from a SGA10 sap flow gauge on August 15, 1991 


the heat supplied was reasonable and in good agreement 
with the energy balance of a gauge installation on Red 
Cedar described by Van Bavel and Van Bavel (1990). 
Another possible reason for the overestimation of 
crop transpiration by sap flow measurements was given 
by Ham et al. (1990). They reported that their experimen- 
tal plot included many cotton plants that were much 
smaller than those used for sap flow measurements. Con- 
verting the sap flow per plant to crop transpiration by 
multiplying by the plant density, resulted in a 63% larger 





value than the value obtained by normalizing sap flow 
measurements by leaf area. Ham et al. (1990) suggested 
that in the case of a variable plant size within the canopy, 
the sap flow per plant and day should be divided by the 
leaf area per plant and then multiplied by the mean leaf 
area index of the plot. However, in this study the soybean 
canopy in the sap flow measurement area was quite ho- 
mogeneous and no pronounced differences in plant size 
were detected. Moreover, the relative good agreement 
among the replications of the sap flow measurements 
indicated that variable plant sizes were not the reason for 
the overestimation of soybean crop transpiration by the 
SGA10 flow gauges in the present work. 

The sensors used in this study were designed for herba- 
ceous plants with a stem diameter of 9 to 12 mm. The 
width of the resistance heater was 10 mm and the energy 
input for all gauges ranged from 0.10 to 0.25 W. Ham and 
Heilman (1990) tested similar sap flow gauges on sun- 
flower and found that transpiration was constantly over- 
estimated during periods of high sap flow. Adjusting the 
power to the stem did not improve gauge performance at 
high flow rates. However, increasing the heater width 
greatly improved gauge performance. The relationship 
beween stem diameter and size of the heated segment is 
important for the gauge design (Ham and Heilman 1990). 
According to Sakuratani (1981), the heater width may be 
two times larger than the diameter of the stem. Possibly, 
the SGA10 sap flow gauges in the present work overesti- 
mated soybean transpiration because the heater was too 
narrow (equal to the stem diameter) and the power distri- 
bution of the gauges may not have been suitable for a 
outdoor installation on soybeans. 


Conclusions 


The stem heat balance technique is a non-destructive 
method, suitable for measuring the xylem sap flow of 
intact herbaceous and woody plants. Sap flow gauges 
responded quickly to changes in the environmental con- 
ditions (solar radiation, water supply). However, soy- 
bean crop transpiration calculated from field-measured 
sap flow data differed markedly from crop water use 
determined by changes in the soil water content and sim- 
ulated transpiration by the locally calibrated SOYGRO 
model. It is possible that the heater width of the Dyna- 
max SGA10 gauges was not appropriate for measuring 
the soybean transpiration rate in the field. Currently, the 
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stem heat balance method has not often been examined 
under field conditions. Further investigations are re- 
quired to detect the critical operation range of this new 
technique. 
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Abstract. This paper evaluates the sensitivity of two dif- 
ferent methods of water stress detection in a simulated 
patch of pasture grown in a greenhouse. The perfor- 
mance of two indices, based respectively on canopy tem- 
perature and soil water content values — the latter gauged 
by means of a time domain reflectometry (TDR) system, 
was assessed against actual evapotranspiration, mea- 
sured by a very accurate weighing system. Both methods 
were able to detect water shortage by the time transpira- 
tion was reduced to some 80% of its potential value. The 
soil-based index, however, relied on the estimate of root 
water extraction rate, which may not be known. It is 
concluded that detection of water shortage by means of 
a canopy temperature-based stress index is to be pre- 
ferred to measuring soil water deficits by time domain 
reflectometry, despite the accuracy of the TDR-based soil 
water content estimate. 





Introduction 


Although there is little doubt that maximum yields are 
only achieved under conditions of non-limiting water 
supply, the question how early a water shortage is to be 
detected in order to avoid a significant decline in produc- 
tivity is still an open one. As far as pasture is concerned, 
Krist and Snijders (1987), summarizing several years of 
experiments with grassland in The Netherlands, conclud- 
ed that yield reduction per mm potential moisture short- 
age could vary from some 18 kg dry matter ha~' on clay 
soil, to more than 30 kg ha~' on sandy soil. It would 
seem, therefore, that optimum water management, at 
least for pasture, requires water stress detection to be as 
timely as is technically feasible. 

As the research described here was a part of the devel- 
opment of a farm management information system, we 
confined ourselves to testing stress detection by means 
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suitable for remote operation. For many years soil water 
content has been measured by either manual or inaccu- 
rate methods, so that soil based stress detection methods 
did not provide a feasible alternative to the canopy-tem- 
perature based stress indices now widely used. Time-do- 
main reflectometry, TDR, (Topp et al. 1980), however, 
has proven to be an easily automated and accurate means 
for measuring soil water content, so that a TDR-based 
remote stress detection is presently feasible. 

The principle underlying canopy temperature-based 
stress detection is simple: plants not adequately supplied 
with water will transpire less than well-watered plants, 
and hence will be warmer than their well-watered coun- 
terparts. Detection of water stress by this means, howev- 
er, implies both that the temperature of an otherwise 
similar, well-watered plot be known and that the devia- 
tion of temperature from this “ideal” value can be detect- 
ed with sufficient accuracy, two conditions not easily 
met. ; 

Soil-based detection of water shortage, on the other 
hand, is based on the observation that it is increasingly 
difficult for a crop to extract water from the soil beyond 
a threshold pressure head. However, recognition of water 
shortage by this means requires that this threshold value 
be known, which is not trivial either. 

Accordingly, both methods have potential shortcom- 
ings, so that it was thought worthwhile to test their accu- 
racy in the present study. 


Temperature-based stress index 


From the coupling of the energy balance equation to the 
transfer equations for both sensible and latent heat, an 
expression for the temperature difference between the 
canopy and the air can be derived (e.g. Monteith 1981): 





K (1) 
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where T is temperature (K) and the subscripts c and a 
refer to the canopy and air, respectively; R,, is the net flux 
density of radiation (Wm~ 7); G is the flux density of 
sensible heat transferred into the soil (W m~ 7); D is va- 
pour pressure deficit of air (Pa); r, and r, (s m~ ') are the 
boundary layer (aerodynamic) resistance and the canopy 
resistance, respectively. Furthermore, @,c, is the volumet- 
ric heat capacity of air (Jm~* K~'); € is the rate of 
change of latent heat content with sensible heat content 
of air saturated at present temperature and y is the ther- 
modynamic psychrometric constant (Pa K~'). 

Jackson (1982), reviewing previous literature, con- 
cluded that the effect of wind on canopy temperature, 
implicit in the boundary layer resistance, is likely to be 
minor (at least for winds not too strong), a contention 
supported for pasture under cloudy humid conditions by 
the experimental results of Feldhake and Edwards 
(1992). Hence, temperature of a vegetated canopy is af- 
fected primarily by three climate factors: air humidity, 
net radiation balance and air temperature; if, and only if, 
a high correlation exists among radiation, temperature 
and humidity, the number of climate variables affecting 
canopy temperature could be reduced, possibly, to one — 
as proposed by Idso (1982) with his non-water-stressed 
baseline approach. In variable climates the influence of 
individual weather factors needs to be taken into ac- 
count, as concluded by Smith et al. (1986), for instance, 
from an experiment in New South Wales. Fuchs (1990) 
confirmed “that the detection of water stress by a canopy 
temperature measurement must account for the effects of 
meterological factors on the energy balance of the fo- 
liage”’. 

A method satisfying this requirement was proposed by 
Jackson et al. (1981), who introduced a crop water stress 
index (CWSI) as follows. The highest possible tempera- 
ture a canopy can achieve (7,,) is given by: 

lim T, 


r 
T, = T, + —*(R,—G). 
M r,— 00 ‘toe ) 


K (2) 
There is, on the other hand, also a lower limit (T,, ) to 
canopy temperature, that is achieved for a minimal 
canopy resistance (r,,, ): 


D 
tat tmp Gg) —? 
c “4 


lim T. 
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Then, varying climatic conditions can be accounted for 

by using the difference T,, — T,, as the reference against 

which deviations of the measured canopy temperature 

from the “unstressed” temperature T,, have to be gauged. 

Namely: 


ea 


ST a oe. 
' ~<. 


—& 
Equation (4) can be calculated from Eqs. (1), (2) and (3), 
by eliminating (R, — G), assuming this to be an indepen- 


dent variable, although strictly speaking it is not. The 
crop water stress index can then be written as an explicit 


function of canopy resistance, or of the ratio of actual to 
potential transpiration: 


r.—r E 
WSI = —————_ = 1-—. — (5 
. ” (1+e)r,+r, E, " 


where E is calculated by the Penman-Monteith equation: 





_ eM a, G) + Q,¢, D/y 


L{Q +6)r, +r, ] eee 


L being latent heat of evaporation of water. Potential 
transpiration (E,) is calculated by Eq. (6) with r, = r,,. 

Applying this “‘big-leaf’’ model to detect stress, obvi- 
ously assumes that the canopy temperature that is mea- 
sured is the one temperature defined by such a model — a 
requirement that may be hard to satisfy for many crops 
that are tall or have incomplete cover (Fuchs 1990), but 
one that grass is likely to meet. Moreover, Jackson’s et al. 
(1981) original definition of CWSI referred to the special 
case that r,, is zero, with the minor consequence that the 
stress index of a non-wet canopy, however well watered, 
is larger than zero; the minimal canopy resistance is 
known to be species-specific (KG6rner et al. 1979). At the 
other extreme, 7,, is also a theoretical upper limit to 
canopy temperature that will never be reached. This be- 
cause actual net radiation would become significantly 
smaller than R,, as r, increased, due to thermal radiation 
from the canopy, as Stanghellini (1987) and Paw U and 
Gao (1988) have shown. Moreover, Jackson and Ezra 
(1985) demonstrated that the albedo of leaves could in- 
crease, and thus net radiation decrease, in consequence of 
water stress. Net radiation is also reduced by drooping or 
curling of leaves, a common plant response to water 
stress, also likely to affect r,. The stress index of plants 
wholly deprived of water, therefore, can still be less than 
one, when calculated by means of Eq. (4). However, as 
the main purpose of this application is early detection of 
water stress, this would appear to be of little conse- 
quence. 


Soil-based stress index 


Shortage of water takes place, by definition, when the 
amount of water extracted from the soil does not match 
potential transpiration. This trivial observation allows a 
soil-based “‘stress index’ to be derived, as follows. Ex- 
traction can be represented (e.g. Feddes et al. 1978) asa 
sink term S, defined by: 


S=a(h) Sy d~' (7) 


where S is water extraction rate cmj,o cm; d~'; ais a 
non-dimensional function of soil water pressure head h, 
contained between 0 and 1; and S,, is the maximum pos- 


sible water extraction by plant roots defined by: 


Su = E,/\z,| "i (8) 


where E, is the potential transpiration rate (cm d~*) and 
z, (cm) the effective rooting depth, defined as the depth 
up to which approximately 75% of the roots are found. 





Obviously, the integral of the sink term over the effective 
rooting depth has to be equal to the actual transpiration 
rate of the crop E (cm d~'): 


en oa. cm d~! (9) 
0 

And, for transpiration to attain its potential value, S 
has to equal S, and «(h) must be unity; on the other 
hand, under non-optimal conditions, S,, is reduced by 
means of the pressure-head-dependent function a. Both 
Sy and the shape of «(h) are commonly assumed not to 
vary throughout the profile; a very useful, although 
crude, simplification which implies that: 


a(h) = E/E, =1— CWSI. — (10) 


Hence, insofar as the function «(h) is known, soil wa- 
ter pressure head is as good a stress indicator as CWSI. 
We will show in the following, that lack of information 
about the function «, and not the accuracy of the soil 
water content assessment, is the factor limiting the reli- 
ability of the soil-based stress indicator, for incipient wa- 
ter shortage. 


Materials and methods 


To assess the suitability of both methods of stress detection, an 
experiment was carried out on a simulated sward of perennial rye- 
grass (Lolium perenne) growing on a weighing lysimeter placed in a 
glasshouse. The glasshouse was single-glass, four span, and E —W 
oriented. The ventilation windows facing north were kept constant- 
ly open, no artificial heating or lighting being provided. The grass 
was grown in a box (0.8 by 0.7 m and 1 m deep) filled with sandy 
soil and completely surrounded by a guard area, where a similar 
sward was grown. The box was connected with a water tank in a 
way that allowed the depth of the ground water table to be detected 
and regulated. The walls of the box and of the tank, as well as the 
pipes connecting them, were insulated in order to reduce heat ex- 
change with the surrounding environment. One face of the box was 
an acrylic plate so that, upon removal of the insulating material, 
rooting depth could be observed. The electronic balance on which 
both the box and the tank were placed, allowed detection of weight 
variations with an accuracy of 2.5 g on a total load of 1 500 kg. The 
grass was mown about once a month, when the mean sward height 
was about 25 cm. 

Climate variables measured were: incoming shortwave radiation 
and net radiation, respectively, by a solarimeter and a net radio- 
meter, both placed about 40 cm above the grass. Air temperature 
and humidity by a thermohygrometer placed at a height of about 
30 cm, with a solid state device for humidity measurement and a 
platinum resistance for temperature measurement. Radiative tem- 
perature of the canopy was determined by an infrared thermometer 
with a field of view of about 2°, placed 70 cm above the crop. All 
measurements were performed at 5-minute intervals by a data-log- 
ger. 

Other parameters estimated were: aerodynamic resistance was 
determined by inverting the heat transfer equation, wherein sensible 
heat flux was calculated through the energy balance equation. This 
procedure is likely to be affected by large errors whenever sensible 
heat flux is small, that is, whenever transpiration and net radiation 
are of similar magnitude, and, indeed, the standard error of the 
estimate of r, approached often 100%. Nevertheless, since meaning- 
ful variations of the boundary layer resistance in a closed environ- 
ment are expected to be small, its mean value of 65 s m~' was used 
throughout. Minimal canopy resistance (r,,) was also taken as 
65sm_~', as proposed for clipped grass by a recent revision of 











T 


: 
15 20 aa 


dielectric constant (K,) 


Fig. 1. Calibration data for the TDR-probe. The thick line repre- 
sents the fitted Eq. (11). The thin line, Topp’s equation for mineral 
soils, with 95% confidence limits (broken lines) 


F.A.O. guidelines (Smith 1991). Heat flux into the soil (G) was 
assumed to be a constant one-tenth of net radiation, as proposed by 
De Bruin and Holtslag (1982); the results proved to be rather inten- 
sitive to this assumption. 

The profile of soil water content in the box was determined 
through a TDR-system, equipped with six three-wires probes, 
placed horizontally at depths of 5, 10, 20, 30, 45 and 60 cm, respec- 
tively. The relationship between measured dielectric constant (K,) 
and volumetric water content of the soil (@,, cm? cm~ *) was deter- 
mined by comparing K, values with gravimetrically determined 0, 
measurements of drying samples of the soil, that was 90.5% sand 
(particle size between 50 and 2 000 pm) and the rest finer material 
(smaller than 50 ym); with 3% organic material. The 275 points 
calibration, shown in Fig. 1, are best-fitted by: 


0, = 3.99 1073 + 1.55 10-2 K, — 1.48 10-5 K? — 3.14 1077 K? 
a 


with a standard error of the estimate of 0, of about 0.5%. Most data 
fall just within the 95% confidence band of Topp’s etal. 1980 
best-fit line for mineral soils. However, Topp et al. (1980) attributed 
most of their variance to differences in soil type, density or salinity, 
that may have played a role here, rather than to experimental errors. 
For instance, both Nadler et al. (1991) and Amato et al. (1993) 
derived calibration curves for mineral soils similar, yet not identical, 
to Topp’s original curve. The relationship between volumetric soil 
water content and pressure head (h), was determined according to 
the standardized Dutch “Staring” series — Wésten et al. (1987): 


\h| = 9.47 10568 — 1.75 10°07 + 1.38 10° 6° — 6.00 10565 + 
1.58 10°0* — 2.57 10463 + 2.52 10°62 — 1.38 1076 + 3.30. 
MPa (12) 


Irrigation took place about once a week, by refilling the amount of 
water losi by evapotranspiration; as the water table was kept con- 
stantly at 72cm depth, A values of the upper soil layer down to 
30 cm, in which the largest part by far of the rooting system was 
found, did not change much: they varied approximately from 
12 kPa soon after an irrigation to 40 kPa before the next one. 
After collecting reference data for about one month, the method 
was tested during a dry period of a couple of weeks, starting from 
May 17. For this purpose both irrigation and ground water table 
control were suspended 15 days after a sward cut. The water table 
dropped slightly, to 72 cm in two weeks. Figure 2 shows the profile 
of 6, the day after the last irrigation and two weeks afterwards. The 
difference in soil water content, calculated as the difference between 
the integral of the two polynomial relationships fitted to the data, 
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Fig. 2. The profile of 0, (volumetric soil water content) the day after 
the last irrigation and two weeks later; points represent data. A 
water content at saturation equal to the highest measured value was 
assumed in order to fit the polynomial relationships shown 


indicated by the shaded area in Fig. 2 was 15.8 kg, whereas decrease 
in weight during the same interval was 15.9 kg. Such an agreement 
has to be fortuitous, given the inaccuracy of the estimate of soil 
water content at saturation and of the profile demarcations. With 
respect to the latter, linear interpolation between data points in 
Fig. 2 would have resulted in an estimated water content decrease 
of 14.7 kg which, nonetheless, sanctions the reliability of the present 
measurements of volumetric soil water content. 


Results 
Temperature-based water stress detection 


Jackson etal. (1981) calculated the CWSI from eight 
measurements around 2 p.m. More data, however, could 
smooth out the effect of intermittent cloudiness on 
canopy resistance, as the latter is assumed to be constant. 
Accordingly, one hundred CWSI values were calculated 
daily from as many measurements centred at 1 p.m. 
(M.E.T.). A reference CWSI set was then assembled by 
averaging concurrent data collected during the 32 days 
preceding, and the three following, the last irrigation. A 
set of one hundred CWSI values was calculated similarly 
for each day of the subsequent dry period and a two-sam- 
ple t-test was performed to detect daily significant devia- 
tions from the reference sample, as shown in Fig. 3. 

Mean daily ratio of actual to potential transpiration 
during the reference period was 0.97, with a standard 
deviation, o = 0.12. Mean CWSI, on the other hand, was 
0.07 with o = 0.06. Figure 4 depicts mean daily value of 
CWSI, for each day of the dry period as a function of 
reduction in transpiration, the solid symbols referring to 
CWSI sets significantly different (confidence level 1%o) 
from the reference set. Means of the above mentioned 
reference set are shown by a larger symbol with standard 
deviation bars. 

If the admittedly arbitrary criterion had been adopted 
that water shortage is proven by two consecutive samples 
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Fig. 3. Mean daily trend of CWSI for the reference period (thick 
line) and the plus and minus one-standard-deviation bands (broken 
lines). CWSI sets of two days are shown as examples: one, May 19, 
open circles, not significantly different from the reference and an- 
other, May 24, closed circles, significantly higher (confidence level 
1%o) 
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1-E/E, 
Fig. 4. Mean daily CWSI as related to the reduction in transpira- 
tion during the dry period. The larger symbol represents the mean 
values for the reference period, with standard deviation bars. Con- 
secutive days are connected, starting from May 17, lower left. Filled 
symbols refer to days whose CWSI set deviated significantly from 
the reference sample, as explained in the text 


of CWSI being significantly larger than the reference, 
then irrigation would have been given on May 24, that is 
10 days after the previous water application, some 22 and 
28 mm actual and potential transpiration, respectively, 
having elapsed. The same result could have been achieved 
both by similarly applying Jackson’s definition of CWSI, 
that is by assuming r,, = 0, and by using a more complex 
(radiation-dependent) model of r,, as shown by Stan- 
ghellini et al. (1992). One may conclude, therefore, that a 
mean daily pattern of reference CWSI, sufficiently ac- 
counts for weather variations, although it might be ar- 
gued that a shorter time interval could have been consid- 
ered as well. 
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Fig. 5. Measured « (that is, E/E,) as related to the soil water pres- 
sure head, h, for each day of the dry period, from May 17 (left) to 
May 28 (right). The full line is that of best-fit, whereas the horizon- 
tal lines represent the mean value of « for the whole reference period 
(dash-point) and its standard deviation (broken line) 


Soil-based water stress detection 


The function «(/) is shown in Fig. 5: the points represent 
estimates of « as the ratio of actual to potential transpira- 
tion, Eq. (10), corresponding soil water pressure heads 
being determined from 6, data (down to 30 cm) by means 
of Eq. (12). One could stipulate that extraction is signifi- 
cantly reduced, that is, «(A) is significantly smaller than 
unity, as soon as the ratio E/E, falls below its mean value 
minus o, that is 0.85, which is the value for |h| equal to 
40 kPa, according to the best-fit line shown in Fig. 5. 
Hence, irrigation would have been given by May 24 also 
by this method, although this would require an indepen- 
dent independent knowledge of the best-fit line. 

The literature would not be of much help in providing 
this information: De Jong and Kabat (1990), for instance, 
could quantify the “limiting point” (that is the pressure 
head at which extraction starts to decrease) no better 
than between 20 and 80 kPa for grass on sandy soil. Our 
best-fit line, moreover, would extrapolate to a relatively 
small absolute value (300 kPa) of the wilting pressure 
head, whereas published values for grass vary between 
400 and 800 kPa; Belmans et al. (1984) having contrib- 
uted the former and De Jong and Kabat (1990), the latter. 
In fact, this lack of precision is not surprising: in all 
experiments, including our own, a profile has to be con- 
trived from the sparse data available within the rooting 
depth, the latter being, for good measure, also coarsely 
estimated. Stanghellini et al. (1992) have shown that, in 
consequence of both these factors, the error on the hori- 
zontal position of data points such as in Fig. 5 becomes 
more relevant as / increases. Altogether, it is unlikely that 
a literature-derived relationship between « and h can be 
reliably applied to particular conditions for the early 
recognition of water shortage, despite the accuracy of the 
estimate of soil water content, demonstrated in this ex- 
periment. 


Conclusion 


The main conclusion to be drawn is that infrared temper- 
ature-based stress indices are more likely to be suitable 
for early detection of water shortage than soil based 
methods. Similarly, Itier et al. (1990), arguing that pre- 
dawn leaf water potential could be used as a good stress 
indicator in tomatoes, claimed that soil water content 
could not. By other means, Penuelas et al. (1992) working 
with greenhouse strawberries subject to two mildly stress- 
ing treatments, concluded that most variations on the 
CWSI theme were more reliable than soil water content- 
based indicators. 

On the other hand, it is known that protracted water 
shortage can induce morphological modifications restor- 
ing the balance between water extraction and transpira- 
tion by less damaging means than stomatal closure 
(Hsiao, 1993). Accordingly, sustained stress may reduce 
the effectiveness of infrared thermometry (De Lorenzi 
et al. 1993), whereas the impact of the inaccuracy of the 
soil-based stress indicator would decrease. 
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Abstract. Barley is a crop that has been classified as toler- 
ant to soil salinity, but under sprinkler irrigation with 
saline water it can readily absorb salts through its leaves 
and develop injury. Experiments using a triple-line- 
source sprinkler system were conducted on barley be- 
tween 1989 and 1991 to determine: (1) the specific effects 
of foliage wetting on the mass of different shoot compo- 
nents; (2) the relative contribution of root and foliar ab- 
sorption processes to foliar Cl accumulation; and (3) the 
extent by which these processes affect Cl partitioning in 
the shoot at the end of the season. Some plants were 
covered with plastic during the irrigation process to pre- 
vent foliar wetting while others remained uncovered. 
Salinity affected the partitioning of dry matter in the 
shoots regardless of whether plants were covered during 
the irrigation process. The organs associated with repro- 
duction, e.g., heads and peduncles, comprised a larger 
fraction of the total shoot biomass under high salinity 
than under low salinity, indicating that plants under 
salinity stress were able to redistribute their dry matter to 
favor reproductive growth. The Cl concentration of the 
young leaves sampled from uncovered plants was linearly 
related (i.e., r? > 0.71) to the Cl concentration of the 
irrigation water. Equivalent leaves from covered plants 
also contained a substantial amount of Cl but concentra- 
tions were weakly correlated (i.e., r? < 0.41) with the 
concentration of Cl in the irrigation water. At low salin- 
ity, there were no differences in leaf Cl concentrations 
between covered and uncovered treatments. In young 
leaves, differences between these treatments progressively 
increased with increasing salinity, indicating that the rel- 
ative contribution of Cl in the leaf from foliar absorbed 
salts increased with increasing Cl in the irrigation water. 
Only in the youngest leaves sampled at the end of the 
season from plants grown at high salinity was the Cl 
concentration in uncovered plants (foliar plus root-ab- 
sorbed Cl) found to be more than twice that in covered 
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plants (only root-absorbed Cl) indicating that most of the 
Cl in young leaves originated from foliar absorption. In 
addition, only in the youngest leaves (e.g., flag leaves) 
was the slope of the relationship between leaf-Cl concen- 
tration and Cl concentration of the sprinkling water of 
uncovered plants more than twice that of covered plants, 
also indicating that foliar-Cl absorption was more sub- 
stantial than root-Cl absorption. At high salinity, the 
difference in leaf Cl concentration between covered and 
uncovered plants was maximum in the youngest leaf (flag 
leaf), but differences became progressively smaller with 
increasing leaf age until ultimately concentrations of 
chloride in leaves older than the flag leaf-2 were highest 
in covered plants. In older tissue, it was difficult to distin- 
guish which process, foliar or root absorption, was most 
responsible for leaf-Cl accumulation. These processes 
may not be entirely independent of one another and 
much of the Cl in the oldest leaves of uncovered plants 
could have been derived from foliar sources during the 
first month of sprinkling, reaching maximal levels, and 
thereby restricting root-absorbed Cl. Furthermore, since 
these leaves at the end of the season are more injured and 
drier than those from covered plants, late-season sprin- 
kler irrigations may have been responsible for leaching 
some of the Cl out of these necrotic leaves. 





Sprinkler irrigation with saline water often decreases the 
yield of the crop more than if the same water is applied by 
surface or drip irrigation methods (Bernstein and Fran- 
cois 1973, 1975; Gornat et al. 1973; Meiri et al. 1982). 
This is not surprising since many crops under sprinkler 
irrigation absorb salts directly through their leaves and 
suffer foliar injury (Bernstein and Francois 1975; Busch 
and Turner 1967; Ehlig and Bernstein 1959; Francois and 
Clark 1979; Gornat et al. 1973; Harding et al. 1958; Maas 
et al. 1982a, b). Although foliar salt accumulation and 
injury has been considered to be a major cause of this 
yield reduction, no information is available to predict 
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yield losses as a function of salt concentration of the 
sprinkling water (Maas 1985). Foliar injury from saline 
sprinkling water, however, may not always translate into 
a yield reduction. Injury produced on many annual crops 
from six weeks of sprinkler irrigation in a greenhouse 
using water with an EC of 1.5 to 6.0 dS/m did not reduce 
yields below those of uninjured controls (Maas et al. 
1982b). In another case, plum trees sprinkled one season 
with saline water produced severe foliar injury but yields 
were not affected (Mantell et al. 1989). In the subsequent 
year, however, yields were reduced because of a decrease 
in fruit-set, even though sprinkling with saline water was 
discontinued. 

Plants that are sprinkle irrigated with saline water can 
accumulate salt in their leaves by two processes: (1) by 
root absorption and subsequent translocation to the 
leaves; and (2) by direct foliar absorption. Root absorp- 
tion is a rather selective process and most crop plants are 
relatively effective at excluding much of the salts in the 
soil solution (i.e., Cl and to a larger extent Na) from 
entering the leaves. The absorption of salts by leaves is 
less selective (Aragiiés et al. 1994) and Na and Cl in the 
sprinkling water can enter the leaf in proportion to their 
concentration (Maas et al. 1982a, b). Therefore the ad- 
vantage of ion-selective processes that the plant root has 
acquired through evolution may be lost under sprinkler 
irrigation (Grattan 1993). 

Susceptibility to foliar injury varies among crops but 
not necessarily in strict relation to the rate of foliar salt- 
absorption (Maas et al. 1982 b; Maas 1985). Tree crops, 
for example, are known to be susceptible to specific-ion 
toxicities when grown in saline soils (Maas 1990), but 
were found to vary in their sensitivity to sprinkler-applied 
salt solutions (Ehlig and Bernstein 1959). Leaves of de- 
ciduous tree crops (almond, apricot, and plum) and or- 
ange absorbed substantial amounts of Na and Cl and 
were severely damaged by saline sprinkling water where- 
as those from avocado, which absorbed Na and Cl more 
slowly, were not damaged. Most annual crops are not 
normally susceptible to specific-ion toxicities but can be 
particularly sensitive if sprinkle irrigated with saline wa- 
ter (Bernstein and Francois 1975; Gornat et al. 1973; 
Maas etal. 1982a,b). Tomato and potato readily ab- 
sorbed Na and Cl into their leaves under sprinkler irriga- 
tion with saline water and were severely injured whereas 
sorghum and cauliflower absorbed salts in their leaves 
more slowly and were only slightly injured. Safflower, on 
the other hand, readily absorbed salts by its leaves but 
was only slightly injured. 

Experiments at the US Salinity Laboratory showed 
that Na and Cl absorption by leaves from many herba- 
ceous crops, including barley, were not only linear func- 
tions of the hours sprinkled (Maas et al. 1982a, b) but 
that increased irrigation frequency contributed more to 
foliar salt accumulation than increased sprinkler dura- 
tion (Maas etal. 1982a). Although these greenhouse 
experiments were useful in distinguishing sensitivities 
among different crops to saline sprinkling water, they 
were not designed to determine the relative contribution 
of root- and foliar-absorption processes to foliar-Cl accu- 
mulation or to determine what impacts these absorption 


processes have individually, and collectively, on crop 
yield in field conditions. 

The experiments described here used a triple line- 
source sprinkler system developed and tested by Aragiiés 
et al. (1992) to: (1) evaluate the specific effects of foliage 
wetting on the biomass of different shoot components; 
(2) determine the relative contribution of root and foliar 
absorption processes to foliar Cl accumulation in field- 
grown barley; and (3) determine the extent to which these 
absorption processes affect Cl partitioning in the shoot at 
the end of the season. 


Materials and methods 


The triple-line-source sprinkler (TLS) system at the Servicio de In- 
vestigacion Agraria (SIA) field experimental site in Zaragoza, Spain, 
described in detail by Aragiiés et al. (1992), was used for our study. 
Three parallel sprinkler laterals were spaced 15 m apart, a distance 
equivalent to the wetted radius of a sprinkler. Sprinklers along the 
lateral were spaced 4.5 m apart. The center line discharged saline 
water while the outer lines discharged low-saline water. This ar- 
rangement produced a linear salinity gradient between the center 
and outer laterals while providing an even distribution of water 
(irrigation uniformity >90%). Minimum and maximum EC values 
and major-ion concentrations in the applied water for each season 
are presented in Table 1. 

Irrigations were scheduled according to evapotranspiration esti- 
mates over a 10-d period from pan evaporation data and site- 
specific crop-coefficients for non-stressed barley. Experimental plots 
were irrigated two to three times per week for 30 to 36 min per 
irrigation. The total number of irrigations throughout the experi- 
mental period in 1989, 1990, and 1991 were 34, 29, and 31 and the 
total quantity of water applied by sprinkling was 400, 417, and 
406 mm, respectively. In addition, plots received 151, 102 and 
127 mm of rain during the 1989, 1990, and 1991 season. All three 
sprinkler laterals were supplied with low salinity water for an addi- 
tional three-minute period at the end of each irrigation to remove 
saline water from the leaf surface and to reduce additional foliar salt 
absorption and injury. This post-irrigation, leaf-rinsing technique 
was also used by Frenkel et al. (1990) as a means of minimizing leaf 
burn. In 1991, a 3-min pre-irrigation with low saline water was also 
given as a means of reducing foliar absorption. 


Table 1. Minimum and maximum EC values (dS/m) and ion con- 
centrations (mmol/l) of the applied water 





Year 





1990 





Minimum values * 


EC : 2:5 1.6 
Cl ’ 10.5 6.6 
Na : 8.9 5.6 
Ca i 5 3.8 


Maximum values * 


EC 16.3 cg 14.7 
Cl 162 t/3 169 
Na 62 65 74 
Ca 57 60 52 





* Time and volume weighted average values that include both 
sprinkling and rain water. Minimum and maximum values reflect 
the composition of the applied water at the farthest and closest 
pluviometer from the center line 





Barley (Hordeum vulgare L.) seeds were sown in rows parallel to 
the sprinkler laterals on 22 Dec. 1988, 20 Jan. 1990 and 19 Nov. 
1990 and were harvested on 29 June 1989, 28 June 1990, and 1 July 
1991, respectively. In 1989, four Winter cultivars were tested; ‘Al- 
bacete’, ‘Barbarrosa’, ‘Criter’, and ‘Igri’. In 1990, four Spring culti- 
vars (‘“Georgie’, ‘Hassan’, ‘Patty’, and ‘Trait D’Union’) were tested 
since heavy rains in Nov. and Dec. 1989 delayed access to the field 
until Jan. 20; a date considered too late to re-plant the same Winter 
cultivars. Since no significant differences in leaf Cl accumulation 
were found among cultivars in either 1989 or 1990, only ‘Albacete’ 
was studied in 1991. The loam soil at the experimental site was 
classified as a mixed, mesic, Typic Torrifluvent. Barley stands were 
established by sprinkler irrigation with low-saline water until plants 
were at a suitable age to begin the saline sprinkling water treat- 
ments. Five locations, each 2.8 m apart, at increasing distances from 
each side of the central lateral were designated as plots to sample 
leaves and plants. Plots at each sampling location were 1.2 m in 
length and 1.4 m wide. Pluviometers were situated in each of the five 
plots on both sides of the central lateral to measure water quantity 
and to monitor EC. In 1989 and 1990, two of the four cultivars 
tested were planted on one side of the central (saline) lateral while 
the other two were planted on the other side. Before sprinkling with 
saline water began, some of the barley plants in each of the designat- 
ed plots were covered with transparent plastic during irrigations to 
avoid foliar wetting during the irrigation process but allowing the 
soil under the covered plants to become salinized. All other plants 
that were tested remained uncovered. Soil borders at the perimeters 
of each plot assured that the sprinkling water applied to the plot was 
contained within the plot. In 1990, however, plants in covered treat- 
ments were covered at all times with a transparent, corrugated 
plastic-roof. Although the plastic roof reduced incoming radiation 
(400-700 nm) by 20%, temperature and humidity was relatively 
unaffected. The plastic roof was sloped from 1.6 to 1.4 m above the 
soil surface. Sprinkling water runoff from the roof was collected in 
containers and then distributed to the respective plots via drip 
irrigation on the soil surface. Soil salinity measurements, taken at 0 
to 30 cm in each of the plots with a four-electrode probe, indicated 
that soil salinity was similar under covered and uncovered sections 
of the plots. In 1990, average Cl, concentrations from 0—30 cm for 
treatments one through five were 86, 98, 119, 140, and 166 mmol/l, 
respectively. These values were calculated from average soil EC, 
values (dS/m) using the expression; Cl (mmol/l) = —2.0+10.5 EC 
(dS/m). Sprinkling with saline water began on 6 Feb., 23 Feb., and 
31 Jan. in 1989, 1990, and 1991 and terminated on 2 June, 4 June, 
and 7 June, respectively. In 1989, leaves were sampled from plants 
on 6 April, 9 May and 23 May. In 1990, leaves were sampled on 
16 April, 7 May, and 29 May. When the sampling date was before 
heading, the youngest most fully expanded leaf was sampled. When 
the sampling date was after heading, the leaf below the flag leaf (flag 
leaf-1) was sampled. At each sampling, leaves were green and 
healthy except for those taken on 29 May, 1990. On this date, many 
of the leaves sampled from uncovered plots at moderate to high 
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levels of salinity were senescent or severely injured. In 1991, samples 
were taken on 21 May while most plants were still healthy. In this 
year, five samples, each containing two to three plants per sample, 
were taken in each of the five plots. The shoots in each sample were 
partitioned into seven separate components: (1) sheaths and culms; 
(2) heads; (3) peduncles; (4) flag leaf blades; (5) flag-leaf blades-1 
(i.e., the leaf blades below the flag leaf); (6) flag leaf blades-2 
(i.e., the leaf below flag leaf-1); and (7) old-leaf blades (i.e., all 
remaining leaves). No distinction was made between tillers and 
main stems before combining similar organs. At low salinity, only 
a small fraction of the “told leaves” were senescent at the time of 
sampling. At the other extreme, nearly all the “old leaves” from 
uncovered plants at high salinity were senescent as well as many of 
the “flag leaf-2” leaves. All other sampled leaves and tissues were 
green and healthy. The sampled plant material was washed two to 
three times with deionized water for 10 to 20s to remove salts on 
the tissue surface, and dried in an oven at 70°C to a constant 
dryness. The dry tissue was weighed and then ground in a blender. 
The chloride concentration in the plant tissue was determined on 
dilute nitric-acetic acid extracts of the ground plant material by the 
Cotlove (1963) coulometric-amperometric titration procedure. 


Results 
Shoot biomass 


In 1991, increased salinity of the irrigation water signifi- 
cantly reduced the biomass of all seven shoot compo- 
nents, but not in the same proportion. The fraction of 
whole-shoot biomass comprised by a particular compo- 
nent depended upon the salinity level. Depending upon 
the particular shoot component, this fraction either in- 
creased, decreased or was not affected by increasing 
salinity (Table 2). The values in the table are slopes, inter- 
cepts, and correlation coefficients (r) from linear regres- 
sion analysis of the fraction of shoot biomass relative to 
the irrigation water Cl concentration. The table shows 
that as salinity increased, the biomass of the head and 
peduncle became a larger fraction of the total shoot mass 
(i.e., positive correlation) while the biomass fraction of 
sheath plus culms and flag leaf became less (i.e., negative 
correlation). The mass fraction of the shoot components 
were affected by salinity the same, for the most part, 
regardless of whether they were covered or uncovered 
during sprinkler irrigation. 


Table 2. Slopes, y-intercepts and correlation coefficients (r) from linear regression analyses of the dry weights of the given organs expressed 


as fractions of shoot dry weight, on the Cl concentration (mmol/l) of the irrigation water for the covered and uncovered treatments. Data 
from the 1991 experiment 





Shoot organ Slope 





Covered Uncovered 


y-intercept 








Covered Uncovered Covered Uncovered 





Flag leaf 

Flag leaf-1 

Flag leaf-2 

Old leaves 

Head 

Peduncle 

Sheath and culm 


—6.74x 1075 
—7.30 x 1075 
—5.96 x 1075 
1.60 x 1075 
4.74x 1074 
1.89 x 1074 
—4.65x10~* 


—4.48 x 1075 
—1.69x10~> 
—5.45x 1075 
8.50 x 1075 
6.34 x 1074 
2.19 x 1074 
—8.34x 1074 


0.021 0.018 
0.043 0.039 —0.67 —0.23 
0.046 0.045 —0.75 —0.47 
0.095 0.099 0.09 0.43 
0.157 0.128 0.82 0.86 
0.104 0.095 0.67 0.72 
0.535 0.576 —0.73 —0.91 


—0.80 —0.53 











April 6.1989 


Y=773+ 2.43 X; r=0.89 


LEAF Cl (mmol/kg) 


Y=720+ 1.24 X; r=0.64 








May 9,1989 


Y =672+3.34 X; r=0.84 


LEAF Cl (mmol/kg) 


Y=684+0.71 X; r=0.38 














May 23,1989 


Y=835+2.77 X; r=0.87 
_8 =o oe 8 


LEAF Cl (mmol/kg) 


Y=786 + 1.81 X; r=0.57 











40 80 120 160 
IRRIGATION WATER Cl (mmol/L) 


Fig. la—c. Leaf Cl concentration of covered (0) and uncovered (m) 
barley (cvs ‘Albacete’, ‘Barbarrosa’, ‘Criter’ and ‘Igri’) sampled at 
various times in 1989 in relation to the chloride concentration of the 
irrigation water 


Tissue Cl concentration 


Leaf Cl concentration data indicate that both root and 
foliar absorption processes are responsible for leaf Cl 
accumulation. The extent to which each of these process- 
es contribute to the overall concentration of Cl in the leaf 
varies among organs and with sampling time. 

Linear regression analysis was conducted on the rela- 
tionship between leaf Cl concentration and the Cl con- 
centration in the applied water for each cultivar in both 
1989 and 1990. No significant differences (i.e., P>0.05) 
in slopes or y-intercepts were found among the different 
cultivars in either year. Consequently, data from the four 
cultivars in each of these years were combined. 

In 1989, the Cl concentration in leaves from uncov- 
ered plants was linearly (r?=0.71 to 0.79, P<0.001) re- 
lated to the concentration of Cl in the sprinkling water 
regardless of the sampling date (Fig. 1a—c). Leaf Cl con- 
centration in covered barley plants (i.e., leaves that were 
not wetted during sprinkler irrigation) was not correlated 
as well with the Cl concentration of the irrigation water 
as those from uncovered plants and r? values ranged 
from 0.14 to 0.41 depending upon the date sampled. 





April 16,1990 
Cov: Y =693 + 1.59 X; r=0.40 
Uncov: Y = 642+ 5.05 X; r=0.88 


LEAF Cl (mmol/kg) 








LEAF Cl (mmol/kg) 





Y=1119+3.02 X; r=0.61 





May 29,1990 


09.0-. 
aes 


LEAF Cl (mmol/kg) 


Cov: Y= 1555+ 13.0 X;r=0.87 g 
Uncov: Y = 2893-3.98 X; r=0.22 














40 80 120 160 
IRRIGATION WATER Cl (mmol/L) 


Fig. 2a—c. Leaf Cl concentration of covered (0) and uncovered (m) 
barley (cvs ‘Georgie’, ‘Hassan’, ‘Patty’ and ‘Trait D’Union’) sam- 
pled at various times in 1990 in relation to the chloride concentra- 
tion of the irrigation water 


As expected, differences in leaf Cl concentration be- 
tween covered and uncovered barley plants were unde- 
tectable at low Cl concentrations in the irrigation water, 
but increased as the Cl in the water increased. Differences 
between the slopes of the relationship between leaf Cl and 
Cl in the irrigation water, between covered and uncov- 
ered treatments, are significant (P<0.05) while there are 
no significant differences between the y-intercept values. 
This indicates that the relative contribution of Cl in the 
leaf from foliar absorbed salts increased with increasing 
Cl in the irrigation water. 

In 1990, leaves from both covered and uncovered 
plants were sampled at three different times throughout 
the season similarly to the previous season. The four 
cultivars tested in 1990 were ‘Georgie’, Hassan ‘Patty’, 
and ‘Trait D’Union’. Leaf Ci concentration data were 
again averaged and plotted (Fig. 2). 

The data from the first two sampling periods (16 April 
and 7 May) were comparable to those from the previous 
year in that: (1) leaf Cl concentration from uncovered 
plants increased linearly (r? =0.77 and 0.88) with increas- 
ing Cl concentration in the sprinkling water; and (2) there 
were no differences in leaf Cl between covered and uncov- 





ered plants at low salinity but differences increased with 
increased Cl in the sprinkling water (i.e., slopes were 
significantly different). One difference between the data 
from this year and those of the previous year was that at 
high salinity (i.e., >120 mmol Cl/l) at the second sam- 
pling, most of the Cl in the leaves originated from foliar 
absorption. That is, the concentrations of Cl in leaves of 
uncovered plants (root plus foliar-absorbed Cl) were 
more than twice those from covered treatments (only 
root-absorbed Cl). The difference may be attributed to 
the cultivars grown. In 1989, winter cultivars were tested, 
whereas in 1990 spring cultivars were grown. Since spring 
cultivars mature earlier, leaves sampled at the same date 
would have been physiologically older than those of the 
winter cultivars. 

Data from the third and latest sampling (29 May) in 
1990 was unlike those from the previous year. In 1990, 
leaf Cl concentrations in covered plants at high salinity 
(> 140 mmol Cl/l), were significantly higher than those 
from uncovered plants. Furthermore, the Cl concentra- 
tion in sampled leaves from uncovered plants at high 
salinity was substantially lower than the previous sam- 
pling. For each cultivar, the concentration of Cl in sam- 
pled uncovered leaves increased with increased Cl in the 
sprinkling water to a maximum, which varied among 
cultivars, then decreased with increased Cl in the sprin- 
kling water. Leaf Cl in covered plants, on the other hand, 
steadily increased as Cl in the sprinkling water increased. 


Distribution of Cl among the organs of the shoot 


In order to better understand the differences in results of 
the third sampling between the 1989 and 1990 years, 
covered and uncovered barley shoots were partitioned 
into various organs at the end of the 1991 season and 
analyzed for Cl. Tissue Cl concentrations in these organs 
are graphically expressed in Fig. 3. 

In 1991, leaf Cl concentration in uncovered flag and 
flag-1 leaves was linearly related (r? 0.78 and 0.86 respec- 
tively) to the concentration of chloride in the sprinkling 
water. Differences in leaf Cl between covered and uncov- 
ered plants progressively increased with concentration of 
Cl of the sprinkling water. This behavior was similar to 
those shown for the first two samplings in 1989 and 1990. 
This similarity is not surprising since all these leaves were 
relatively young. Thus, it is inferred that most of the Cl 
that accumulated in the flag leaves of uncovered plants 
from locations where Cl in the irrigation water exceeded 
100 mmol/I was by foliar absorption. This assumes, how- 
ever, that root and foliar absorption processes operate 
independently. 

The 1991 tissue Cl data shown in Fig. 3 indicate that 
Cl was partitioned differently in covered than in un- 
covered barley shoots. Young tissues (e.g., the peduncles 
and head) contained the lowest concentration of Cl re- 
gardless of whether they were exposed to saline spray. 
This result was expected since these new organs had rela- 
tively little time to accumulate Cl as compared to older 
tissue. However, at the highest salinity level, heads and 
peduncles from uncovered plants contained significantly 
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more Cl than those that were covered during sprinkling 
irrigation. The sheaths and culms, which contributed to 
approximately half the shoot biomass, contained a sub- 
stantial amount of Cl and were responsible to a large 
extent for sequestering root-absorbed Cl from actively 
photosynthesizing leaf blades (Wolf et al. 1991). 

Leaf Cl concentrations are of particular interest. At 
low salinity levels, leaf Cl of covered and uncovered 
plants was not significantly different but leaf Cl concen- 
trations steadily increased from young to older leaves 
where concentrations reached a maximum within the flag 
leaf-2. At high salinity, the difference in leaf Cl concen- 
tration between covered and uncovered plants was maxi- 
mum in the youngest leaf (flag leaf), but the difference 
became progressively smaller with increasing leaf age and 
concentrations in old leaves were highest in covered 
plants. In covered plants, leaf Cl increased with increas- 
ing Cl in the irrigation water, especially with increased 
leaf age. The slope of this relationship progressively in- 
creased (i.e., 0.87, 2.99, 4.03 and 5.71) as leaf age in- 
creased (i.e., flag, flag-1, flag-2, and old leaves). In con- 
trast, leaf Cl in uncovered plants sampled at the end of 
the season, did not increase with increasing leaf age at the 
higher salinity levels. Instead, the slope of the relation- 
ship between leaf Cl and irrigation water Cl concentra- 
tions progressively decreased (i.e., 8.39, 6.71, 4.19, and 
3.02) with increasing leaf age. 

The effects of the interaction between salinity level 
and leaf age on leaf Cl concentration (mmol/kg) can be 
seen in Fig. 4. In this figure, the ratio of leaf Cl concentra- 
tion in uncovered to that in covered plants at low, medi- 
um and high irrigation water Cl concentration is plotted 
in relation to the leaf position (leaf age) on the shoot. In 
all cases, the ratio decreases as leaf age increased. Howev- 
er, the decrease in the ratio was greater the higher the 
salinity level. As salinity increased, the ratio increased 
in younger leaves. In the flag leaf, the ratio was greater 
than two at high salinity indicating that most of the Cl in 
the leaf was from foliar absorbed Cl. However, in older 
leaves (i.e., flag leaf-2 and older), moderate to high 
salinity (72 to 188 mmol Cl/l) decreased the ratio of chlo- 
ride concentration of uncovered leaves to that in covered 
leaves. In the oldest leaves, the ratio was less than unity 
at all salinity levels, emphasizing that the concentration 
of Cl in covered leaves was greater than that of leaves 
uncovered during sprinkler irrigation. 


Discussion 


Yield data, reported by Aragiiés et al. (1994), indicate 
that saline sprinkling water reduced shoot biomass of 
barley more when plants were subjected to foliar wetting 
than when they were covered during irrigation. Further- 
more, they demonstrated that leaf Cl concentration was 
linearly related (i.e., negative slope) to yield. 

In these experiments, barley plants that were uncov- 
ered during irrigation exhibited substantially more foliar 
injury (i.e., chlorosis and necrosis), particularly in older 
leaves, than those covered during irrigation. This is not 
surprising since older leaves have a longer opportunity 
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Fig. 3a—g. Tissue Cl concentration of various shoot components of 
covered (0) and uncovered (m) barley sampled at the end of the 1991 
season, in relation to the Cl concentration of the sprinkling water. 
Flag leaf-1 indicates the leaf just below the flag leaf, flag leaf-2 the 
next leaf and old leaves all remaining leaves. Individual points are 
mean values with +S.E. (P=0.05) 
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FLAG FLAG -1 FLAG -2. OLD LEAVES 


LEAF POSITION ON THE SHOOT 


Fig. 4. Leaf Cl concentration ratio of those plants uncovered (foliar 
plus root-absorbed Cl) to those covered (root-absorbed Cl only) 
during sprinkler irrigation. Leaf Cl ratio’s are for different aged 
leaves on the plant each at low, medium, and high levels of salinity. 
Data from the 21 May, 1991 sampling 


to accumulate salt by foliar and/or root-absorption pro- 
cesses (depending upon the treatment) than do younger 
leaves. 

Salinity affected the partitioning of dry matter in the 
shoot whether plants were covered during the irrigation 
process or not. In 1991, the organs associated with repro- 
duction (e.g., heads and peduncles) comprised a larger 
fraction of the total shoot biomass under high salinity 
conditions than under low salinity. The biomass of other 
tissue, such as young leaves and sheaths plus culms, 
was a correspondingly smaller fraction of the total shoot 
weight under high salinity than under low salinity. This 
suggests that under salinity stress the plants were able 
to redistribute their dry matter to favor reproductive 
growth. 

Leaf Cl concentration of young leaves sampled from 
plants that were uncovered during irrigations was in 
nearly all cases, linearly related to the concentration of Cl 
in the irrigation water. Equivalent leaves from covered 
plants also contained a substantial amount of Cl but leaf 
Cl concentrations were not as well correlated with Cl 
concentration in the irrigation water as salinity increased. 

The extent by which foliar absorption or root absorp- 
tion processes were contributing to overall leaf Cl con- 
centration was dependent upon the shoot tissue and sam- 
pling date. The contribution of each process can be eval- 
uated by two approaches. To assess the fraction of the 
leaf Cl that was derived from the soil via root absorption 
in comparison to that derived from the leaves via foliar 
absorption, it is appropriate to compare leaf Cl concen- 
trations between covered and uncovered leaves at a given 
salinity level. If leaf Cl from plants that were uncovered 
(i.e., root plus foliar absorption) is more than twice that 
from the covered (i.e., root absorption only) treatment, 
then more than half of the Cl may be assumed to be 
derived from foliar absorption. On the other hand, a 
comparison of the slopes of the relationship between irri- 
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gation water and tissue Cl concentration in covered and 
uncovered treatments may indicate the relative impor- 
tance of foliar and root absorption processes. Thus, if the 
slope of the uncovered treatment is more than twice that 
from the covered treatment, then the foliar absorption 
process dominates over the root absorption process in 
terms of foliar Cl accumulation. These comparisons are 
valid when y-intercepts are similar and assume that foliar 
absorption and root absorption processes are, for the 
most part, independent of one another in terms of leaf Cl 
accumulation. However, over the duration of the experi- 
ment, it is obvious that foliar absorption and root ab- 
sorption will not always be entirely independent of one 
another. For example, due to differences in foliar injury, 
transpiration rates in uncovered plants towards the end 
of the season will likely be less than those that were cov- 
ered during irrigations. Similarly, the cumulative transpi- 
ration of plants that were covered throughout the 1990 
season may have been less than those that were uncov- 
ered. Since transpiration affects the transport of root-ab- 
sorbed Cl within the plant, large differences in cumulative 
transpiration could affect root Cl accumulation. 

We found that in the uncovered treatment, these 
slopes were more than twice those from the covered treat- 
ment indicating that foliar absorption was the dominant 
process in certain young leaves (i.e., those sampled 9 May 
1989; 16 April and 7 May 1990; and 21 May 1991). Our 
explanation of this is as follows: Root-absorbed Cl accu- 
mulates in stems and older leaves of barley shoots while 
the young organs are protected against uptake and accu- 
mulation (Wolf et al. 1991). Therefore the concentration 
of Clin young tissue from root absorption and transloca- 
tion processes, is low relative to older tissue. This ion 
discrimination does not occur under foliar absorption 
processes and young tissue exposed to saline spray can 
readily absorb Cl in its tissue. It is hypothesized that salts 
enter the foliage through the cuticle rather than the stom- 
ata (Franke 1967; Haynes and Goh 1977). Differences 
in cuticular composition or thickness may exist between 
young and old leaves, thus resulting in a greater or lesser 
degree of salt absorption through the foliage. The tissue 
Cl concentration in peduncles and heads (both relatively 
young tissue) also increased linearly (r?=0.78 and 0.84, 
respectively) with increased Cl in the irrigation water. 
The slopes of uncovered treatments were both more than 
twice their respective covered treatments indicating that 
direct foliar absorption by these organs was the dominant 
process responsible for Cl accumulation. 

Some grasses appear to preferentially partition Cl to 
their sheaths relative to their blades (Greenway 1962) 
though this effect was found to be more pronounced in 
sorghum than barley (Boursier et al. 1987). In our study, 
sheaths were not separated from culms, yet their com- 
bined concentration of Cl along with their relatively high 
percentage of biomass, suggest that much of the total 
shoot Cl is present in these tissues. 

There were two samplings where Cl in the leaves from 
uncovered treatments was more than twice that of the 
respective covered treatments. At the 7 May, 1990 sam- 
pling, chloride in the leaves from uncovered plots was 
more than twice that of covered plots at irrigation water 





154 


Cl levels in excess of 117 mmol/l. Similarly, flag leaves 
from uncovered plots sampled on 21 May, 1991 con- 
tained over twice the leaf Cl of those from covered treat- 
ments at irrigation water Cl levels in excess of 105 mmol/l 
(Fig. 4). In both these instances where leaf Cl in uncov- 
ered leaves was more than twice those in covered leaves, 
it is reasonable to conclude that a larger fraction of the Cl 
in the uncovered leaves originated from foliar absorption 
and that this process is mostly responsible for Cl accumu- 
lation. However, in all other instances where leaf Cl in 
uncovered leaves is less than twice that in covered leaves, 
one can not safely conclude that most of the Cl originated 
from root absorption processes. Since Cl influx into root 
tissue was found to be dependent, to some extent, upon 
internal Cl concentrations (Cram 1983), it may be possi- 
ble that foliarly-absorbed Cl affects the accumulation of 
root-absorbed Cl in the leaf. Furthermore, leaf Cl accu- 
mulation in barley was found to be linearly related to the 
number of hours sprinkled in greenhouse studies for up 
to five weeks (Maas et al. 1982b). Additional sprinkling 
did not increase leaf Cl further. It is possible, therefore, 
that much of the Cl in older leaves of uncovered plants 
at higher salinity was derived from foliar-absorption 
sources during the first month of sprinkling, reaching 
maximal levels, and thereby restricting root-absorbed Cl 
that would otherwise be accumulating in the leaf if the 
plant was not wetted by saline water. Therefore, the as- 
sumption that foliar and root Cl absorption processes are 
independent of one another may be valid during early 
stages but not necessarily later on. This hypothesis war- 
rants further testing. 

Exposure of barley plants to saline spray caused leaves 
to prematurely senesce. Older leaves from uncovered 
plants at moderate to high salinity levels were severely 
injured. They were often much drier at the final sampling 
than those from covered plants and contained less Cl 
(Fig. 3, old leaves). This response was also observed on 
younger tissue at the final sampling in 1990 (Fig. 2, 
29 May) and in an experiment conducted at this station 
in 1987 (data not shown) when uncovered leaves were 
substantially drier than covered leaves. Maas et al., 
(1982a) also noted that for the initial five weeks, foliar 
Na and Cl concentrations increased in the leaves of pep- 
per plants that were sprinkle irrigated with saline water, 
but salt concentrations often decreased thereafter when 
leaves became injured and senescent. Although they sug- 
gested that Cl was more likely to be the cause of injury 
than Na, there was no consistent correlation between 
foliar injury and leaf Cl concentration. Although unlike- 
ly, during the desiccation process some of the Cl in the 
injured and often dry leaves from these studies may be 
translocated out of the leaves to other tissue. More likely, 
the leaf Cl concentrations are lower in older, often dry 
leaves than those from covered plants due to leaching of 
leaf Cl during late season sprinklings. In a very limited 
test in the laboratory, leaching of Cl from injured leaves 
during two, 30-s rinse periods with deionized water to 
remove leaf-surface salts, on the other hand, was found 
to be small (Grattan et al. 1981). Further examination of 
the leaching of ions from dry and senescent leaves is 
warranted. 


These experiments demonstrated that barley shoots 
wetted by saline sprinkling water, readily accumulated Cl 
in their young leaves by direct foliar absorption. Further- 
more, foliar absorption was more important for leaf Cl 
accumulation than was root absorption. However, leaf 
Cl concentration data from older leaves of plants that 
were covered and uncovered can not be used to determine 
the fraction of leaf-Cl that accumulated from root ab- 
sorption versus that from foliar absorption. Controlled 
solution culture studies with labeled Cl in the growth 
media and non-labeled Cl in the sprinkling water is need- 
ed to advance our understanding of the interactions be- 
tween these two absorption processes and determine the 
source of the leaf Cl (i.e., foliar versus root absorption) 
among different aged leaves. 
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Abstract. The main objective of the study was to evaluate 
the sensitivity of Iraqi local barley cultivar (Black) to soil 
water deficit as compared to other barley cultivars, name- 
ly: CM-72 and Arivat. The local cultivar proved to be 
susceptible during germination and emergence, but it de- 
veloped resistance during the vegetative growth and yield 
formation stages. Growth analysis of individual leaves 
clearly showed that Black barley produced leaves of long 
growth duration which could affect the growth processes 
especially cell division. Thus, the reduction in the number 
of cells composing those leaves was small in the cultivar 
Black. Leaf growth rates as soil matric potential de- 
creased, were reduced considerably, although, no clear 
trends were observed between cultivars. On the other 
hand, the smaller cell volume of local cultivar may have 
a role in its ability to tolerate water deficit. Some evidence 
that the local cultivar is post-flowering resistant is dis- 
cussed. 





The need to develop crops with higher tolerance to envi- 
ronmental stresses has increased tremendously in the last 
twenty years. A genetic approach has been adopted to 
develop some crop species that can tolerate environmen- 
tal stresses (Epstein 1983; Krizek 1984). There is consid- 
erable evidence that genetic variability in sensitivity to 
water stress exists within crop species such as barley, 
wheat, sorghum, and pearl millet (Omara 1987; Bidinger 
et al. 1987; Winter et al. 1988). This can be considered as 
experimental materials for comparative studies in stress 
physiology, and may lead to discovery of structural, 
physiological and biochemical features consistently asso- 
ciated with the ability to tolerate drought (Morgan et al. 
1986; Wenzel 1987; Morgan 1988; Rascio et al. 1990). 
Moreover, such studies may provide plant breeders with 
information to screen genotypes that can withstand water 
stress. However, hypothesis concerning specific plant 
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characteristics which may be associated with drought re- 
sistance can be developed mainly through knowledge of: 
(1) growth variables (rate and duration) and the contri- 
bution of growth processes (cell division and expansion) 
to plant growth (Yasseen 1983; Yasseen et al. 1987), and 
(2) performance of crop plants under stress conditions 
(Winter et al. 1988; Al-Maamari 1989). This study was 
conducted to investigate the importance of growth vari- 
ables and the contribution of growth processes to plant 
growth and yield, and the possibility of using these 
parameters in screening methods for drought resistance 
in barley. Also, the study included a comparative evalua- 
tion of resistance to water stress of Black barley which is 
considered to be tolerant to water stress and salinity (Yas- 
seen et al. 1987; Al-Maamari 1989). 


Materials and methods 


Plant material and soil 


Three barley cultivars were chosen for the study, California Mari- 
out-72 (CM-72), and Arivat, (Hordeum vulgare L.), and an Iraqi 
local cultivar, Black, (Hordeum distichum L.). Some agronomic, 
genetic and physiological characteristics of these cultivars are given 
in Table 1. Seeds of CM-72 were obtained from R. C. Sherwood 
Grain Co. Los Banos, California 93635 U.S. A. in November 1989, 
while seeds of Arivat and Black were obtained from Testing and 
Certification Branch, The Ministry of Irrigation and Agriculture, 
Mosul, Iraq, in October 1989. The viability test for all these seeds 
showed that 97 to 99% of the seeds used were viable. 

The soil, a silty clay loam, used in the study was obtained from 
Rabeaa area (about 120 km north of Mosul) from the Northern 
Irrigation Jazera Project. The soil chemical and physical properties 
are listed in Table 2. 


Treatments and rewatering method 


Four moisture levels were determined from a standard curve relat- 
ing soil matric potential to soil moisture content. The relationship 
was derived according to the method described by Richards (1954), 
using the pressure cooker technique to determine high matric poten- 
tials and the pressure membrane technique to determine low matric 
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Table 1. Some agronomic, genetic and physiological characteristics of three barley cultivars 





Cultivar Tiller number Yield 


per plant 


Origin Length of growing 


season (days) 


(kg ha~*) 


Number of 
grains per spike 


Salt 
tolerance 


Protein content 
in the grain (%) 


Spike 
type 





CM-72 
Arivat 
Black 


Short * 
Short (155-160) 
Long (180-210) 


Few 
Few 
Many 


High * 


High (3000) 
Low (1200) 


Six-row * 
Six-row 60 
Two-row 26 


Tolerant 
Sensitive 
Tolerant 





* No documented data in Iraq 


Table 2. Chemical and physical properties of the soil 





Properties Values 





pH (water extract) 7.8 
Sat. Ext. ECe (mS cm~') 0.62 
Sand (g kg~‘ dry soil) 27 
Silt (g kg~! dry soil) 

Clay (g kg~* dry soil) 

Exch. Na (meq kg~? dry soil) 

Exch. K (meq kg~! dry soil) 

Ca (meq dm~') 

Mg (meq dm~*) 

S (meq dm~'*) 





Table 3. Effect of soil matric potential on seedling emergence, as a 
percentage of control, of three barley cultivars after 14 days from 
sowing 





Soil matric potential (MPa) Black Arivat CM-72 





Field capacity (—0.03) 100 100 
—0.2 96 92 
—0.6 101 104 
—1.0 100 67 





Cultivars: 
Treatments: 
Interaction: 


P< 0.001 
P <0.001 
P <0.001 


potentials. Each level was prepared by adding appropriate amount 
of water to the soil, mixed thoroughly, to ascertain complete homo- 
genization as well as avoiding aggregation of soil particles. These 
levels were —0.03 (field capacity), —0.2, —0.6 and —1.0 MPa. 
Seeds were sown in plastic pots, 15 x 15 x 20 cm in size, each hold- 
ing 4kg dry soil. All treatments were replicated three times, and 
distributed randomly inside a glasshouse. Seedlings were thinned to 
four per pot after emergence had stopped in each treatment; 14 days 
after sowing in —0.03, —0.2, and —0.6 MPa; and 21 days at 
—1.0 MPa. The experiments_were conducted during the period 
January to April 1990 under natural conditions of light with photo- 
synthetically active radiation (PAR) of 26 to 35 Wm ~? and with 
average relative humidity ranging from 52 to 68%. The glasshouse 
was heated during the winter by electric heaters to maintain the 
temperature in the range of 18 to 23°C. 

Maintaining the water content of soil, and consequently the 
specified matric potential during germination, was accomplished by 
covering the pots with pieces of nylon to prevent evaporation of 
water from the soil surface until seedling emergence had stopped. 
During this period there was no need to rewater plants since 
changes in the weight of pots were not observed. However, during 
subsequent stages of growth and development the water content of 
the soil was maintained to the specified levels by adding distilled 
water every two days to the pots in holes of 10 cm in depth and 2 cm 
in diameter, in the center of the pots to hold them at their original 


weights, and to ensure water diffusion equally throughout the soil. 
Moreover, additional amount of water was added to the soil, in the 
same way, representing what the fresh weight of the growing plants 
can add to the original weight of pots. Estimation of fresh weights 
of growing plants was done weekly on extra plants grown under the 
same treatments. 


Measurements 


Seedling emergence. Twenty five seeds were sown in each pot, and 
the emergence of seedlings was followed daily and determined after 
fourteen days under the different levels of soil matric potential. 


Leaf growth. Leaf length was determined daily, by a ruler, from the 
beginning of emergence of the leaf tip from the sheath until appear- 
ance and exposure of its ligule. This was done for the third and fifth 
leaves. The growth of these leaves was analyzed according to the 
equation and method described by Gallagher (1979). The equation 
used was: F = (R x D) + 0.1 F, where F is final length (cm), R is 
growth rate (cm day~') and D is duration of growth (days). This 
equation has been used to describe leaf growth in several barley 
cultivars (Yasseen et al. 1987). 

The area of individual leaves (third, fifth and flag), and the 
remaining leaves on the main stem were determined according to 
the method described by Watson (1937), Aase (1978) and Yasseen 
et al. (1988). 


Cell number and volume. The number and volume of cells were 
determined using the technique described by Sunderland (1960), 
and, modified by Yasseen (1983). Cells were counted using a double 
chambered Hawksley haemocytometer slide. The cell volume was 
estimated by the following formula, assuming that the cell density 
is equal to the density of water: 


Fresh mass/Cell density 
Cell number 





Cell volume = 


Yield. Plants were harvested at the end of the season. Total number 
and weight of seeds, and the weight of 1000 seeds were determined. 


Results 
Growth and yield 


Seedling emergence, leaf area and yield were reduced as 
the soil matric potential decreased, but the cultivars dif- 
fered significantly in their response to water shortage at 
different stages of growth. For example, Arivat and 
CM-72 were more tolerant than the local cultivar during 
germination and seedling stages (Table 3). The reduction 
in the percentage of seedling emergence of Black was 
substantial especially at the lower levels of soil matric 
potential. During subsequent growth stages, however, 
the local cultivar seemed to gradually develop resistance 
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Fig. 1a, b. Effect of soil matric potential on the leaf area of the main stem, as a percentage of control, 50 a and 100 b days after sowing 
in three barley cultivars. The figures on the top of bars represent the standard errors of the means 


Table 4. Effect of soil matric potential on the number and weight of seeds of three barley cultivars 





Soil matric potential 
(MPa) 


Number of seeds/plant 
Black Arivat CM-72 


Weight of seeds (g plant™ ') 
Black 


1000-seeds weight (g) 


Arivat CM-72 Black Arivat CM-72 





Field capacity (—0.03) 
—0.02 

—0.6 

—1.0 


45.3 
38.3 
18.0 
14.0 


33.6 
22.6 

9.6 
10.6 


30.3 
19.3 
10.3 

6.0 


1.70 
1.47 
0.66 
0.49 


1.47 
1.03 
0.30 
0.32 


1.52 
0.85 
0.33 
0.19 


43.6 
45.6 
31.4 
31.0 


50.3 
45.3 
BS 
32.9 





Cultivars 
Treatments 
Interaction 


:P < 0.001 
:P < 0.001 
:N.S. 


to water stress. The reduction with water stress in the 
total leaf area on the main stem was considerable in Ari- 
vat and CM-72 compared to the local cultivar (Fig. 1). 
This large difference between the three cultivars became 
more pronounced at the end of the season. The data of 
yield components showed almost the same trend found in 
the vegetative growth (Table 4). The cultivar Black had 
more seeds and these seeds were heavier than the other 
two cultivars at lower soil matric potentials. Also, it is 
interesting to note that the 1000-seeds weight showed 
that the cultivar CM-72 produced the biggest seeds with 
no stress followed by Arivat, while the cultivar Black had 
the lowest seed weight. However, under stress conditions, 
the reduction in the 1000-seeds weight was substantial in 
CM-72 followed by Arivat, while the reduction in Black 
cultivar was the least. 


:P < 0.001 
:P < 0.001 
:N.S. 


:N:S. 
:P < 0.001 
:P<0.01 


Growth processes 


The reduction in leaf area due to water stress can be 
attributed mainly to the reduction in cell number and cell 
volume. The reduction in both growth processes was 
highly significant especially at the lower soil matric po- 
tentials. There were obvious differences among the three 
cultivars in the area of individual leaves as well as the 
number of cells in those leaves (Tables 5 and 6). At low 
soil matric potentials the reduction in the number of cells 
was substantial in the cultivars CM-72 and Arivat com- 
pared to Black, whereas no clear differences were found 
between these cultivars in the volume of cells (Table 7), 
with one exceptional case that the cell volume in flag leaf 
of CM-72 at — 1.0 MPa was reduced considerably (about 
42%). 
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Table 5. Effect of soil matric potential on the final area (cm?) of leaf 3 and 5 and flag leaf of three barley cultivars 





Soil matric potential 
(MPa) 


Leaf 3 


Black Arivat 


CM-72 


Leaf 5 
Black Arivat 


CM-72 


Flag Leaf 
Black 


Arivat 





Field capacity (—0.03) 
—0.2 
—0.6 
—1.0 


9.8 
10.7 
8.0 
7.6 


12.0 
10.7 
9.4 
6.2 


12.6 
10.2 
of 
3.0 


9:7 
10.4 
7.4 
1.2 


12.0 
12 
7 | 
2.4 


Re 
=| 
1.6 
0.9 


17 
12 
0.7 
0.6 





Cultivars 
Treatments 
Interaction 


NCS. 
:P < 0.001 
NGS. 


:P < 0.001 
:P < 0.001 
:P<0.01 


:P < 0.001 
:P < 0.001 
:P <0.01 


Table 6. Effect of soil matric potential on cell number (10°) in leaf 3 and 5 and flag leaf of three barley cultivars 





Soil matric potential 


Leaf 3 


(MPa) 


Black 


Arivat 


CM-72 


Leaf 5 


Black Arivat 


CM-72 


Flag Leaf 
Black 


Arivat 





Field capacity (—0.03) 


—0.2 
—0.6 
—1.0 


3.914 
4.100 
2.907 
2.840 


3.643 
3.872 
3.291 
2.220 


2.752 
1.900 
1.594 
0.826 


3399 
3.637 
2.359 
2.514 


3.456 
3.765 
0.827 
0.876 


2.592 
1.910 
0.399 
0.209 


0.632 
0.448 
0.406 


0.700 
0.529 
0.262 


0.435 


0.180 





Cultivars 
Treatments 
Interaction 


:P < 0.001 
[P<O0l 
N.S. 


:P < 0.001 
:P < 0.001 
:N.S. 


:P < 0.001 
:P < 0.001 
:N.S. 


Table 7. Effect of soil matric potential on cell volume (mm?) in leaf 3 and 5 and flag leaf of three barley cultivars 





Soil matric potential 
(MPa) 


Leaf 3 


Black Arivat CM-72 


Leaf 5 Flag Leaf 


Arivat CM-72 Black Arivat CM-72 





Field capacity (—0.03) 
—0.2 
—0.6 
—1.0 


96 
92 
63 
57 


67 
52 
52 
58 


61 
52 
56 
52 


53 
43 
47 
59 





Cultivars 
Treatments 
Interaction 


:P < 0.001 
:P < 0.001 
-NCS. 


Table 8. Effect of soil matric potential on the growth rate 
(cm day~') of leaf 3 and 5 of three barley cultivars 


:N.S. 
:P<0.01 
:P < 0.05 


Table 9. Effect of soil matric potential on the duration of growth 
(days) of leaf 3 and 5 of three barley cultivars 





Soil matric 
potential 
(MPa) 


Leaf 3 
Black Arivat CM-72 


Leaf 5 
Black Arivat CM-72 





Soil matric 
potential 
(MPa) 


Leaf 3 
Black Arivat CM-72 


Leaf 5 
Black Arivat CM-72 





Field capa- 2330. 196: “24 
city (—0.03) 
—0.2 
—0.6 


—1.0 


2435 29D. TT 


1.85 
1.74 
1.79 


Be 
2:23 
1.88 


Field capa- 8.0 
city (—0.03) 
—0.2 
—0.6 


—1.0 


9.6 19 10.3 8.4 6.3 


10.4 
9.4 
9.0 


9.9 
8.5 
8.5 


9.1 
43 
1D 


11.0 
99 
11.5 


pa 
5.6 
4.7 


Safi 
a 
4.3 





Cultivars 
Treatments 
Interaction 


:P< 0.05 
:P < 0.001 
:P <0.01 


A reduction in growth processes can be explained by 
a reduction either in the rate and/or duration of those 
processes. In the present study the main reduction in the 
length and area of leaves due to soil water deficit was 
attributable to the reduction in both rate and duration in 





Cultivars 
Treatments 
Interaction 


:P <0.01 
:P<0.01 
:N.S. 


:P < 0.001 
:P <0.01 
:P<0.01 


CM-72 and Arivat, (Tables 5, 8, and 9), whereas the 
growth duration of leaves in Black did not change with 
water deficit. A reduction in the growth rate of leaves was 
the main reason for the reduction in the length and area 
of the local cultivar with water deficits. 





Discussion 


The present findings indicate that cultivars differ in their 
reaction to soil water deficit during various stages of their 
life cycle. CM-72 is tolerant during germination and 
seedling emergence, while the local cultivar proved to be 
susceptible during those stages. The latter cultivar 
seemed to gradually build up a degree of tolerance to 
water deficit with age, so it had a larger area of leaves 
under all levels of soil matric potential compared to the 
other cultivars. Moreover, its superiority over the other 
cultivars may become more obvious during the post-an- 
thesis period and at maturation (Davidson and Birch 
1978; Al-Maamari 1989). The local cultivar was sensitive 
to salinity during germination and early seedling stage, 
but tolerant during the vegetative growth period (Yasseen 
et al. 1987 and 1988). Changes in the colloidal-chemical 
properties of the cytoplasm with respect to osmotic stress 
could be a reason for changes in sensitivity to water 
deficit during different stages of growth and development 
(Levitt 1980). CM-72 has been considered to be salt 
tolerant during germination and seedling stage (Ayers 
et al. 1952; Norlyn 1980), and may be tolerant to water 
deficit conditions. However, CM-72 proved to be less 
tolerant than the local cultivar during the vegetative 
growth period in the present experiments. Therefore, 
there is no particular stage in the crop life cycle that can 
be exploited for reliable screening tests, since the sensitiv- 
ity to osmotic stress changed dramatically during the dif- 
ferent stages of crop growth and development. Unfortu- 
nately much of successful work with breeding pro- 
grammes and screening techniques done during the last 
decade were based on a particular stage of the crop life 
cycle (Epstein et al. 1980; Blum et al. 1980; Epstein 1983; 
Omara 1987). Moreover, large scale screening techniques 
for breeding drought resistance in winter wheat have 
been used (Winter et al. 1988) with limited success for the 
determination of physiological characteristics that can be 
considered as important criteria in selecting wheat culti- 
vars suitable for cultivation in the semi-arid areas, and 
for breeding programmes in the development of cultivars 
resistant to drought. Therefore, much work remains to be 
done to develop and improve screening techniques for 
drought resistance, since, breeding for such conditions 
has been limited by a lack of simple screening methods, 
and a lack of evidence that the methods employed would 
result in improved yield stability under drought. One pos- 
sible method which seems worthwhile is by increasing the 
length of the growing season, and if leaves are able to 
expand rapidly under drought and salinity stresses, may 
lead to considerable maintenance of yield under these 
conditions (Monteith 1977). It has been shown that long- 
season cultivars of wheat are generally more tolerant to 
salinity during the vegetative growth and yield formation 
periods than short-season cultivars (Torres-B and Bing- 
ham 1973; Torres-B et al. 1974). Leaf growth in several 
cultivars under salinity stress was described quantitative- 
ly using logistic equations and the derived variables were 
calculated from such description (Yasseen 1983). Faster 
rates and longer durations of growth of individual leaves 
were found in long-season cultivars than in short-season 
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cultivars. This was considered an important criterion in 
selecting wheat cultivar suitable for cultivation in saline 
soils, and for breeding programmes to develop salt toler- 
ant cultivars. Data of this study confirm the importance 
of growth variables (rate and duration) in the evaluation 
of barley cultivars in their tolerance to water deficit. It is 
likely that longer duration of growth processes in leaves 
of the local cultivar led to a larger leaf area, and to more 
and larger seeds under conditions of water deficit com- 
pared to the other cultivars. Moreover, these leaves had 
more cells, indicating that the process of cell division in 
the local cultivar is relatively less sensitive to water deficit 
(Uchimiya and Takahashi 1973; Yasseen et al. 1987; Al- 
Maamari 1989). However, cell volume in these cultivars 
was reduced in a similar manner with stress, but the small 
volume of cells in leaves the local cultivar may have a role 
in resisting water stress (Cutler et al. 1977). 

Finally, it is interesting to note that Arivat was almost 
as sensitive as CM-72 to water deficit. This unusual ob- 
servation needs further investigation considering its ad- 
aptation to the environmental conditions in Iraq. 


Acknowledgement. The authors wish to offer special thanks to pro- 
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Book reviews 


Alscher RG, Cumming JR (eds): Stress responses in plants: 
adaptation and acclimation mechanisms. Wiley-Liss, New 
York, 1990, xii, 407 pp, US $ 103.00 

This book highlights the ephemeral nature of science and 
rapid progress that is being made in understanding adap- 
tation and acclimation to environmental stresses. Ten 
years ago, this reviewer and P. J. Kramer edited a book, 
produced by the same publisher, entitled Adaptation of 
Plants to Water and High Temperature Stress. None of the 
authors in that book is an author in the present volume 
and only six of the papers published in the earlier book 
are cited in the book under review. In part this is because 
the scope of the present book is broader than its predeces- 
sor which was restricted to water and high temperature 
stresses. The authors in the present volume also address 
anthropogenic stresses such as ozone and heavy metals, 
and include cold temperature stress as well as water and 
high temperature stress. 

The book has 17 chapters contributed by 30 authors, 
principally located in North America, but including a 
few from Germany, Sweden and the United Kingdom. 
The book commences with two chapters by Amthor and 
McCree and by Setter on the carbon balance and assimi- 
late partitioning, respectively, of plants subjected to 
stress. While photosynthesis plays a key role in both of 
these processes, the role of stomata is not discussed until 
Chapter 11 in which Mansfield and Atkinson review the 
responses of stomata to water deficits with particular 
reference to the role of abscisic acid. The third chapter 
deals with the mechanisms that plants utilise to cope with 
desiccation, a state that the author, Leapold, points out 
is only transient. The fourth chapter by Winston discuss- 
es the formation of free radicales in cells in response to a 
range of stresses, but it is hard for this reviewer to see 
how this chapter and the ones by Hassan and Scandalios 
(Chapter 8) on superoxide dismutases and Smith and 
co-workers (Chapter 9) on glutathione as an antioxidant 
are relevant to adaptation and acclimation mechanisms 
in plants. Young and Britton in Chapter 5 review the 
role of carotenoids in light harvesting and in the protec- 
tion of chloroplasts against photooxidation. However, 
carotenoids are also subject to photooxidation leading to 
photoinhibition and the destruction of chlorophyll and 
carotenoids when other stresses, such as water stress, oc- 
cur at high light intensities. Since carotenoids are thought 
to be precursors of abscisic acid, the chapter by Morgan 
(Chapter 6) on the role of biotic stresses on plant hor- 
mones is a logical subsequent chapter. Morgan emphasis- 
es the role of abscisic acid and cytokinins, suggesting that 
abscisic acid increases as a result of almost all stresses, 
whereas cytokinins decrease in those stresses in which 
their responses have been studied. Morgan also briefly 
discusses auxins and giberellins, but leaves the role of 
ethylene to Wang and co-workers to review in the follow- 
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ing chapter. Interestingly, these authors discuss the role 
of ethylene as an air pollutant as well as its production by 
the plant in response to chilling, waterlogging, water 
shortage and chemical stresses. 

The book then moves on to review the role of super- 
oxide dismutases and glutathione under stress condi- 
tions, as mentioned above, before Flores (Chapter 10) 
points out that polyamines, notably putrescine, increase 
significantly in response to a series of stresses. Although 
polyamines are thought to be involved in the regulation 
of cell division and cell differentiation, Flores concludes 
that “‘at this point we cannot conclude if this (increase in 
putrescine) is part of a protective response or part of the 
syndrome of toxicity observed under such conditions”. 
Chapters 12 and 13 by Hallgren and Oquist and by 
Burke, respectively, deal with the adaptation and accli- 
mation of plants to low and high temperature stress, re- 
spectively, with Burke introducing readers to the concept 
of thermal kinetic windows. This is the temperature range 
over which enzymatic activity is maximal; while this 
varies with species, plant enzymes do not appear to accli- 
mate to changes in temperature. Burke reviews the influ- 
ence of high temperatures on photosynthesis and carbo- 
hydrate metabolism, but does not discuss the role of heat 
shock proteins which are comprehensively reviewed by 
Vierling in Chapter 16. Chevone etal. in Chapter 14 
point out that water deficits reduce the influence of dam- 
age by the air pollutant, ozone, and that increasing the 
adaptation of plants to water deficits can make the plants 
more susceptible to ozone. This is a good example of the 
interaction between stresses that must always be taken 
into account when studying environmental stresses. 
Chapters by Cumming and Taylor (Chapter 15) and Stef- 
fens (Chapter 17) on heavy metal tolerance and the role 
of photochelatins as metal-binding polypeptides pro- 
duced in response to heavy metal stress conclude the 
book. 

The book is well produced and only one typographical 
error was detected by this reviewer. It has a useful subject 
index, but no author index. As may be clear from the 
comments above, this reviewer had difficulty with the 
order in which the editors placed the chapters. Some 
overlap and redundancy between chapters has occurred, 
as may be expected in a multi-author volume, but this is 
not excessive. The book will be of greatest benefit to 
graduate students and researchers working in stress phys- 
iology. A recurring theme by almost all authors is the lack 
of knowledge in the field under review and the need for 
further research. For example, Vierling (Chapter 16) con- 
cludes that “‘the exact function of any heat shock protein 
is not yet known” and Steffens (Chapter 17) points out 
that “‘the mechanisms of plant adaptation to metal stress- 
es are complex, and a great deal of biochemical, physio- 
logical, and molecular work needs to be done before we 
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reach a complete understanding of this system”. Despite 
the efforts by the editors to develop a series of papers that 
enable general principles on the adaptation and acclima- 
tion to be developed, this reviewer was disappointed 
that no such general principles were identified. Morgan 
(Chapter 6) comes the closest to identifying a role for 
hormones in acclimation of plants to a range of stresses, 
but even he concludes that “‘it is probably safe to say that 
there are few cases where a regulatory role for a hormone 
in a particular stress has been unquestionably estab- 
lished”’. 

Is this a book that will interest readers of Irrigation 
Science? There is little that deals directly with irrigation. 
Burke in Chapter 13 demonstrates that irrigation helps to 
maintain leaf temperatures within the thermal kinetic 
window during periods of high ambient temperature. 
Clearly, the chapters dealing with water deficits are of 
interest to those involved with irrigation practice in re- 
gions where water is in limited supply or expensive, and 
the role of waterlogging as a stress must be of concern to 
students and practitioners of irrigation. However, the 
chapters dealing with heavy metals, air pollutants (except 
where they interact with water shortage or excess), and 
photooxidants are likely to be of only marginal interest. 
It is a worthwhile record of the current (1990) state of 
knowledge in aspects of stress physiology and as such it 
is worthy of a place in your library, but I cannot recom- 
mend that irrigation scientists should have a personal 
copy on their desk. 


Neil C. Turner, CSIRO Division of Plant Industry, Lab- 
oratory for Rural Research, Private Bag, P.O., Wembley, 
W.A. 6014, Australia 


United Kingdom Register of Research and Irrigation, 
Drainage and Flood Control. International Commission 
on Irrigation and Drainage, British Section. Hydraulics 


Research, Wallingford, 1992, 93 pp 


The initiative to compile the United Kingdom Register 
came from the British National Committee of the Inter- 
national Commission on Irrigation and Drainage 
(ICID), and the disciplines covered reflect, to a large 
extent, the interest of ICID members. The reception given 
to the first (1986), second (1988) and third (1990) issues 
of the Register has shown that it is an accepted reference 
work within the profession. In spite of the availability of 
technical journals and a number of national and interna- 
tional conferences, practicing planners, agriculturalists 
and engineers are not always aware of the research being 
undertaken and are therefore not able to fully utilise the 
results. The Register helps to meet the need for dissemi- 
nation of such information. 

The fourth (1992) edition sets out relevant UK re- 
search during the two years period January 1990 to June 
1992. The 196 research entries (compared to 228 entries 
in the third edition) have been organized under three 
headings which are sub-divided into eleven subject areas: 
1. Irrigation — hydrology, agriculture, engineering, man- 
agement. 


2. Drainage and Flood Control — agriculture, river engi- 
neering, policy and management. 

3. Related Topics — monitoring and evaluation, sociology 
and economics, health, environment. 

Each entry comprises project details, publications, 
sponsoring agency and names and addresses of re- 
searchers to contact for further information. Each entry 
is indexed under subject, researcher, location and institu- 
tion. A section of institutional entries from 13 principal 
institutions conducting or sponsoring research in these 
areas is also included describing the organization and the 
type of work undertaken. 

As with previous editions, the Overseas Development 
Administration, which funds a large number of projects, 
has helped to finance the publication and distribution of 
this edition. The Register is compiled and published by 
the Overseas Development Unit of Hydraulics Research, 
Wallingford, and copies are available from: Research 
Register Coordinator, Overseas Development Unit, HR 
Walingford, Howbery Park, Wallingford, Oxon OX10 
8BA, UK. 


Amir Kassam, 88 Gunnersbury Avenue, Ealing, London 
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Shalhevet J., Liu Changming, Xu Yuexian (eds) Water use 
efficiency in agriculture. Priel Publishers, Rehovot, 1992, 
297 pp. 


This volume is the proceedings of a binational workshop 
held in Bijing, China during April 1991. The proeedings 
consist of papers by seven Israeli and nineteen Chinese 
scientists. The primary topics concentrated on the envi- 
ronmental aspects of irrigation including salinity and cli- 
mate, crop response to water, and irrigation technology. 
The objectives of the workshop were to: (1) address the 
importance of water conservation in China, (2) identify 
methods that have been successful for improving water 
use efficiency in both China and Israel, and (3) review 
Israel’s experiences in alleviating its own water deficien- 
cies. 

The seven Israeli papers summarized irrigated agricul- 
ture in Israel today. One paper described how Israel’s 
agriculture has been transformed from ‘“‘an extensive, 
largely self-dependent cropping system based on rainfall 
into an intensive, open system relying on fossil fuel-based 
inputs of water and agrochemicals to produce high-value 
horticultural crops for export”. During this transition 
Israeli crop yields doubled per unit of applied water and 
water use per unit of irrigated land decreased 30%. As 
impressive as these achievements are, the author sur- 
mised that the current water use efficiency has to double 
by the end of the century for Israel to remain competitive 
on a global basis. This goal for water use efficiency is to 
be achieved while eliminating the negative consequences 
of intensive irrigated agriculture on the quality of water 
resources. 

The other Israeli papers dealt with specific issues that 
must be addressed to achieve this goal. One paper de- 





scribed strategies to manage saline soils and waters. The 
potential use of sewage effluent to supplement water re- 
sources was assessed with the major concerns described 
being bacterial contamination, heavy metals, salinity, 
sediment, organic substances, and nutrients. An excellent 
presentation on uniform irrigation applications suggest- 
ed that uniformity be treated as a control variable rather 
than as a predetermined value. Although complete irriga- 
tion uniformity is unattainable and generally not eco- 
nomically desirable, low uniformity was shown to lead to 
significant yield loss and have serious environmental con- 
sequences. Models describing sea water intrusion, a ma- 
jor issue along any coast, were presented along with sev- 
eral practical applications. A relatively new technique to 
measure the flow of sap up the stems of plants has been 
developed to estimate transpiration. The ratio of transpi- 
ration, measured as sap flow, and potential transpiration 
was found to be 0.3 for citrus, 0.5 for pine, and 0.8 for 
both corn and cotton. 

The majority of the Chinese papers addressed efficient 
water resource utilization on the North China Plain. This 
Plain constitutes one-third of China’s flatland and ac- 
counts for over 1.8 million hectares of cultivated land. 
Since 1949, 7 100 reservoirs have been built to control the 
runoff of rivers transversing the North China Plains and 
1.4 million tube wells have been installed to pump 60% of 
the groundwater recharge. The environmental issues 
raised because of excess use of water resources include 
soil salination, groundwater depletion, sea water intru- 
sion, and sediment. Modelling efforts described for the 
North China Plains included the conjunctive use of sur- 
face and groundwater to increase water use efficiency and 
control salinity, the most effective cropping patterns and 
rotations to match water demand with water resources, 
and water management strategies to improve crop water 
use efficiency. 

Wheat and cotton are the dominate crops on the 
Plains so that several papers dealt with improving water 
management of these crops. The potential for saving wa- 
ter at the initiation of peach fruit enlargement by with- 
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holding irrigation and then irrigating fully as fruit devel- 
opment proceeded was reported to increase production. 
Another specific topic under study in China was the use 
of mulches and plastic films to conserve water and main- 
tain production. Similar to studies elsewhere, wheat pro- 
duction was found to increase linearly with applied water 
until 200 mm of irrigation had been applied in three irri- 
gations. With respect to total crop water use (irrigation, 
rainfall, and stored soil water), wheat yields increased 
linearly up to 500 mm. One report indicated that drip 
irrigation systems designed specially for Chinese condi- 
tions were successful for grain, cotton, and oil crops. 

One of the major contributions of the proceedings was 
the documentation of environmental concerns on the 
North China Plains. Reports indicated that about 2 mil- 
lion hectares of salt-affected soils now exist in the Plains; 
the rate of water table decline in areas like Cangzhou are 
as much as 5 m per year with more than 7 000 km? affect- 
ed by declining groundwater; the density of tube wells 
near the sea coast is more than 19 wells per km? with the 
water table now being 5 to 10m below sea level; and 
along the lower reaches of the Yellow River more than 
200 million Mg of sediment are deposited annually. Due 
to one or more of these factors, the irrigated area has 
decreased recently by more than 140 000 ha, the average 
water table elevation during the last 20 years has declined 
6m, and the amount of stream flow reaching the sea is 
now only 6% of the flow 40 years ago. 

The reports in this volume serve to highlight the ur- 
gent need for enhanced water use efficiency in both China 
and Israel. This issue is not unique to these two countries. 
In every country in the world practicing irrigation, these 
same issues are paramount. I strongly recommend the 
proceedings of this workshop to all concerned with food 
production, water resource protection, and enhancement 
of the environment. 


G. J. Hoffman, Department of Biological Systems Engi- 
neering, University to Nebraska, Lincoln, NE 68583-0726, 
USA. 
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Saturated Flow 


and Soil Structure 


A Review of the Subject and 
Laboratory Experiments on the 
Basic Relationship 


1993. XI, 235 pp. 141 figs. 4 tabs. Hardcover 
DM 148,- ISBN 3-540-55791-1 


Because soil and water pollution have a serious 
impact on the environment, soil scientists and 
hydrologists need fundamental help for the esti- 
mation of the consequences. 

The experiments described in this volume deal 
with the quantification of the morphology of 
interaggregate voids and of the flow through 
such voids as well as around impermeable 
inclusions. The diagrams given in the appendix 
can be used as references for such measure- 
ments. 
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Contaminants in 
Terrestrial Environments 


1993. XVI, 439 pp. 75 figs. 3 fold-outs. 
Hardcover DM 228,- ISBN 3-540-55277-4 


O. Franzle 


Defining ecology as a system-theory oriented 
synthesis of both earth and life sciences, this 
book aims at a novel coherent understanding of 
chemical impact on the lower atmosphere and 
characteristic types of terrestrial ecosystems. To 
this end comprehensive flux-analytical and hier- 
archical modelling approaches were developed 
and consulted, which include a thorough con- 


sideration of the specific physical and geograph- 


ic boundary conditions of the processes 
involved. 


Volume 10 
CJ.M.G. Pieri 
Fertility of Soils 


A Future for Farming 

in the West African Savannah 
Translated from the French by P. Gething 
1992. XVIII, 348 pp. 87 figs. Hardcover 

DM 257,- ISBN 3-540-53283-8 

The results of 30 years of investigation in the 
West African savannah are re-interpreted from 


the point of view of modern agronomical and 
environmental sciences. 


Prices are subject to change without notice. In EC countries the local VAT is effective. 


For information on prices in Austrian schillings and Swiss francs please consult the German book directory 
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Volume 11 
R. Young, A. Young 


Sandstone Landforms 


1992. VIII, 163 pp. 63 figs. 2 tabs. Hardcover 
DM 154,- ISBN 3-540-53946-8 


As very few studies on geomorphological 
aspects of sandstones exist, this book is impor- 
tant in that it demonstrates the significance of 
landforming processes from geological, geo- 
chemical and dynamic points of view. 


Volume 12 
G. Gro& 


Numerical Simulation 
of Canopy Flows 


1993. VIII, 167 pp. 127 figs. Hardcover 
DM 158,- ISBN 3-540-52520-3 


Starting with the description of meteorological 
variables in forest canopies and its parameter 
variations, a numerical three-dimensional 
model is developed. Its applicability is demon- 
strated, first, by wind sheltering effects of hedg- 
es and, second, by the effects of deforestation 
on local climate in complex terrain. 
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